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Abstract 



At present, the W boson is investigated experimentally at LEP2 via its pair production. In 
order to achieve precise theoretical predictions for the measurement of the W-boson mass 
and the non-abelian triple-gauge-boson couplings, the inclusion of radiative corrections 
is required. Since the W bosons decay very rapidly into light fermion pairs, the actual 
processes under investigation are e + e~ — > 4 fermions. 

The full 0(a) corrections to these processes are not available at present. Since the 
main contributions originate from diagrams with two resonant W-boson propagators, an 
expansion of the amplitude around the poles of the two resonant W bosons is a reasonable 
approach, which is also gauge-invariant. In the double-pole approximation, the contribu- 
tions are classified in factorizable and non-factorizable corrections. The amplitudes of the 
factorizable corrections are composed of those for the on-shell W-pair production and the 
on-shell W-decays multiplied by the two propagators of the resonant W bosons. All other 
corrections are called non-factorizable, because they do not factorize into a simple product 
of W-pair production and W decays. 

As a first step of this work, the non-factorizable corrections of the processes e + e~ — > 
4 fermions are calculated in double-pole approximation. The non-factorizable corrections 
are implemented into an existing Monte Carlo program and various distributions are stud- 
ied. It turns out that the non-factorizable corrections are negligible with respect to the 
experimental accuracy of LEP2; however, they should become relevant for a future linear 
collider with higher luminosity. 

A further building block of the radiative corrections to four-fermion production are the 
bremsstrahlung processes e + e~ — > 4 fermions + 7. These processes are of physics interest in 
their own right. For instance, the radiative processes can be used to obtain information on 
the quartic-gauge-boson couplings 77WW, 7ZWW, and 77ZZ, which are part of the tree- 
level amplitude. The tree-level helicity amplitudes for the processes e + e~ — > 4 fermions and 
e + e~ — > 4 fermions + 7 for all possible final-state fermions are calculated. A multi-channel 
Monte Carlo program for both classes of processes is constructed. This is a particularly 
difficult task owing to the very complex peaking structure of the differential cross section. 

In order to include radiative corrections into the Monte Carlo program, the infrared and 
collinear singularities must be extracted from the bremsstrahlung process. This is done 
by applying the dipole-subtraction method to four-fermion production. This subtraction 
method has already been worked out in massless QCD for dimensional regularization. In 
the Electroweak Standard Model it is more convenient to regularize the amplitude with 
an infinitesimal photon mass and small fermion masses. Hence, the subtraction method 
is reformulated for mass regularization. 



Ill 



Finally, the doubly-resonant virtual corrections are implemented into the Monte Carlo 
program for the processes e + e~ — > 4 fermions(+7). Therefore, the results of the non- 
factorizable corrections and the already existing results of the on-shell W-pair production 
and on-shell W-decay are used. For the real corrections the complete bremsstrahlung 
process e + e~ — > 4fermions + 7 is taken into account. All results are combined in a four- 
fermion generator, which is the first Monte Carlo generator including the complete 0(a) 
corrections to the processes e + e~ — ► W + W~ — > 4 fermions in double-pole approximation. 
This generator is used to produce numerical results for the total cross section, angular, 
and invariant-mass distributions. 
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Zusammenfassung 



Am Beschleuniger LEP2 wird zur Zeit die W-Paar-Produktion experimentell untersucht. 
Um genaue theoretische Vorhersagen fur die Bestimmung der W-Boson-Masse und die 
Untersuchung der Drei-Eichboson-Kopplungen zu erhalten, miissen Strahlungskorrekturen 
beriicksichtigt werden. Weil die W-Bosonen sehr schnell in je zwei leichte Fermionen 
zerfallen, werden in Wirklichkeit die Prozesse e + e~ — > 4 Fermionen untersucht. 

Die vollstandigen Strahlungskorrekturen in 0(a) sind zu diesen Prozessen noch nicht 
bekannt. Der Hauptbeitrag stammt von Diagrammen mit zwei resonanten W-Boson- 
Propagatoren. Deshalb besteht ein naheliegender und auch eichinvarianter Ansatz darin, 
die Amplitude um die Pole der beiden W-Bosonen zu entwickeln. In dieser Doppelpol- 
naherung konnen die Strahlungskorrekturen in faktorisierbare und nicht-faktorisierbare 
Korrekturen klassifiziert werden. Die Amplitude der faktorisierbaren Korrekturen setzt 
sich aus den Amplituden der W-Paar-Produktion, den beiden Amplituden der W-Boson- 
Zerfalle und den zwei Propagatoren der resonanten W-Bosonen zusammen. Alle iibrigen 
Korrekturen werden mit nicht-faktorisierbar bezeichnet, weil sie nicht aus einem einfachen 
Produkt von Beitragen zur Produktion und zu den Zerfallen geschrieben werden konnen. 

In einem ersten Schritt wurden in der vorliegenden Arbeit die nicht-faktorisierbaren 
Korrekturen zu den Prozessen e + e~ — > 4 Fermionen in Doppelpolnaherung berechnet. Sie 
wurden in ein existierendes Monte-Carlo-Programm eingebaut, und verschiedene Verteil- 
ungen wurden studiert. Dabei zeigte sich, dass diese Korrekturen vernachlassigbar gegen- 
iiber dem experimentellen Fehler von LEP2 sind. Jedoch werden sie fur einen zukunftigen 
Linearbeschleuniger mit hoherer Luminositat voraussichtlich wichtig. 

Ein weiterer notwendiger Bestandteil der Strahlungskorrekturen zur Vier-Fermion-Pro- 
duktion sind die Bremsstrahlungsprozesse e + e~ — ► 4 Fermionen + 7. Mit diesen Prozessen 
konnen auch die Vier-Eichboson-Kopplungen 77WW, 7ZWW und 77ZZ studiert wer- 
den, die auf Born-Niveau enthalten sind. Die Helizitats-Amplituden zu den Prozessen 
e + e~ — > 4 Fermionen und e + e~ — > 4 Fermionen+7 fur alle Endzustande wurden berech- 
net, und ein Monte-Carlo-Programm fur beide Klassen von Prozessen geschrieben. Die 
Schwierigkeit lag dabei in dem sehr komplexen und stark variierenden Verhalten des dif- 
ferentiellen Wirkungsquerschnittes. 

Um Strahlungskorrekturen mit Hilfe eines Monte-Carlo-Programms zu berechnen, miis- 
sen die infraroten und kollinearen Singularitaten vom Bremsstrahlungsprozess abgespalten 
werden. Dazu wurde die Dipol-Subtraktionsmethode auf die Vier-Fermion-Produktion 
angewandt. Diese Subtraktionsmethode existierte in der Literatur fiir masselose QCD 
und dimensionale Regular isierung. Fiir das Elektroschwache Standardmodell werden nor- 
malerweise die Singularitaten mit einer infinitesimalen Photonmasse und kleinen Fermion- 



V 



massen regularisiert. Daher wurde die Subtraktionsmethode fiir die Massenregularisierung 
umgeschrieben. 

Alle doppeltresonanten virtuellen Korrekturen wurden in das vorher erwahnte Monte- 
Carlo-Programm fiir die Prozesse e + e~ — > 4 Fermionen(+7) eingebaut. Dazu wurden die 
nicht-faktorisierbaren Korrekturen und die schon existierenden Ergebnisse fiir die W-Paar- 
Produktion und W-Zerfalle verwendet. Fiir die reellen Korrekturen wurde der komplette 
Bremsstrahlungsprozess e + e~ — > 4 Fermionen+7 beriicksichtigt. Dieses Programm ist der 
erste Monte-Carlo-Generator, der alle Strahlungskorrekturen in 0(a) zur Vier-Fermion- 
Produktion in Doppelpolnaherung beinhaltet. Mit ihm wurden numerische Ergebnisse fiir 
den totalen Wirkungsquerschnitt, Winkelverteilungen und invariante-Massen-Verteilungen 
erzeugt. 
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Chapter 1 
Introduction 



The Glashow-Salam- Weinberg model |1|, known as the Electroweak Standard Model 
(SM), is very successful in describing electroweak phenomena. Since the SM is a spon- 
taneously broken gauge theory, it is renormalizable [^] and hence observables can be, in 
principle, calculated to any finite order in perturbation theory. 

An important feature of the SM lies in the appearance of elementary gauge-boson-self- 
interactions resulting from the non-abelian structure of the gauge group. The 7WW and 
ZWW vertices can be studied in detail at the e + e~ collider LEP2 ||. Beside the investi- 
gation of the triple-gauge-boson-couplings, LEP2 also allows for a precise determination 
of the W-boson mass. Two methods are used @, ||: the measurement of the total cross 
section near threshold and the reconstruction method, where the Breit-Wigner resonance 
shape is reconstructed from the decay products of the W bosons. 

LEP2 is operating above the W-pair production threshold and produces about 10 4 
W pairs. Hence, the typical experimental accuracy is of the order of one to a few per 
cent. The accuracy of the W-boson mass measurement is expected to be < 50MeV at 
LEP2 Q and about 15MeV for a future linear collider ||. This experimental accuracy 
should be matched or better exceeded by the precision of the theoretical predictions. Since 
W bosons decay very rapidly into light fermions, the actual reaction under investigation 
is e + e~ — > W + W~ — > 4 fermions. 




Figure 1.1: Diagrams with two resonant W-boson propagators contributing to e + e — > 4/ 
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In the LEP2 energy region, the lowest-order cross section is dominated by the diagrams 



that involve two resonant W bosons, as shown in Fig. |1.1| . All other lowest-order diagrams 
are typically suppressed by a factor T w /M w ps 2.5%, but may be enhanced in certain 
phase-space regions. Since all these contributions are required at the one-per-cent level, 
the complete lowest-order matrix element has to be taken into account. 

Furthermore, the implementation of the finite width of unstable particles, such as 
the W bosons, has to be done properly. A finite width is necessary in the phase-space 
region where the unstable particle becomes resonant, i.e. nearly on shell, otherwise the 
cross section has a non-integrable unphysical singularity. The finite width is naturally 
introduced via Dyson summation, where the width arises from the imaginary parts of the 
resummed self-energy diagrams: 



oo 



k 2 ~ M 2 to 



-Sy(fc 2 ) 

k 2 - Ml 



1 (l.o.r 



k 2 - M 2 + Z v (k 2 



where M symbolizes the bare mass of the unstable particle, and £y denotes the one- 
particle- irreducible self-energy. 

However, since only a part of higher-order corrections are included in this way, the 
whole result is gauge-dependent, and wrong results can be obtained in certain phase-space 
regions |7|, p|, |9j jTOfl . The reason is that Ward identities are violated and, hence, gauge 
breaking terms can be amplified in the presence of small scales, or unitarity cancellations 
do not take place properly. Ward identities are, in particular, crucial for processes with 
nearly on-shell virtual photons or for the production of longitudinal polarized gauge bosons 
at high energies. 

For the calculation of the tree-level processes e + e~ — > 4/ and e + e~ — > 4/7, three 
schemes for the implementation of the finite W-boson width are compared in Chapter |2|: 
the constant-width scheme, where the imaginary part of the self-energy is replaced by a 
constant, the running-width scheme with the naive running T(k 2 ) = Ty8(k 2 )k 2 /My, and 
the complex-mass scheme, where the boson masses are replaced by complex masses in all 
couplings and propagators and, in particular, in the definition of the weak mixing angle. 
In general, the first two schemes violate Ward identities || [Kj], while the third scheme 
fulfils all Ward identities. 

A gauge- invariant approach for the introduction of the finite width is the pole expansion 
Tl, [T2| outlined in the following. The complete matrix element for a process, with an 



unstable particle in the intermediate state, can be written as 

where the residue is denoted by r, and the mass of the unstable particle by M. The 
symbol n summarizes all terms that are regular at k 2 = M 2 . Both the location of the 
complex pole, which corresponds to the mass of the unstable particle, and the residue are 
gauge-invariant quantities. For stable particles the mass lies on the real axis, for unstable 
particles in the complex plane. 

In perturbation theory, the complex mass M is determined by the location of the pole 
after Dyson summation ( |1.0.1|) : 

M 2 = M 2 - S V (M 2 ) = M v - iM v r v , (1.0.3) 
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where Ty and My denote the finite width and the real mass of the unstable particle, 
respectively. After expanding the self-energy about M 2 , the inverse propagator reads 

k 2 - Ml + S V (A; 2 ) = [1 + £ V (M 2 )P 2 — M 2 ) + O ((k 2 - M 2 ) 2 ) , 
and the matrix element can be rewritten in the form 



k 2 - M 2 + E y (k 2 ) K ' ' [1 + £ V (M 2 )](P — M 2 ) 
R(k 2 ,9) R{M 2 ,6) 



+ 



k 2 - Ml + E y (k 2 ) [1 + £' V (M 2 )](£; 2 -M 2 ) 



+ N(k*,0), (1.0.4) 



where N includes the non-resonant diagrams, and 6 summarizes all kinematic variables, 
except for the invariant mass k 2 of the unstable particle. Since the complete amplitude 
is gauge- invariant, the single terms of the Laurent expansion about the squares of the 
complex mass are also gauge-invariant. The first term on the right-hand side of ( |1.0.4| ) 
corresponds to the leading term in a Laurent expansion about k 2 = M 2 , and dominates the 
cross section in the resonance region. The remaining terms are finite in the limit k 2 — > M 2 . 
Therefore, a reasonable and gauge-invariant approximation is to neglect the non-leading 
terms and to keep only the resonant term. This simplifies the calculation considerably, 
since all obviously non-resonant diagrams can be left out from the calculation from the 
beginning. 

In Chapter [| and |], the pole scheme is applied to the radiative corrections of four- 
fermion production. Radiative corrections are, in general, required in the theoretical 
predictions for four-fermion production, in order to match the accuracy of LEP2 of about 
one per cent. The full 0(a) calculation involves 10 3 -10 4 diagrams and is, therefore, ex- 
tremely complicated. Since the non-doubly-resonant radiative corrections are of the order 
of aT*w\a(. . )/{ , kMyj) ~ 0.1%, the restriction to the doubly-resonant corrections is a 
reasonable approach. 

After the introduction of the finite W-boson width, the matrix element reads 



M - eg) R + (klkle) R_{k 2 +1 k 2 _^, 

M ~ {k\ -M 2 )(k 2 -M 2 ) + 4-M 2 + k 2 _-M 2 



where k± are the momenta of the resonant virtual W bosons and the factor [1 + £ V (M 2 )] 
is already included in the definition of R+-, R+, and -R-Q In double-pole approximation, 
the matrix element is expanded about the squares of the W-boson masses and all non- 
doubly-resonant terms are neglected: 

R+ -(M\M\e) 

{k% — M 2 )(fcl — M 2 ) [ 1 

The double-pole approximation simplifies the calculation considerably, since only dia- 
grams with two resonant W-boson propagators have to be calculated. Furthermore, a 



1 We use the the on-shell renormalization scheme of Ref. [jl3), where Re{Ey(My)} = 0. In this renor- 
malization scheme, the deviation of the self-energy at the complex pole Ey(M 2 ) yields only 0(a 2 ) cor- 
rections. Hence, the factor [1 + E v (Af 2 )] is neglected in the following chapters. 
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large part of the doubly-resonant radiative corrections are included in the factorizable cor- 
rections. The factorizable corrections are composed of the on-shell W-pair production, the 
two resonant W-boson propagators, and the two on-shell W-boson decays. The remaining 
radiative corrections are called non-factorizable and are explicitly calculated and discussed 
in Chapter |3|. 

The double-pole approximation is a good approximation if the doubly-resonant con- 
tributions dominate the cross section. If these are suppressed, as close to the W-pair 
production threshold, the other contributions become important. This is also true if non- 
doubly-resonant terms are enhanced, as e.g. by nearly on-shell photons. These terms can 
be suppressed by applying appropriate cuts on the phase space. 

Note that the pole expansion only works if the on-shell limit exists. Since the non- 
factorizable corrections, given in Section |3.3| , involve on-shell-divergent terms, like ha.(k± — 
M 2 ), the pole expansion of these corrections has not such a simple form as in ( |1.0.6|) . 
Thus, the non-factorizing corrections are calculated for off-shell W bosons, while the limit 
k± — > M 2 is performed whenever possible. 

With the results of Chapter [2] and |3], a Monte Carlo generator is constructed in 
Chapter [| that includes all doubly-resonant 0(a) radiative corrections to four-fermion pro- 
duction. More precise, this generator includes the complete tree- level matrix element, the 
virtual corrections in double-pole approximation, and the complete bremsstrahlung pro- 
cess. The cancellation of the soft and collinear singularities is achieved within the subtrac- 
tion method as discussed in Section |4.2| . Although the real non-factorizable corrections are 
calculated in Chapter ^| in double-pole approximation, the complete bremsstrahlung pro- 
cess is taken into account in Chapter [|. In this way, the problem of overlapping resonances 
is avoided (see Section [4.1.3|) . The Monte Carlo program of Chapter ^ is the first generator 
that includes the complete 0(a) corrections to the processes e + e~ — > W + W~ — > 4fermions 
in double-pole approximation. 
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Chapter 2 

Tree-level processes e+e — — ► 4/( +7) 



While the most important process at LEP2 for the studies of the gauge sector in the 
Electroweak Standard Model is certainly e + e~ — > W + W~ — > 4/, many other reactions 
have now become accessible. Besides the 4-fermion-production processes, including single 
W-boson production, single Z-boson production, or Z-boson-pair production, LEP2 and 
especially a future linear collider allow us to investigate another class of processes, namely 
e + e~ — ► 4/7. 

The physical interest in the processes e + e~ — > 4/7 is twofold. First of all, they are an 
important building block for the radiative corrections to e + e~ — > 4/, and their effect must 
be taken into account in order to get precise predictions for the observables that are used for 
the measurement of the W-boson mass and the triple-gauge-boson couplings. On the other 
hand, those processes themselves involve interesting physics. They include, in particular, 
triple-gauge-boson-production processes such as W + W~7, ZZ7, or Z77 production and 
can therefore be used to obtain information on the quartic gauge-boson couplings 77 WW, 
'jZWW, and 77ZZ. While only a few events of this kind are expected at LEP2, these 
studies can be performed in more detail at future linear e + e~ colliders |L4] . 

Some results for e + e~ — > 4/7 with an observable photon already exist in the literature. 



In Refs. p5| , |16|| the contributions to the matrix elements involving two resonant W bosons 
have been calculated and implemented into a Monte Carlo generator. This generator has 
been extended to include collinear bremsstrahlung |T7[ and used to discuss the effect of 



hard photons at LEP2 [18]. The complete cross section for the process e + e~ — > ude~z/ e 7 
has been discussed in Ref. [|B|. In Ref. P0fl , the complete matrix elements for the processes 
e + e~ — > 4/7 have been calculated using an iterative numerical algorithm without referring 
to Feynman diagrams. We are, however, interested in explicit analytical results on the 
amplitudes for various reasons. In particular, we want to have full control over the imple- 
mentation of the finite width of the virtual vector bosons and to select single diagrams, 
such as the doubly-resonant ones. No results for e + e~ — > 4/7 with e + e~ pairs in the final 



state have been published in the past. The results of this section are published in Ref. [21 



and agree very well with the recent calculations of Ref. |22fl , where finite- mass effects due 
to nearly collinear photon emission are discussed for the process e + e" — > ndji^u^y. 

In order to perform the calculation as efficient as possible we have reduced all processes 
to a small number of generic contributions. For e + e~ — > 4/, the calculation is similar 
to the one in Ref. [J23|, and the generic contributions correspond to individual Feynman 
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diagrams. In the case of e + e~ — > 4/7 we have combined groups of diagrams in such a 
way that the resulting generic contributions can be classified in the same way as those for 
e + e~ — > 4f. As a consequence, the generic contributions are individually gauge-invariant 
with respect to the external photon. The number and the complexity of diagrams in the 
generic contributions for e + e~ — > 4/7 has been reduced by using a non-linear gauge-fixing 
condition for the W-boson field [24]. In this way, many cancellations between diagrams are 
avoided, without any further algebraic manipulations. Finally, for the helicity amplitudes 
corresponding to the generic contributions concise results have been obtained by using the 
Weyl-van der Waerden formalism (see Ref. E| and references therein). 

After the matrix elements have been calculated, the finite widths of the resonant par- 
ticles have to be introduced. We have done this in different ways and compared the 
different treatments for e + e~ — > 4/ and e + e~ — > 4/7. In particular, we have discussed a 
"complex-mass scheme", which preserves all Ward identities and is still rather simple to 
apply. 

The matrix elements to e + e" — > 4/ and e + e~ — > 4/7 exhibit a complex peaking be- 
haviour owing to propagators of massless particles and Breit-Wigner resonances, so that 
the integration over the 8- and 11-dimensional phase spaces, respectively, is not straightfor- 
ward. In order to obtain numerically stable results, we adopt the multi-channel integration 
method [fj^, |26| and reduce the Monte Carlo error by the adaptive weight optimization 
procedure described in Ref. [^7| . In the multi-channel approach, we define a suitable map- 
ping of random numbers into phase-space variables for each arising propagator structure. 
These variables are generated according to distributions that approximate this specific 
peaking behaviour of the integrand. For e + e~ — > 4/ and e + e~ — ■> 4/7 we identify up to 
128 and 928 channels, respectively, which necessitates an efficient and generic procedure 
for the phase-space generation. 



2.1 Analytical results 



2.1.1 Notation and conventions 

We consider reactions of the types 

e + (p+,o-+) +e~(p_,er_) fi(ki,<Ti) + f 2(^2, 02) +/3O3, 03) + f 4 (k 4 ,a 4 ), (2.1.1) 
e+ (p+>°"+) +e"(p_,o-_) -> /1 (fci, 01) + f 2 (k2,a 2 ) + h(h,a z ) + f 4 (k 4 ,a 4 ) + 7(^5, A). 

(2.1.2) 

The arguments label the momenta p±, ki and helicities Oi = ±1/2, A = ±1 of the cor- 
responding particles. We often use only the signs to denote the helicities. The fermion 
masses are neglected everywhere. 

For the Feynman rules we use the conventions of Ref. [0, 58[] . In particular, all fields 
and momenta are incoming. It is convenient to use a non-linear gauge-fixing term [24] of 
the form 

2 



C 



fix 



&>W+ + ie(A^ - — Z")W- 



iM 



- i (d»Z,-M zX ) 2 --(d»A p 



(2.1.3) 
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where and x are the would-be Goldstone bosons of the W and Z fields, respectively. 
With this choice, the (fy^W^A vertices vanish, and the bosonic couplings that are relevant 
for e + e~ — > 4/7 read 




= -iegvww [9u P {k- - k+) M - 2g fiV k v ,p + 2g w k VjU ] , 



V v w- 

with V = A, Z, and the coupling factors 

9aww = 1, 9zww = — — ■ (2.1.5) 

s w 

Note that the gauge-boson propagators have the same simple form as in the 't Hooft- 
Feynman gauge, i.e. they are proportional to the metric tensor g^ v . This gauge choice 
eliminates some diagrams and simplifies others owing to the simpler structure of the 
photon-gauge-boson couplings. 

The vector-boson-fermion-fermion couplings have the usual form 




(2.1.6) 



fj 

where u± = (1 ± 7s)/2. The corresponding coupling factors read 

9am = -Qh 9%^ = ~—Qi + r~^-' dwfifi = 779c - < ^ <T ~' ( 2 - L7 ) 

where Qj and Zj^ = ±1/2 denote the relative charge and the weak isospin of the fermion 
fi, respectively, and // is the weak-isospin partner of /j. The colour factor of a fermion 
is denoted by Nf., i.e. iVf epton = 1 and A^ uark = 3. 



2.1.2 Classification of final states for e + e — ► 4 / 



The final states for e + e 



4/ have already been classified in Refs. ^B], [H5). We 

It is based on 



introduce a classification that is very close to the one of Refs. p9 , [30 
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the production mechanism, i.e. whether the reactions proceed via charged-current (CC), 
or neutral-current (NC) interactions, or via both interaction types. The classification 
can be performed by considering the quantum numbers of the final-state fermion pairs. 
In the following, / and F denote different fermions (/ ^ F) that are neither electrons 
nor electron neutrinos (f,F ^ e~,v e ), and their weak-isospin partners are denoted by 
/' and F', respectively. We find the following 11 classes of processes (in parenthesis the 
corresponding classification of Ref. is given): 



(i) CC reactions: 




(a) e+e" - ff'FF', 


(CC11 family), 


(b) e+e" - u e e+ff', 


(CC20 family), 


(c) e+e" -> ff'e~P e , 


(CC20 family), 


ii) NC reactions: 




(a) e+e- - ffFF, 


(iVC^ family), 


( b ) e+e- -> //./•/. 


(NC4-16 family) 


(c) e+e" -> e~e + ff, 


(iVC^S family), 


(d) e+e" — > e _ e+e _ e+, 


(NC4-36 family) 



(iii) Mixed CC/NC reactions: 

(a) e+e" -> ///'/', family), 

(b) e+e" -> !/ e P e //, (iVC2i family), 

(c) e+e" — > v Q v Q v Q v^ (NC4- 9 family), 

(d) e+e" — > z/ e z/ e e"e+, (mix56 family). 



The radiation of an additional photon does not change this classification. 

2.1.3 Generic diagrams and amplitudes for e + e _ — ► 4 / 

In order to explain and to illustrate our generic approach we first list the results for e+e" — > 
4/. All these processes can be composed from only two generic diagrams, the abelian and 
non-abelian diagrams shown in Fig. |2.1| . All external fermions f a ,...j are assumed to be 
incoming, and the momenta and helicities are denoted by p a ,...j and cr a ,...j, respectively. 
The helicity amplitudes of these diagrams are calculated within the Weyl-van der Waerden 
(WvdW) formalism following the conventions of Ref. []25| (see also references therein) . 

Leptonic and semi-leptonic final states 

We first treat purely leptonic and semi-leptonic final states. In this case, none of the gauge 
bosons in the generic graphs of Fig. [2.1| can be a gluon, and the colour structure trivially 
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a) abelian graph 

fa(Pa, O - a ) 





fb(pb,cr b ) 



fc(Pc, <?c) 

fdiPd, <?d) 

fe(Pe,(Te) 
ff(Pf>°f) 



b) non-abelian graph 



fa(Pa, <7a) 




fb(Pb,<7a) 



fc(Pc, &c) 
fdiPd, (7a) 

fe(Pe,(7e) 
ffiPf^f) 



Figure 2.1: Generic diagrams for e + e — > 4/ 



leads to a global factor N^Nj 3 , which is equal to 1 or 3, after summing the squared 
amplitude over the colour degrees of freedom. The results for the generic amplitudes are 



M VlV2 {Pa,Pb,Pe,Pd,Pe,Pf) 

°b 



-, a S 



4e 5 aa - ab 5 ac - ad 5 ae - af 9v 1 f afg 9 v lf gfb 9 v lf cfd 9v 2 f eff 



x 



Pvi {Pc + Pd)Pv 2 (Pe + Pf) A a a 



(Pb+Pe+Pf)' 



A 2 a ' c ' e (p OJ p b ,p c ,Pd,p e ,Pf), 



(2.1.? 



M vWW C,ad,ae (Pa ,Pb,Pe,Pd,Pe,Pf) 



-4e 5 Ca - ab 5 act+ 5 ad _5 aet+ 5 af _ (Q c - Qd)gvww9vf a f b 9 W f c f d 9 W f e f f 

X Pv(Pa+Pb)Pw(Pc+Pd)Pw(Pe+Pf) ^ (Pa, Pb, Pe, Pd, Pe, P f) , (2.1.9) 



where the vector-boson propagators are abbreviated by 

1 



Pv(p) 



V = A,Z,W,g, 



M A M, 0. 



(2.1.10) 



(The case of the gluon is included for later convenience.) The auxiliary functions J 4^ a ' <Tc ' <Te 
and A^ a are expressed in terms of WvdW spinor products, 

At ++ (Pa,Pb,Pc,Pd,Pe,Pf) = (PaPc}{PbPf}*((PbPd}*(PbPe) + (PdP /)* (PePf)) , 
At + ~(Pa,Pb,Pc,Pd,Pe,Pf) = Aj ++ (p a ,p b ,p c ,p d ,p f ,p e ), 
At~ + (Pa,Pb,Pc,Pd,Pe,Pf) = A^ ++ (Pa,Pb,Pd,Pc,Pe,Pf), 
A 2~~(Pa,Pb,Pc,Pd,Pe,Pf) = A^ ++ (p 

a, Pb, Pd, Pc Pf, Pe), 

A2' aC ' ai (Pa,Pb,Pc,Pd,Pe,P}) = (A+^'- Ud (p a ,p b ,p c ,p d ,p e ,p } ))\ (2.1.11) 

Aj(p a ,Pb,Pc,Pd,Pe,Pf) = (PbPd}*(PbPf}*(PaPb}(PcPe) + (PbPd)* (PdPf)* (PaPe) (PePd) 

+ {PbPf}*{PdPf}*(PaPe}(PePf}, 
As(Pa,Pb,Pc,Pd,Pe,Pf) = Aj(p b ,p a ,p c ,p d ,p e ,p f ). (2.1.12) 
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The spinor products are defined by 



(pq) 



£ AB PAqB = %^/PoCfo 



Or, 



" p cos — sin — - e h 



9 n 



r , 



•a 1 • V 

q cos — Sill — 



2 2 2 2 

where pa, Qa are the associated momentum spinors for the momenta 

— Po(l) si 11 P cos P , sin 9 P sin p , cos 9 P ), 

— Po(l) sin 6^ cos <f) q , sin 6 q sin^, cos# g ). 



(2.1.13) 



(2.1.14) 



Incoming fermions are turned into outgoing ones by crossing, which is performed by 
inverting the corresponding fermion momenta and helicities. If the generic functions are 
called with negative momenta —p, —q, it is understood that only the complex conjugate 
spinor products get the corresponding sign change. We illustrate this by simple examples: 



A(p, q) = (pq) 
B(p, q) = (pq)* 



A(p, -q) = A(-p,q) = A(-p,-q), 
-B(p, -q) = -B(-p, q) = B(-p, -q). 



;2.i.i5) 



We have checked the results for the generic diagrams against those of Ref. |23| and found 
agreement. 

Using the results for the generic diagrams of Fig. |2.1| , the helicity amplitudes for all 
possible processes involving six external fermions can be built up. It is convenient to 



construct first the amplitudes for the process types CC(a) and NC(a) (see Section [2.1.2|) 
in terms of the generic functions ( |2.1.8| ) and ( |2.1.9|) , because these amplitudes are the basic 
subamplitudes of the other channels. The full amplitude for each process type can be built 
up from those subamplitudes by appropriate substitutions and linear combinations. 
We first list the helicity amplitudes for the CC processes: 



(7+ ,cr- ,(T1 ,CT2 ,CT3 ,(T4 

— J vl ww 



(p + ,p_,fa,k 2 ,fa,fa) 

'-^- a2 '-^- ai ( P+ ,P-, -fa, -fa, -fa, -h) 



-r 2^ [ JV Www 

V=j,Z 






-, ~fa, —fa 


-fa, -h) 














1 A / j— fl, — CT2,CT+,(T_, — CT3,- 




fa, 


-k 2 ,p+,p- 


-k 3 , -k 4 ) 




_j_ KA— C3)— <X4,<7+,(7_,— IJ\,- 
~r JVl VW 




fa, 


-fa,p + ,p- 


-fa, -k 2 ) 




_j_ \J—0'l,—0'2t—<73, — <74i<7+ 
"f JVl W y 




ki, 


-k 2 , -k 3 , - 


fa,P+,P-) 




, kyf — C3,— C4, — Cl,— <T2,CT4- 
T JVl W y 




fa, 


—k^, —fa, — 


fa,P+,P-) , 


(2.1.16) 



1 y cr + ,(7_ ,tri ,(T2 ,0-3 ,tT4 

Jvl CCb 



(p+,P-,fa,fa,fa,fa) 

-,CT- 

CCa 



M^-'^'^ip^P-, fa, fa, fa, fa) 



(2.1.17) 



M a c + c:-' ai ' a2 ' a3 ' ai (p + ,p-,fa,k 2 ,fa,fa) 

M^:-^^\p + ,p.,fa,k 2 M,ki) 



-fa,p~, fa, k 2 , -p+, fa) 



(2.1.18) 
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The ones for the NC processes are given by 
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l CCa 

+ M c ^- ana2 ^ a -(-ku ~h, h, -p_, -p +l fe) 

+ -Mcca" CT2 *" ,T+ * <T4 * CTl ' _ "(-^3, -*a, fci, -P-, -P+, fed) 

" ^-^-4,-+^,^,- - _^ fcij _ p _ 5 (2 L25) 

• M cc/°NCd' <T2 ' ,73 ' fT4 (P+'P- ) k i,k 2 ,h } h) 

= M^-' ai,a2,a3 ' a4 (p + ,P-, h, kz, h, h) 

- MxQ a -' ai ' n '~' T+ ' /r *{-ka 1 p- 1 k 1 , k 2 , -p+, fa) 

- Mc + c :-^> am (p + ,p„, fa, fa, fa, fa) 

+ M° c + 67 A ^-°-^°\p + i -fa, fa, -p_, fc 3 , A*) 
+ Mc a 6 s f-' ai " T4 '~ a+ " T2 (-k3,P-, fa, fa, - P+ , k 2 ) 

- M^r ' ' " k,. -fa, fa, -p_, -p + , *?„). (2.1.26) 

The relative signs between contributions of the basic subamplitudes Alcca and A^NCa to 
the full matrix elements account for the sign changes resulting from interchanging external 
fermion lines. 

For the CC reactions, the amplitudes Afcca are the smallest gauge- invariant subset of 
diagrams JTI]]. In the case of NC reactions, the amplitudes AfNCa are composed of three 



separately gauge-invariant subamplitudes consisting of the first two lines, the two lines in 
the middle, and the last two lines of ( j2.1.19| ). 



Hadronic final states 

Next we inspect purely hadronic final states, i.e. the cases where all final-state fermions 
are quarks. This concerns only the channels CC(a), NC(a), NC(b), and CC/NC(a) given 
in Section |2.1.2| . The colour structure of the quarks leads to two kinds of modifications. 



Firstly, the summation of the squared amplitudes over the colour degrees of freedom can 
become non-trivial, and secondly, the possibility of virtual-gluon exchange between the 
quarks has to be taken into account. More precisely, there are diagrams of type (a) in 
Fig. ^71] in which one of the gauge bosons Vi j2 is a gluon. The other gauge boson of V^ 2 
can only be a photon or Z boson, since this boson has to couple to the incoming e + e~ 
pair. Consequently, there is an impact of gluon-exchange diagrams only for the channels 
NC(a), NC(b), and CC/NC(a), but not for CC(a). This can be easily seen by inspecting 
the generic diagrams in Fig. |2.1| : the presence of a gluon exchange requires two quark- 
antiquark pairs qq in the final state. 

We first inspect the colour structure of the purely electroweak diagrams. Since the 
colour structure of each diagram contributing to the basic channels CC(a) and NC(a) is 
the same, the corresponding amplitudes factorize into a simple colour part and the "colour- 
singlet amplitudes" Alcca and A^NCa, given in ( |2.1.16| ) and ( |2.1.19| ), respectively. The 
amplitudes for NC(b) and CC/NC(a) are composed from the ones of CC(a) and NC(a) in 
a way that is analogous to the singlet case, but now the colour indices q of the quarks 
have to be taken into account. Indicating the electroweak amplitudes for fully hadronic 
final states by "had,ew", and writing colour indices explicitly, we get 

A^cCa,had!ew^c 1 !c 2 X 1 C4(P+'-P-> fa, fa, fa, fa) 
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= Mc + c :-'^' a3 ^(p + ,p„, k x , k 2 , k 3 , k 4 ) 5 C1C2 6 C3C4 , (2.1.27) 

^NCa,^X 2 ,SCc 3 ,C4(P+'P-' fc l' fc 2, *3, 

= A<N + ca^* CT1 ' ,72 ' a3 * CT4 (P + ,P-, fe, **, fc 4 ) Wc 3C4 , (2.1.28) 

^NCJb,had 1 ,ew,'ci^C2*c3,c 4 (P+> P-> ^1, k 3 , £4) 

- A ^NCa"' <Tl ' ,72 ' fT3 ' <T4 (P+^- 5 fc 2, fo, ^4) 4 1C2 5 C3 C4 

- A^N + cr* CT3 * CT2 ' ,T1 ' a4 (P + ,P-, h, h, h, h) S CaC2 6 ClC4 , (2.1.29) 

^CC/ ff NCaXd?w?ci, C 2,C3, C 4(P+' P~' ^1) 

= -A^NCa _ '' T1 ' ,T2 ' ,7:i '' T4 (P+ ) P- ) fc 2, fcs, 5 C1C2 5 C3 C4 

- M a c + c :-> ai ' am ( P+ ,p„, h, h, h, k 2 ) S C1C4 5 C3C2 . (2.1.30) 

In the calculation of the gluon-exchange diagrams we can also make use of the "colour- 
singlet" result ( |2.1.8j ) for the generic diagram (a) of Fig. |2J], after splitting off the colour 
structure appropriately. Since each of these diagrams involves exactly one internal gluon, 
exchanged by the two quark lines, the corresponding matrix elements can be deduced in 
a simple way from the diagrams in which the gluon is replaced by a photon. The gluon- 
exchange contributions to the channels NC(b) and CC/NC(a) can again be composed from 
the ones for NC(a). Making use of the auxiliary function 

MZ^-^^( P+}P _ } k u fc 2 , k 3 , k 4 ) 
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2 



E [ My^^'-^i-h, -k 2 , P+ , P _, ~k 3 , ~h) 



v=-r,z 



+m v ; 3 


— (J^ ,U-\-,U 


- , — 01,— (72 


-fa, 


-fa, p. 


\-,P 


_, -A;i, -k 2 ) 




— CT2, — 0-3, 


— 04,o+,o_ 


-fa, 


-fa, ~ 


fa, 


-fa,P+,P-) 


+ M~? 


— <T4, — (Tl, 


-02,CT + ,CT_ 


-fa, 


-fa, - 


fa, 


-fa,P+,P-) 



(2.1.31) 



where g s = i/47ra s is the strong gauge coupling, the matrix elements involving gluon 
exchange explicitly read 

\A<y+,°- ,0-1,0-2 ,0-3,0-4 / 7 7 7 7 \ 

- /W NCa,had.gluon,ci,C2,C3,c 4 l/ : '+'-/ ; '-' ft 2 5 ^3; «-4j 



NCa, had, gluon, ci,C2,C3,C4 

^ + ' ff -' ffl,oa '° 3 ' ff4 (p + ,P-, fa, k 2 , k 3 , fa) J A^ 1C2 A^ 3C4 , (2.1.32) 



^NCb,had 1 ,gluon 3 ,ci 4 C2,C3,c 4 (P+?P-5 ^1) ^2, fa, faj) 

= A** + ' ff -' ff ™' ff4 (p + ,p_, fci, fa, k 3 , k 4 ) \ K lC2 K iCi 

- M^-^^( P+ ,p., fa, k 2 , fa, fa) \ X a C3C2 X a ClC4 , (2.1.33) 

»>0- + ,0_,01,02,03,04 / 7 7 7 7 \ 

- /Vt CC/NCa,had,gluon,ci,C2,c 3 ,C4^+'-P-' ft 2, «-3, «-4j 

^ + >->-i^3^ (p+)P _ )A . i)A . 2)A , 3)A . 4) 1 A« C2 A« 3C4 . (2.1.34) 

The colour structure is easily evaluated by making use of the completeness relation A^ A^ = 
— \8ij8u + 28u5jk for the Gell-Mann matrices A",-. 
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The complete matrix elements for the fully hadronic channels result from the sum of 
the purely electroweak and the gluon-exchange contributions, 

. >(T+,Cr_ ,(T1, (72, Cr3,(T4 . y<T-|-,<7_ ,(71,(72 ,(T3, (T4 , I / (7+ ,(7_ ,(71 ,(72 ,(73 ,(74 

J ...,had,ei,C2,C3,C4 •> 1 ...,had,ew,ci,C2,C3,C4 "T - ^ ...,had,gluon,ci,C2,C3,C4 

The gluon-exchange contributions are separately gauge-invariant. 

For clarity, we explicitly write down the colour-summed squared matrix elements 
for the fully hadronic channels. Abbreviating Ai r7+ ' cr -' cra,crb,crc,crd (p + ,p-, k a , kb, k c , ka) by 
A4...(a, b, c, d) we obtain 



(2.1.35) 



^ |-McCa,had (1, 2, 3, 4)| 2 = 9|M C Ca(l, 2, 3, 4)| 2 , (2.1.36) 

olour 

E l-^NCa,had(l, 2, 3, 4) | 2 = 9|A* NCa (l, 2, 3, 4) | 2 + 2|A4 g (l, 2, 3, 4) | 2 , (2.1.37) 



colour 



/j lA^NCb.had 

(1,2,3,4)| 2 

colour 

= 9|A4 NCa (l, 2, 3, 4) | 2 + 9|A^NCa(3, 2, 1, 4) | 2 - 6 Re {M NG a(l, 2, 3, 4)A^ Ca (3, 2, 1, 4)} 
+ 2|A< g (l,2,3,4)| 2 + 2|>l g (3,2,l,4)| 2 + |Re{M g (l,2,3,4)M*(3,2,l,4)} 
- 8 Re {Af N c a (l, 2, 3, 4)A1*(3, 2, 1, 4)} - 8 Re {Ai N Ca(3, 2, 1, 4)A<*(1, 2, 3, 4)}, 

(2.1.38) 

E l-McC/NCa,had(l,2,3,4)| 2 
colour 

= 9|A4 NCa (l, 2, 3, 4) | 2 + 9| A<cca(l, 4, 3, 2)| 2 - 6 Re {Al NC a(l, 2, 3, 4)A^ Ca (l, 4, 3, 2)} 
+ 2|A4 g (l, 2, 3, 4) | 2 - 8 Re {AW 1 , 4 , 3, 2)A4*(1, 2, 3, 4)}. (2.1.39) 

Owing to the colour structure of the diagrams, a non-zero interference between purely elec- 
troweak and gluon-exchange contributions is only possible if the four final-state fermions 
can be combined into one single closed fermion line in the squared diagram. This im- 
plies that fermion pairs must couple to different resonances in the electroweak and the 
gluon-exchange diagrams, leading to a global suppression of such interference effects in 
the phase-space integration (see Section [2.3.5|) . 



2.1.4 Generic functions and amplitudes for e + e — > 4f + 7 

The generic functions for e + e~ — > 4/7 can be constructed in a similar way. The idea is 
to combine the contributions of all those graphs to one generic function that reduce to 
the same graph after removing the radiated photon. These combined contributions to 
e + e~ — > 4/7 are classified in the same way as the diagrams for the corresponding process 
e + e~ — > 4/, i.e. the graphs of Fig. |2.1| also represent the generic functions for e + e~ — ► 4/7. 
Finally, all amplitudes for e + e~ — > 4/7 can again be constructed from only two generic 
functions. Note that the number of individual Feynman diagrams ranges between 14 and 
1008 for the various processes. We note that the generic functions can in fact be used to 
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construct the amplitudes for all processes involving exactly six external fermions and one 
external photon, such as e~e~ — > 4/7 and e ± 7 — > 5/. 

As a virtue of this approach, the so-defined generic functions fulfill the QED Ward iden- 
tity for the external photon, i.e. replacing the photon polarization vector by the photon 
momentum yields zero for each generic function. This is simply a consequence of elec- 
tromagnetic charge conservation. Consequently, in the actual calculation in the WvdW 
formalism the gauge spinor of the photon drops out in each contribution separately. 

Assuming the external fermions as incoming and the photon as outgoing, the generic 
functions read 

My 1 vf'" C '' Td,ae '' T, ' X (Qa, Qb, Qc, Qd, Qe, Qf, Pa, Pb, Pc, Pd, Pe, Pf, k) 

^ o aa - ab o ac - ad o ae - af yv 1 f a f g yv2f g f b yv 1 f c f d yv 2 feff 

X A^' a ^'\Q a ,Q b ,Q c ,Q d ,Q e ,Q f , Pa ,p b ,p c ,p d ,p e , Pf ,k), (2.1.40) 
My^^^iQa, Q b , Q c , Qd,Qe,Qf, Pa, Pb, Pc, Pd,Pe,Pf, k) 

= -4V2e 5 5 aa - r7b 5 (7c , + 5 ad -5 >7ei+ 5 af - (Q c - Qd)9vww9 v b f afb 9wf c f d 9wf e f f 

x A^ a ' x (Q a ,Q b ,Q c ,Q d ,Q e ,Q f ,p a ,p b ,p c ,p d ,p e ,p f ,k), (2.1.41) 

with the auxiliary functions 

^2 (Qa, Qb, Qc, Qd, Qe, Qf, 

Pa,Pb,Pc,Pd,Pe,Pf,k) = -(PaPc) 

Pv 1 (Pc+Pd)Pv 2 (Pe+Pf) 

'{PbPfY ( Qc-Qd (PaPc) 



X 



i(PbPd)*(PbPe) + (PdPf)*{PePf)) 



+ 
+ 



(p a k) \(p b +Pe+Pf) 2 (p c k) 

+ (^T^F !^ {{PaPdY{PaPe) + M)) 

Qb((PaPd)* (PaPe) + (PcPd)* (PcPe)) ((PbP f)* (PaPb) ~ (Pfk)*(p a k)) 

(Pa +Pc+ Pd) 2 (Pak) (p b k) 
(Qa + Qc~ Qd)(PbPf)*(PcPd)*(PaPc)((Pbk)*(p b Pe) ~ (Pfk)*(p e Pf)) 



(Pa+Pc+Pd) 2 (Pb+Pe+Pf) 2 (Pak) 

Qd - (Qc - Qd)2{Pd ■ k)P Vl (gg + pa) . M . v .„ 

■ ■ P Vl (p c +Pd- k)P V2 (Pe + Pf) (PbPf) 

(Pb+Pe+Pf) 2 

(PcPd)((PbPd)*(PbPe) + (PdPf)*(PePf)) + (Pck)((p b k)*{p h p e ) - (g fk)* (p e p /)) 

(p c k) (pdk) 



X 



(Pa+Pc+Pd) 2 

((PbPf)*(PePf) + (Pbk)* (p e k))((p a Pd)* (PaPc) + (PcPd)* {PcPe)] 

(Pek)(p f k) 

A 2 + ~ + (Qa, Qb, Qc, Qd, Qe,Qf, Pa, Pb, Pc, Pd, Pe,Pf, k) 

= A+ +++ (Q a , Q b , Q c , Q d , -Q f , -Qe,p a , Pb, Pc, Pd,Pf, Pc, k), 
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A\ ++ {Qa,Qb,Qc,Qd,Qe,Qf, Pa,Pb,Pc,Pd,Pe,Pf,k) 

= A2 +++ (Q a , Qb, — Qd, ~Qc, Qe, Qf, Pa,Pb,Pd,Pc,Pe,Pf,k), 
^2 ^{Qa,Qb,Qc,Qd,Qe,Qf, Pa,Pb,Pc,Pd,Pe,Pf,k) 

= A+ +++ (Q a , Q b , -Q d , -Q c , -Qf, -Q e ,p a ,p b ,p d ,p c ,p f ,p e , k), 

A 2 (Qa, Qbi Qc, Qdi Qe, Qf, 

Pa,Pb,Pc,Pd,Pe,Pf,k) 
= A+ +++ (Q b , Q a , -Q e , -Q f , -Q c , -Qd,Pb,Pa,Pe,Pf,Pc,Pd, k), 

A 2 + ~ + (Qa, Qb, Qc, Qd, Qe, Qf, Pa, Pb, Pc, Pd, Pe,Pf, k) 

= A% +++ (Q b , Q a ,Qf, Qe, -Qc, -Qd, Pb, Pa,Pf, Pc, Pc, Pd, k), 

^2 \Qa, Qb, Qc, Qd, Qe, Qf, 

Pa,Pb,Pc,Pd,Pe,Pf,k) 

= A^ +++ (Q b , Q a , -Q e , -Q f , Q d , Q c , p b , p a , p e ,p } , p d , p c , k), 

^2 {Qa, Qb, Qc, Qd, Qe, Q f, Pa,Pb,Pc,Pd,Pe,Pf,k) 

= At +++ (Q b , Qa, Qf, Qe, Qd, Qc, Pb, Pa, Pf, Pe, Pd, Pc, k) , 
A a a ,a c ,a d ,- ^ Q ^ Q ^ Q ^ Q ^ Q ^ ^ ^ ^ ^ ^ ^ ^ (2.1.42) 

= (A2 Ca ~ ac ~' Ji,+ {Qa, Qb, Qc, Qd, Qe, Qf,Pa, Pb, Pc, Pd, Pe,Pf, kj) 



and 



At + {Qa, Qb, Qc, Qd, Qe, Qf, Pa,Pb,Pc,Pd,Pe,Pf,k) 

= PviPa + Pb)Pw{Pc + Pd)Pw{Pe + Pf) ^~ k^p^k)^ 

X ( (PbPd}*(PbPf}*(PaPb}(PcPe) + (PbPd}*(PdPf}*(PaPe)(PcPd} 
+ {PbPf}*{PdPf}*(PaPc}(PePf}) 

+ Pv(Pa +Pb- k)P W (p c + Pd)P\v(Pe + Pf) J- 

X{ (PdPf}*[ {PaPc){PcPd)((PbPd}*{PaPb} - (Pdk}*(p a k)) 

+ (PaPc}(PcPf}((PbPf}*(PaPb) - (Pfk}*(p a k))} 

+ (PcPe}((PbPd}*(PaPb) - (Pdk)* (Pak})((p b Pf)* (PaPb) - (p f k}*(p a k})} 

pf . \ TD ( , 1\TD ( , \Qd~ (Qc-Qd)2(Pd- k)P W (Pc+Pd) 
+ PvVPa + Pb)Pw{Pc + Pd- k)P W {p e + Pf) k)(p d k) 

X{ (PbPf)*[ (PaPc}(PcPf}((PdPf}*(PcPd) - (Pfk}*(p c k)) 
+ (PcPe}(PaPb}((PbPd}*(PcPd) + (Pbk)* (Pck))} 

+ (PaPc}((PdPf}*(PcPd) - {Pfk}*(Pck})((pbPd)*(PcPd} + {Pbk}* (Pck))} 

+ Mp. + + Pd )P w (r< + p, - k) Q ' + <Qf ~^^ Pw<P " +Pl> 

X { (PbPd}*[ (PcPe}(PaPb}((PbPf}*(PePf) + (Pbk)*(p e k}) 
+ {PaPe}{PcPd)({PdPf}*(PePf} + (Pdk)* (p e k))} 
+ (PaPc}((PbPf}*(PcPf) + (Pbk)* (Pek))({p d Pf}* (PcPf) + (Pdk}*(p e k))}, 
As + (Q a , Qb, Qc, Qd, Qe, Qf, Pa, Pb, Pe, Pd, Pe,Pf, k) 

= ^3 + (~Qb, -Qa, Qc, Qd, Qe, Qf, Pb, Pa, Pc, Pd, Pe, Pf, k), 
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^3°' (Qa,Qb,Qc,Qd,Qe,Qf,PaiPb,PciPdiPeiPf,k) (2.1.43) 
Az aa,+ (Q a , Q b , -Q d , -Q c , -Q f , -Q e ,Pa,Pb,Pd,Pc,Pf,Pe, k) 



Pv,w{p)^Pvw^p) 



The replacements Py — > P v after the complex conjugation in the last lines of (|2.1.43|) 



and ( p.l.44| ) ensure that the vector-boson propagators remain unaffected. Note that the 
vector-boson masses do not enter explicitly in the above results, but only via Py. In 
gauges such as the 't Hooft-Feynman or the unitary gauge this feature is obtained only 
after combining different Feynman graphs for e + e~ — > 4/7; in the non-linear gauge ( |2.1.3|) 
this is the case diagram by diagram. 

The helicity amplitudes for e + e~ — > 4/7 follow from the generic functions AiviVz an d 
M.yww of C P-1 -41|) in exactly the same way as described in Section |2.1.3| for e + e~ — > 



4/. This holds also for the gluon-exchange matrix elements and for the colour factors. 
Moreover, the classification of gauge-invariant sets of diagrams for e + e~ — > 4/ immediately 
yields such sets for e + e~ — > 4/7, if the additional photon is attached to all graphs of a set 
in all possible ways. 

We have checked analytically that the electromagnetic Ward identity for the external 
photon is fulfilled for each generic contribution separately. In addition, we have numeri- 
cally compared the amplitudes for all processes with amplitudes generated by Madgraph 
|32|| for zero width of the vector bosons and found complete agreement. We could not 



compare our results with Madgraph for finite width, because Madgraph uses the unitary 
gauge for massive vector-boson propagators and the 't Hooft-Feynman gauge for the pho- 
ton propagators, while we are using the non-linear gauge (|2.1.3|) . Therefore, the matrix 
elements differ after introduction of finite vector-boson widths. While the calculation 
with Madgraph is fully automized, in our calculation we have full control over the matrix 
element and can, in particular, investigate various implementations of the finite width. 

A comparison of our results with those of Refs. |15|, [16| , which include only the matrix 
elements that involve two resonant W bosons, immediately reveals the virtues of our 
generic approach. 

2.1.5 Implementation of finite gauge-boson widths 

We have implemented the finite widths of the W and Z bosons in different ways: 

• fixed width in all propagators: Py{p) = [p 2 — My + iM v r v ] _1 , 

• running width in time-like propagators: Py{p) = [p 2 — M v + ip 2 (Fv/Mv)^(p 2 )] -1 , 



complex-mass scheme: complex gauge-boson masses everywhere, i.e. J M v — iMyTy 
instead of My in the propagators and in the couplings. This results, in particular, 
in a constant width in all propagators, 

Pvip) = [P 2 -M 2 + iM v T v }-\ (2.1.44) 

and in a complex weak mixing angle: 

c 2 =l-s 2 = M w - iM w r w (2 l 45) 

Cw Sw M|-iM z r z ' [ b) 
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The virtues and drawbacks of the first two schemes have been discussed in Ref. [10 



Both violate SU(2) gauge invariance, the running width also U(l) gauge invariance. The 
complex-mass scheme obeys all Ward identities and thus gives a consistent description of 
the finite-width effects in any tree-level calculation. While the complex-mass scheme works 
in general, it is particularly simple for e + e~ — > 4/7 in the non-linear gauge ( |2.1.3| ). In 
this case, no couplings involving explicit gauge-boson masses appear, and it is sufficient to 
introduce the finite gauge-boson widths in the propagators [cf. (|2.1.44 )] and to introduce 
the complex weak mixing angle ( 2.1.45Q in the couplings. We note that a generalization 
of this scheme to higher orders requires to introduce complex mass counterterms in order 
to compensate for the complex masses in the propagators f33fl . We did not consider the 
fermion-loop scheme |7], |], |10 , Q , which is also fully consistent for lowest-order predictions, 
since it requires the calculation of fermionic one-loop corrections to e + e~ — > 4/7 which is 
beyond the scope of this work. 



2.2 The Monte Carlo program 

The cross section for e + e~ — > 4/(7) is given by 
(2vr) 4 - 3 ™ " 



da 



=1 

|2 




2s 

x \M(p + , P -,h,...,k n )\ 2 , (2.2.1) 



where n = 4, 5 is the number of outgoing particles. The helicity amplitudes M. for e + e~ — > 
4/(7) have been calculated in Sections [2.1.3| and |2.1.4j . The phase-space integration is 
performed with the help of a Monte Carlo technique, since the Monte Carlo method 
allows us to calculate a variety of observables simultaneously and to easily implement cuts 
in order to account for the experimental situation. 

The helicity amplitudes in (|2.2.1|) exhibit a complicated peaking behaviour in different 
regions of the integration domain. In order to obtain a numerically stable result and to 
reduce the Monte Carlo integration error we use a multi-channel Monte Carlo method 



f23| , |26| , which is briefly outlined in the following. 

Before turning to the multi-channel method, we consider the treatment of a single chan- 
nel. We choose a suitable set <E» of 3n — 4 phase-space variables to describe a point in phase 
space, and determine the corresponding physical region V and the relation /cj($) between 
the phase-space variables $ and the momenta k\, . . . , k n . The phase-space integration of 
( PXID reads 

I n = Jda = j^d* />(**(*)) /(**(*)), (2.2.2) 



/o \4-3n 
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where p is the phase-space density. For the random generation of the events, we further 
transform the integration variables $ to 3n — 4 new variables r = (rj) with a hypercube 
as integration domain: <fr = h(r) with < < 1. We obtain 



J n = / d$p (*,(*)) /(*,(*)) = / dr }- ( , (2.2.3) 

•/V ' Jo g(ki(h(r))) 



where g is the probability density of events generated in phase space, defined by 

<9h(r 



—r— \~ = p(h(& 



dr 



(2.2.4) 

r=h- 1 (*) 



If / varies strongly, the efficiency of the Monte Carlo method can be considerably enhanced 
by choosing the mapping of random numbers r to $ in such a way that the resulting density 
g mimics the behaviour of |/|. For this importance sampling, the choice of $ is guided 
by the peaking structure of /, which is determined by the propagators in a characteristic 
Feynman diagram. 

We choose the variables $ in such a way that the invariants corresponding to the 
propagators are included. Accordingly, we decompose the ra-particle final state into 2^2 
scattering processes with subsequent 1 —>■ 2 decays. The variables $ consist of Lorentz 
invariants defined as the squares of time- and space-like momenta, respectively, 

and of polar and azimuthal angles 0j and (pi, defined in appropriate frames. A detailed 
description of the parameterization of an n-particle phase space in terms of invariants 
and angles can be found in Appendix 0. The parameterization of the invariants Sj and 
U in $ = h(r) is chosen in such a way that the propagator structure of the function / is 
compensated by a similar behaviour in the density g. More precisely, if / contains Breit- 
Wigner resonances or distributions like s^ u , which are relevant for massless propagators, 
appropriate parameterizations of Sj are given by: 

• Breit-Wigner resonances: 

Sl = i\4 + M v r v tanfe/! + (?/ 2 - yi)rj (2.2.5) 

with yx 2 = arctan 



' q . — M 2 * 

°mm,max -'"V 



v M v r v 

propagators of massless particles: 
v ^ 1 : Si 



si-V + s 1 T"('l - r - 



1/(1-") 



max * 

v = 1 : Si = exp [ln(s max )rj + ln(s min )(l - n)} . (2.2.6) 

For the choice of v see Appendix 0. The remaining variables in $ = h(r), i.e. those for 
which / is expected not to exhibit a peaking behaviour, are generated as follows: 

Si = s max ri + s mm (l - ri), 4>i = 2nri, cos 0* = 2r { , - 1 . (2.2.7) 
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The absolute values of the invariants tj are generated in the same way as Sj. The result- 
ing density g of events in phase space is obtained as the product of the corresponding 
Jacobians, as given in ( 2.2.4j) . In the Appendix |Cj, we provide explicit examples for the 
generation of events with a specific choice of mappings and h(r), and for the calcu- 

lation of the corresponding density g. 

The differential cross sections of the processes e + e~ — > 4/ and especially e + e~ — > 4/7 
possess very complex peaking structures so that the peaks in the integrand /(«&) in ( |2.2.3[ ) 
cannot be described properly by only one single density g(*&)- The multi-channel approach 
2"B| , pB[ suggests a solution to this problem. For each peaking structure we choose a suitable 
set and accordingly a mapping of random numbers into <&k = hfc(r) with 
< Tj < 1, so that the resulting density g^ describes this particular peaking behaviour 
of /. All densities g k are combined into one density g to t that is expected to smoothes 
the integrand over the whole phase-space integration region. The phase-space integral of 
( gggfl reads 



M 



k=l 



l n = E / d& k Pk(H®k)) QkiH® 



with 



M 

gtot(ki(®k)) = E 



1=1 



gtot[ki{®k 



1 

gi(ki(& k ) 



M 

E 

fe=i 



1 , /(^(h,(r 
dr 







Sw(^(h fc (r)) 



, (2.2.8) 



Pi(ki(® k 



dhi(r) 



Or 



r=h-\^ k ) 



(2.2.9) 



The different mappings h^(r) are called channels, and M is the number of all channels. 

In order to reduce the Monte Carlo error further, we adopt the method of weight 
optimization of Ref. [27 and introduce a-priori weights a k ,k = 1, . . . , M (a k > and 
Z)fc=i a k = 1)- The channel k that is used to generate the event is picked randomly with 
probability a k , i.e. 



M 



fe=l 



/ dr J2 9 ( r o-/3 k -i)e(/3 k -r ) 

M 

dr E d ( r o ~ Pk-iWk - r ) 



h=i 



g to t(ki(®k 
d$ k pkfkii&k^gkfkii&k] 
f(ki(h h (r)) 



dr 







S'tot(^(hfe(r)) 



(2.2.10) 



where (3q = 0, (3j 



M-1,/3 M 

M 



1, and 



;2.2.n) 



is the total density of the event. 

For the processes e + e~ — > 4/ we have between 6 and 128 different channels, for e + e~ — > 
4/7 between 14 and 928 channels. Each channel smoothes a particular combination of 
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propagators that results from a characteristic Feynman diagram. We have written phase- 
space generators in a generic way for several classes of channels determined by the chosen 
set of invariants Sj, tj. The channels within one class differ in the choice of the mappings 
Q2.2.5D , ( |2.2.60 , an d Q2.2.7Q and the order of the external particles. We did not include 



special channels for interference contributions. 
The a^-dependence of the quantity 



1 N 

W(a) = -Y,Hri^)} 2 , (2.2.12) 

iV 3=1 



where w = f /gtot is the weight assigned to the Monte Carlo point (r J , r- 7 ) of the jth event, 
can be exploited to minimize the expected Monte Carlo error 



a. = Jm^z, (2 . 2 , 3 ) 



with the Monte Carlo estimate of 



1 N 



I n = ^E w ( r o^ 3 ) (2-2.14) 
by trying to choose an optimal set of a-priori weights. We perform the search for an 



optimal set of by using an adaptive optimization method, as described in Ref. |27 . 
After a certain number of generated events a new set of a-priori weights a£ ew is calculated 
according to 



,.new 



1 N gJUhJri))) \w(rLrj)} 2 m 



a k oc oik 



Based on the above approach, we have written two independent Monte Carlo programs. 
While the general strategy is similar, the programs differ in the explicit phase-space gen- 
eration. 

2.3 Numerical results 

If not stated otherwise we use the complex-mass scheme and the following parameters: 

a = 1/128.89, a, = 0.12, 

M w = 80.26 GeV, T w = 2.05 GeV, 

M z = 91.1884 GeV, T z = 2.46 GeV. (2.3.1) 



In the complex-mass scheme, the weak mixing angle is defined in ( 2.1.45Q , in all other 
schemes it is fixed by c w = M w /M z , = 1 — c^. 
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The energy in the centre-of-mass (CM) system of the incoming electron and positron 
is denoted by yfs. Concerning the phase-space integration, we apply the canonical cuts of 
the ADLO/TH detector, 

6(1, beam) > 10°, 6(1,1') > 5°, 6(l,q)>5°, 

0(7, beam) > 1°, 6(j,l) > 5°, 6(j,q) > 5°, 

£ 7 > 0.1 GeV, Ei > 1 GeV, E q > 3 GeV, 

m(q, q') > 5 GeV, (2.3.2) 

where 8(i,j) specifies the angle between the particles i and j in the CM system, and /, 
q, 7, and "beam" denote charged leptons, quarks, photons, and the beam electrons or 
positrons, respectively. The invariant mass of a quark pair qq' is denoted by m(q, q'). The 
cuts coincide with those defined in Ref . PU| , except for the additional angular cut between 



charged leptons. The canonical cuts exclude all collinear and infrared singularities from 
phase space for all processes. 

Although our helicity amplitudes and Monte Carlo programs allow for a treatment of 
arbitrary polarization configurations, we consider only unpolarized quantities. 

All results are produced with 10 events. The calculation of the cross section for 
e + e~ — > e + e~ fi + fi~ requires about 50 minutes on a DEC ALPHA workstation with 500 
MHz, the calculation of the cross section for e + e~ — > e + e~ /i + yU~7 takes about 5 hours. 
The numbers in parentheses in the following tables correspond to the statistical errors of 
the results of the Monte Carlo integrations. 

2.3.1 Comparison with existing results 

In order to compare our results for e + e~ — > 4/ with Tables 6-8 of Ref. we use the 
corresponding set of phase-space cuts and input parameters, i.e. the canonical cuts de- 
fined in (|2.3.2|) , a CM energy of y/s = 190 GeV, and the parameters a = a(2M w ) = 



1/128.07, a s = 0.12, M w = 80.23 GeV, T w = 2.0337 GeV, M z = 91.1888 GeV, and 
Tz = 2.4974 GeV. The value of s w , which enters the couplings, is calculated from 
a(2M w )/(2sl) = G M M^/(7rV2) with = 1.16639 x lO^GeV" 2 . 

In Table [2.1| , we list the integrated cross sections for various processes e + e~ — > 4/ with 
running widths and constant widths, and for the corresponding processes e + e~ — > 4/7 
with constant widths. For processes involving gluon-exchange diagrams we give the cross 
sections resulting from the purely electroweak diagrams and those including the gluon- 
exchange contributions. The latter results include also the interference terms between 
purely electroweak and gluon-exchange diagrams. In Table |2j] we provide a complete list 
of processes for vanishing fermion masses. All processes e + e~ — > 4/(7) not explicitly listed 
are equivalent to one of the given processes. 

For NC processes e + e~ — > 4/ with four neutrinos or four quarks in the final state 
we find small deviations of roughly 0.2% between the results with constant and running 



widths. Assuming that a running width has been used in Ref. |3"5| , we find very good 
agreement. 

Unfortunately we cannot compare with most of the publications []T7|, [T^, |T^, ^] for 
the bremsstrahlung processes e + e~ — > ^f^- In those papers, either the cuts are not 
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cr/fb 


e+e" -> 4/ 


e+e" -> 4/ 


e + e — > 4/7 




running width 


constant width 


constant width 


z/ e z/ e e _ e + 


256.7(3) 


257.1(7) 


89.4(2) 


z/ M /z + e~P e 


227.4(1) 


227.5(1) 


79.1(1) 




228.7(1) 


228.8(1) 


81.0(2) 




218.55(9) 


218.57(9) 


76.7(1) 


e~e + e~e + 


109.1(3) 


109.4(3) 


38.8(4) 


e~e + /i~/i + 


116.6(3) 


116.4(3) 


43.4(4) 




5.478(5) 


5.478(5) 


3.37(1) 




11.02(1) 


11.02(1) 


6.78(3) 




14.174(9) 


14.150(9) 


5.36(1) 




17.78(6) 


17.73(6) 


6.63(2) 




10.108(8) 


10.103(8) 


4.259(9) 


v c v c v c v c 


4.089(1) 


4.082(1) 


0.7278(7) 




8.354(2) 


8.337(2) 


1.512(1) 




4.069(1) 


4.057(1) 


0.7434(7) 




8.241(2) 


8.218(2) 


1.511(1) 


ude~z/ e 


693.5(3) 


693.6(3) 


220.8(4) 


ud/x"P /i 


666.7(3) 


666.7(3) 


214.5(4) 


e~e + u u 


86.87(9) 


86.82(9) 


32.3(2) 


e~e + dd 


43.02(4) 


42.95(4) 


16.17(8) 


U U fi~ fi + 


24.69(2) 


24.69(2) 


12.70(4) 


d d fi~ fi + 


23.73(1) 


23.73(1) 


10.43(2) 


z/ e z/ e uu 


24.00(2) 


23.95(2) 


6.84(1) 


z/ e z/ e d d 


20.657(8) 


20.62(1) 


4.319(6) 


u u z/ M z/ M 


21.080(5) 


21.050(5) 


6.018(9) 


d d v^v^ 


19.863(5) 


19.817(5) 


4.156(5) 


uudd 


2064.1(9), 2140.8(9) 


2064.3(9), 2141(1) 


615(1), 672(1) 


udsc 


2015.2(8) 


2015.3(8) 


598(1) 


uuuu 


25.738(7), 71.28(4) 


25.721(7), 71.30(4) 


9.78(2), 42.1(1) 


dddd 


23.494(6), 51.35(3) 


23.448(6), 51.32(3) 


5.527(7), 28.68(4) 


uucc 


51.61(1), 144.72(9) 


51.57(1), 144.75(9) 


19.61(4), 86.1(2) 


uuss 


49.68(1), 126.52(8) 


49.62(1), 126.52(8) 


15.17(2), 75.1(2) 


dds s 


47.13(1), 104.79(6) 


47.02(1), 104.74(6) 


11.10(2), 59.2(1) 



Table 2.1: Integrated cross sections for all representative processes e + e~ — > 4/ with run- 
ning widths and constant widths and for the corresponding processes e + e~ — > 4/7 with 
constant widths. If two numbers are given, the first results from pure electroweak diagrams 
and the second involves in addition gluon-exchange contributions. 
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a/fb 


8 ■ E ■ 


1 GeV 

X VJC V 


5 OeV 


10 OeV 

1 W VJC V 


1 5 OeV 


leptonic 

process 


1 ° 
1 


0o.04^o ) 


97 ^7fQA 


lo.yo^z ) 


1 l.ZZ^Z j 




32 65(4") 


16 98(31 

1U. Jul U J 


10 48(21 


6 94(1 1 


i n° 

x w 


23 48f3l 


12 30(21 


7 61 (21 


5 04(1 1 

U.Ut:! X J 


1 5° 


1 8 (13(2 1 


9 ^1 (21 


5 QDM 1 


3 90(1 1 


semi-leptonic 
process 


1 


141. y ^Zj 


( l.yU(o J 


43.00(^0 J 


no r>£; / /I A 

Z8.zo^4J 


5° 


86.8(1) 


44.25(6) 


26.78(4) 


17.40(3) 


10° 


62.29(7) 


31.92(5) 


19.40(4) 


12.61(3) 


15° 


47.42(5) 


24.50(4) 


14.97(3) 


9.77(2) 



Table 2.2: Comparison with Table 2 of Ref. jl5| : Cross sections resulting from diagrams 
involving two resonant W bosons for purely leptonic and semi-leptonic final states and 
several photon separation cuts 



(completely) specified, or collinear photon emission is not excluded, and the corresponding 
fermion-mass effects are taken into account. Note that the contributions of collinear 
photons dominate the results given there. 

We have compared our results with the ones given in Refs. |15|, |16[ , where the total cross 
sections for e + e~ — > Af^j have been calculated for purely leptonic and semi-leptonic final 
states. As done in Refs. [15], [LJl only diagrams involving two resonant W bosons have been 



taken into account for this comparison. Table |2^ contains our results corresponding to 
Table |2.2j of Ref. [|15j . Based on Refs. [^, [16| , we have chosen = 200 GeV and the input 
parameters a = 1/137.03599, M w = 80.9 GeV, T w = 2.14 GeV, M z = 91.16 GeV, T z = 
2.46 GeV, s w obtained from a/(2s 2 w ) = G^M^/(nV2) with G M = 1.16637 x lO^GeV -2 , 
and constant gauge-boson widths. The energy of the photon is required to be larger 
than E 1>m \-n , and the angle between the photon and any charged fermion must be larger 
than 9 1)min . A maximal photon energy is required, E^ < 60 GeV, in order to exclude 
contributions from the Z resonance. Our results are consistent with those of Refs. [^, [IT 
within the statistical error of 1% given there. In some cases we find deviations of 2%.[] 



2.3.2 Comparison of finite-width schemes 

As discussed in Refs. M, 91 [LT], [54]] , particular care has to be taken when implementing the 
finite gauge-boson widths. Differences between results obtained with running or constant 
widths can already be seen in Table EO, where a typical LEP2 energy is considered. In 
Table |273| we compare predictions for integrated cross sections obtained by using a constant 
width, a running width, or the complex-mass scheme for several energies. We consider 



1 Note that the input specified in Refs. |H| [T(| is not completely clear even if the information of both 
publications is combined. 
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cr/fb 


V ° 


189 GeV 


500 GeV 

.j v / v / * — a v . v 


2TeV 


lOTeV 


e + e~ — ► ud/i"^ 


^vjlloudlllj W ILiljll 


703 5C3i 


237 4C1 1 

jLO 1 .11 1 I 


1 3 QQC91 


(S94C31 


rnnnin o* wiHf n 


703 4C3") 


238 9C1 ") 

iUU. Jl 1 1 


34 39 C3") 


498 8C1 1 

t: JU. ul 1 1 


complex-mass scheme 


703.1(3) 


237.3(1) 


13 98C2) 


624C3) 


e + e~ — ► ud/i~^7 


LUllO LjCLIIu VV 1U. u 11 


224 0C41 


83 4C31 


6 98C51 


457C61 


T'l 1 1 1 1 1 1 1 1 o' yx n H f" n 

I Lilllllllfci W1U.L11 


994 fiC4"l 


84 9CV) 


1 Q 9C1 1 


dUOl U 1 


complex-mass scheme 


223.9(4) 


83.3(3) 


6.98(5) 


460(6) 


e + e~ — » ude~z/ e 


pnn fi wini"n 

LUllO IjdllL VV 1U. ull 


730 PCS') 


395 3C9"l 


91 1 0C9"l 


39 38C6) 


Tim Timer winfn 

1 Lilllllllfci VVlU-ull 


729 8C3"l 


396 9C2"l 


231 5C2"l 

Zvd l.ul Z J 


530 9C6"l 


POTTl TU P^T-TYl ^ cm QPrl PTT1 P 
v. V.y 1 ± 1 1 V 1 V~A IIICIOO O^llV^lllV^ 


729 8C3") 


395 1 C2l 


210 9C2") 


32 37C6") 


e + e~ — > ude~P e 7 


constant width 


230.0(4) 


136.5(5) 


84.0(7) 


16.8(5) 


running width 


230.6(4) 


137.3(5) 


95.7(7) 


379(6) 


complex-mass scheme 


229.9(4) 


136.4(5) 


84.1(6) 


16.8(5) 



Table 2.3: Comparison of different width schemes for several processes and energies 



two semi-leptonic final states for e + e~ — > 4/(7). The numbers show that the constant 
width and the complex-mass scheme yield the same results within the statistical accuracy 
for e + e~ — > 4/ and e + e~ — > 4/7. In contrast, the results with the running width produce 
totally wrong results for high energies. The difference of the running width with respect 
to the other implementations of the finite width is up to 1% already for 500 GeV. Thus, 
the running width should not be used for linear-collider energies. As already stated above, 
our default treatment of the finite width is the complex-mass scheme in this chapter. 



2.3.3 Survey of photon-energy spectra 

In Fig. |2.2| we show the photon-energy spectra of several processes for the typical LEP2 
energy of 189 GeV and a possible linear-collider energy of 500 GeV. The upper plots 
contain CC and CC/NC processes, the plots in the middle and the lower plots contain 
NC processes. Several spectra show threshold or peaking structures. These structures 
are caused by diagrams in which the photon is emitted from the initial state. The two 
important classes of diagrams are shown in Fig. p. 31 



The first class, shown in Fig. |2.3| a, corresponds to triple-gauge-boson-production sub- 
processes which yield dominant contributions as long as the two virtual gauge bosons 
V\ and V2 can become simultaneously resonant. If the real photon takes the energy -Ky, 



defined in the CM system, only the energy v s', with 

s' = s-2v^£ 7 , (2.3.3) 
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Figure 2.2: Photon-energy spectra for several processes and for ^/s 
500 GeV 



189 GeV and 



26 




b) 




Figure 2.3: Diagrams for important subprocesses, where V\, V2 = W, Z, 7, and V3 = 7, g 
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Table 2.4: Photon energies EY lV2 corresponding to thresholds 



is available for the production of the gauge-boson pair V1V2. If at least one of the gauge 
bosons is massive, and if the photon becomes too hard, the two gauge bosons cannot be 
produced on shell anymore, so that the spectrum falls off for above the corresponding 
threshold E^ lV2 . Using the threshold condition for the on-shell production of the V1V2 
pair, 



Vs* > M Vl + M 



l V 2 : 



(2.3.4) 



the value of EY 1 2 is determined by 



7 



s - (M Vl + M 



V2, 



2^1 



(2.3.5) 



The values of the photon energies that cause such thresholds can be found in Table |2j4 . 
The value E^ 1 corresponds to the upper endpoint of the photon-energy spectrum, which 
is given by the beam energy yfs/2. Since ^fs' is fully determined by s and _E 7 , the 
contribution of the V\ ^-product ion subprocess to the E y spectrum qualitatively follows 
the energy dependence of the total cross section for V1V2 production (cf. Ref. p5| , Fig. 1) 
above the corresponding thresholds. The cross sections for 77 and 7Z production strongly 
increase with decreasing energy, while the ones for ZZ and WW production are comparably 
flat. Thus, the 77 and 7 Z-production subprocesses introduce contributions in the photon- 
energy spectra with resonance- like structures, whereas the ones with ZZ or WW pairs 
yield edges. 

The second class of important diagrams, shown in Fig. |2.3| b corresponds to the produc- 
tion of a photon and a resonant Z boson that decays into four fermions. These diagrams 
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are important if the gauge boson V3 is also resonant, i.e. a photon or a gluon with small 
invariant mass. In this case, the kinematics fixes the energy of the real photon to 



.7 s-MS , 



which corresponds to the 7Z threshold in Table |2.4f This subprocess gives rise to resonance 
structures at E^ z , which are even enhanced by a s /a in the presence of gluon exchange. 

In the photon-energy spectra of Fig. [2.2| all these threshold and resonance effects are 
visible. The effect of the 7Z peak can be nicely seen in different photon-energy spectra, 
in particular in those where gluon-exchange diagrams contribute (cf. also Fig. |2.5|) . The 



effect of the WW threshold is present in the upper two plots of Fig. [2.2| . In the plot 
for y/s = 189 GeV the threshold for single W production causes the steep drop of the 
spectrum for the pure CC processes above 70 GeV. Note that the CC cross sections are 
an order of magnitude larger than the NC cross sections if the WW channel is open. The 
ZZ threshold is visible in the middle and lower plots for ^Js = 500 GeV. The 7Z threshold 
(resulting from the graphs of Fig. |2.3| a) is superimposed on the 7Z peak (resulting from 
the graphs of Fig. |2.3| b) and therefore best recognizable in those channels where the 7Z 
peak is absent or suppressed, i.e. where a neutrino pair is present in the final state or 
where at least no gluon-exchange diagrams contribute. Processes with four neutrinos in 
the final state do not involve photonic diagrams and are therefore small above the ZZ 
threshold. The effects of the triple-photon-production subprocess appear as a tendency 
of some photon-energy spectra to increase near the maximal value of for two charged 
fermion-antifermion pairs in the final state. 



2.3.4 Triple-gauge-boson-production subprocesses 

In Fig. [2.4| we compare predictions that are based on the full set of diagrams with those 
that include only the graphs associated with the triple-gauge-boson-production subpro- 
cesses, i.e. the graphs in Fig. |2.3| a. In addition we consider the contributions of the ZZ7- 
production subprocess alone. For CC processes, the photon-energy spectra resulting from 
the W + W~7-production subprocess are close to those resulting from all diagrams at LEP2 
energies, but large differences are found for higher energies and e in the final state. 
Note that the spectra are shown on a logarithmic scale. Even at LEP2 energies the dif- 
ferences between the predictions for different final states may be important, as can be 
seen, for instance, in Table ^TT] by comparing the cross sections of e + e~ — > and 
e + e~ — > ude~z/ e 7. In the case of NC processes, already for 189 GeV the contributions from 
ZZ7, Z77, and 777 production are not sufficient: in the vicinity of the 7Z peak sizeable 
contributions result from the 7Z-production subprocess (Fig. |2~3|b) even for the /i + /i~uu7 
final state. For e + e~ — > e~e + uu7 other diagrams become dominating everywhere. The 
contribution of ZZ7 production is always small and could only be enhanced by invariant- 
mass cuts. Note that the triple-gauge-boson-production diagrams form a gauge-invariant 
subset for NC processes, while this is not the case for CC processes. 
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Figure 2.4: Photon-energy spectra resulting from the triple-gauge-boson-production sub- 
processes compared to those resulting from all diagrams (V1V2 includes ZZ, 7Z, and 77) 
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Tablp 2.5: Full lowpst ordpr cross spction (pw and gluon) and contributions of purely 
elpctrowcak diagrams (pw) , of gluon-pxchangp diagrams (gluon) , and thpir intprfprpnee for 
500 GpV 

2.3.5 Relevance of gluon-exchange contributions 

In thp analytical calculation of the matrix plpmpnts for p + p~ — > 4/(7) in Spction ^4] wc 
havp sppn that NC proccssps with four quarks in thp final statp involvp, bpsidps purply 



plpctrowcak, also gluon-pxchangp diagrams. Tablp |2.5| illustratps thp impact of thpsp dia- 
grams on thp intpgratpd cross sections for a CM pnprgy of 500 GpV. Thp results for the 
intprfprpnee are obtained by subtracting the purely electroweak and the gluon contribution 
from the total cross section. For pure NC processes the contributions of gluon-exchange 
diagrams dominate over the purely electroweak graphs. This can be understood from 
the fact that the gluon-exchange diagrams are enhanced by the strong coupling constant, 



and, as discussed in Section |2.3.3| , that the diagrams with gluons replaced by photons 
yield a sizeable contribution to the cross section. For the mixed CC/NC processes the 
purely electroweak diagrams dominate the cross section. Here, the contributions from the 
W + W~7-pro duct ion subprocess are large compared to all other diagrams, even if the latter 
are enhanced by the strong coupling. At 500 GeV the gluon-exchange diagrams contribute 
to the cross section at the level of several per cent. The interference contributions are 



relatively small. As discussed at the end of Section 2.1.3, this is due to the fact that inter- 



fering electroweak and gluon-exchange diagrams involve different resonances. Note that 



the interference vanishes for e e — > u u c c 7, and the corresponding numbers in Table |2T5 
are only due to the Monte Carlo integration error. 



In Fig. |2.5| we show the photon-energy spectra for the processes e + e — > uudd7 and 



e + e~ — > u u u u 7 together with the separate contributions from purely electroweak and 
gluon-exchange diagrams. The pure electroweak contributions are similar to the ones 
for e + e~ — > ud/i"^7 and e + e~ — > e~e + uu7 in Fig. |2T2| For the NC process e + e~ — > 
UUUU7, the photon-energy spectrum is dominated by the gluon-exchange contribution, 
which shows a strong peak at 72.5 GeV owing to the 7Z-production subprocess. For the 
CC/NC process e + e _ — > uudd7, the electroweak diagrams dominate below the WW 
threshold, whereas the gluon-exchange diagrams dominate at the 7Z peak and above. The 
interference between purely electroweak and gluon-exchange diagrams is generally small. 
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Figure 2.5: Electroweak and gluon-exchange contributions to the photon-energy spect 
for e + e~ — > u u d d 7 and e + e~ — > u u u u 7 at y/s = 189 GeV 
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Chapter 3 



Non-factorizable photonic 
corrections to 



In this chapter we define and explicitly calculate the non-factorizable photonic corrections. 
It has been shown that they vanish in inclusive quantities, i.e. if the invariant masses of 
both W bosons are integrated out fl37Jf. However, for non-inclusive quantities these cor- 



rections do not vanish in general. The non-factorizable photonic corrections have already 
been investigated by two groups. Melnikov and Yakovlev |37| have given the analytical 
results only in an implicit form and restrict the numerical evaluation to a special phase- 
space configuration. Beenakker, Berends and Chapovsky have provided both the complete 
formulae and an adequate numerical evaluation p8| , |39[] , but do not find agreement with 
all results of Ref. |37j] . For this reason, it is worth-while to present the results of a third 
independent calculation. The material of this chapter has been published in Refs. KO, 



We start out by discussing the definition of the virtual and real photonic non-factoriz- 
able corrections in double-pole approximation in detail. Since only soft photons are rele- 
vant in double-pole approximation, the virtual non-factorizable correction is just a factor 
to the lowest-order cross section. For the corresponding real correction, the situation is 
similar, but in addition an integration over the photon momentum has to be performed. 
This requires a specification of the phase-space parameterization, which includes, in par- 
ticular, the invariant masses of the W bosons. Usually these are defined via the invariant 
masses of the respective final-state fermion pairs and are chosen as independent variables 
P?| , |3"8| , Experimentally, however, the invariant mass of a W boson is identified with 
the invariant mass of the associated jet pair that necessarily includes soft and collinear 
photons. Therefore, the influence of the choice for the invariant masses of the W bosons 
on the non-factorizable corrections should be investigated in order to provide sound pre- 
dictions for physical situations. 

Besides the non-factorizable doubly-resonant corrections, the most important effect of 
the instability of the W bosons is the modification of the Coulomb singularity. Since the 
off-shell Coulomb singularity results from a scalar integral that also contributes to the 
doubly-resonant non-factorizable corrections, it seems to be natural to approximate this 
integral in such a way that both effects are simultaneously included. This requires going 
beyond the strict double-pole approximation. 
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3.1 Definition of the approximation 

3.1.1 Conventions and notations 

We discuss corrections to the process 

e + (p + )+e-(p_) - W + (A; + )+W-(A;_) - frfa) + f 2 (k 2 ) + f 3 (k 3 ) + fr(h). (3.1.1) 

The relative charges of the fermions /j are represented by Qi with % — 1, . . . , 4. The masses 
of the external fermions, mf = k 2 and m 2 = p±, are neglected, except where this would 
lead to mass singularities. The momenta of the intermediate W bosons are defined by 

k + = k\ + k 2 , k_ = fc 3 + fc 4 , (3.1.2) 

their complex mass squared and their respective invariant masses are denoted by 

M' 2 M'i- - iM w r w , M± = y/k[, (3.1.3) 
respectively, and we introduce the variables 

K + = k 2 + - M 2 , K_ = k 2 _- M 2 . (3.1.4) 
Furthermore, we define the following kinematical invariants 



and 



t=(p±- 


k±) 2 , u = 


ip± - k T ) 2 




t+t = (p+ - 


k%) 2 , u_i = 


(P- " h) 2 , 


2 = 1,2, 


t-i = (p- - 


k%) 2 , it+i = 


(p + - ki) 2 , 


2 = 3,4, 


s = 


(p++p-) 2 = (k. 






Sij = 


(ki + kj) , 






&ijk 


(ki + kj + fcfc) , 


hJ,k = 


1,2,3,4, 



(3.1.5) 



(3.1.6) 



which obey the relations 



S = k\ + k 2 _ + Si 3 + S14 + S 23 + S 2 4, Sl2 = S 34 = ^, 

Sijk = s i: j + s ik + s jk , k = 1,2,3,4. (3.1.7) 



3.1.2 Doubly-resonant virtual corrections 

The aim of this chapter is to evaluate the non-factorizable corrections to the process ( |3.1.1|) 
in double-pole approximation (DPA). The DPA takes into account only the leading terms 
in an expansion around the poles originating from the two resonant W propagators. 
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(a) type (mf') (b) type (mf') 




Figure 3.1: Examples of non-factorizable photonic corrections in 0(a). The shaded blobs 
stand for all tree-level graphs contributing to e + e~ — ► W + W~. Whenever Feynman dia- 
grams with intermediate would-be Goldstone bosons ^ instead of W ± bosons are relevant, 
the inclusion of such graphs is implicitly understood. 



In DPA, the lowest-order matrix element for the process ( [3.1.1|) factorizes into the ma- 
trix element for the on-shell W-pair production, A^^ n ^ w+W (p+,p^ } k + , the (trans- 
verse parts of the) propagators of these bosons, and the matrix elements for the decays of 
these on-shell bosons, J\A Boz ^ (k + , k 2 ) and .M^ rn ~^ 3 ^ 4 (/c_, k 3 , k A ): 

Ayfe+e-^W+W" A/f W+ ^/i/2 i^W'-^fah 
KA - /W Born ^Born - /w Born fn -i s \ 

-MBorn = 2^ ' (3.1.8) 

A+,A_ K + K ~ 

The sum runs over the physical polarizations \± of the W ± bosons. 

The higher-order corrections to fl3.1.1|) can be separated into factorizable and non- 
factorizable contributions || [TJ, In the factorizable contributions the production of 
two W bosons and their subsequent decays are independent. The corresponding Feynman 
diagrams can be split into three parts by cutting only the two W-boson lines. The cor- 
responding matrix element factorizes in the same way as the lowest-order matrix element 

(MM- 

The non-factorizable corrections comprise all those contributions in which W-pair pro- 
duction and/or the subsequent W decays are not independent. Obviously, this includes 
all Feynman diagrams in which a particle is exchanged between the production subprocess 
and one of the decay subprocesses or between the decay subprocesses. Examples for such 
manifestly non-factorizable corrections are the diagrams (a), (b), and (c) in Fig. |3.1| . If 
the additional exchanged particle is massive, the corresponding correction has no double 
pole for on-shell W bosons. However, if a photon is exchanged between the different sub- 
processes, this leads to a doubly-resonant contribution originating from the soft-photon 
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region. This can be directly seen from the usual soft-photon approximation (SPA), which 
yields contributions proportional to the (doubly-resonant) lowest-order contribution. 
The doubly-resonant contributions can be extracted on the basis of a simple power- 



counting argument. For instance, the loop integral corresponding to diagram (c) in Fig. |3J 
is of the following form: 

H 4 N(q, kj) 

V - A 2 )[(g - k 3 ) 2 - ml][(q - kJ) 2 - MU(Q + K) 2 ~ M&][(q + h) 2 - m| 

d 4 N(q, h) 

q {q 2 ~ A 2 )(g 2 - 2qk 3 ){q 2 - 2qk_ + k 2 _ - M^){q 2 + 2qk + + k\ - M^){q 2 + 2qk 2 ) ' 

(3.1.9) 

where we have introduced an infinitesimal photon mass A to regularize the infrared (IR) 
singularity. The function N(q, kj) involves the numerator of the Feynman integral, i.e. a 
polynomial in the momenta q and fcj, and possible further denominator factors originating 
from propagators (hidden in the blob of the diagrams) that are regular for q = and 
k\ = Myj. For on-shell W bosons (k± = M^), the integral has a quadratic IR singularity. 
For off-shell W bosons, part of the IR singularity is regularized by the off-shellness, k 2 ± — 
M^j ^ 0, such that the usual logarithmic IR singularity remains. Vice versa, the off-shell 
result develops a pole if either W boson becomes on shell, and is thus doubly-resonant. 
Therefore, the quadratic IR singularity in the on-shell limit is characteristic of the doubly- 
resonant non-factorizable contributions. All terms that involve a factor q in the numerator 
are less IR-singular and therefore do not lead to doubly-resonant contributions and can 
be omitted. Similarly, q can be neglected in all denominator factors included in N(q, kj). 
In summary, q can be put to zero in N(q, kj) in DPA. We have checked this for various 
examples explicitly. As a consequence, we are left with only scalar integrals, and the 
non-factorizable virtual corrections are proportional to the lowest-order matrix element. 
We call the resulting approximation extended soft-photon approximation (ESPA). It differs 
from the usual SPA only by the fact that q is not neglected in the resonant W propagators. 



In ESPA, diagram (c) in Fig. |3.1| gives the following contribution to the matrix element: 
2 „ „ . , f d 4 g 4:k 2 k 3 



M hh = ie 2 Q 2 Q 3 M B omJ 



(2vr) 4 (q 2 - A 2 )(g 2 - 2qk 3 )(q 2 + 2qk 2 ) 

(k 2 _ — M 2 v )(fc 2 — Myj) 
[ ( q _ k _ ) 2_ M 2 v][(q + k+) 2_ M ^y I • • ) 



The q 2 terms in the last four denominators are not relevant in the soft-photon limit and 
were omitted in Refs. [57, |3B], |3S|. In fact, using the above power-counting argument it 



can easily be seen that the differences of doubly-resonant contributions with and without 
these q 2 terms are non-doubly- resonant. We have chosen to keep the q 2 terms, because we 



want to use the standard techniques for the evaluation of virtual scalar integrals [43]. In 
DPA, i.e. if we perform the limit k\ — > after evaluating the integral, we should obtain 
the same result. 

In order to arrive at physical results, we have to incorporate the finite width of the 
W bosons. In DPA this can be done in at least two different ways: 
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As a first possibility, we perform the integrals for zero width and afterwards put k± = 
M w where this does not give rise to singularities. In all other places, i.e. in the resonant 
propagators and in logarithms of the form ln(k± — M w + ie), we replace k± — M w + ie by 
K± = k± — M w + iMwTw- Since the width is only relevant in the on-shell limit, it is clear 
that the (physical) on-shell width has to be used. 

Alternatively, we introduce the width in the W propagators before integration. This 
has to be done with caution. If we introduce the finite width by resumming W-self-energy 
insertions, the width depends on the invariant mass of the W boson and thus on the 
integration momentum. Fortunately, the contribution we are interested in results only 
from the soft-photonic region where the virtual W bosons are almost on shell. Therefore, 
we can insert the on-shell width inside the loop integral. After performing the integral, 
we put k± = M w and IV = where this does not lead to singularities. In DPA this gives 
the same results as the above treatment. 

In the following we write M 2 = M w — iMwIV instead of M w in the loop integrals. It 
is always understood that M 2 and k± are replaced by M w where possible after evaluation 
of the integrals. 

If we implement the width into the integrand, it is clear that only the part of the inte- 
gration region with |g | ?S T\v contributes in DPA. If | q Q | ^> T w , one of the W propagators 
must be non-resonant and the contribution becomes negligible. 

Once the width is introduced, it becomes evident that the relative error of the DPA is of 
the order of Tw/scale. Let Eqm be the center-of-mass (CM) energy and AE = Ecu. — 2Mw 
be the available kinetic energy of the W bosons. Then, for AE > Mw the scale is given 
by Mw, for Tw ^5 AE < Mw it is given by AE and for AE < Tw it is given by Tw- This 
shows that the DPA is only sensible several Tw's above threshold. This is simply due to 
the fact that, close to threshold, the phase space where both W propagators can become 
doubly- resonant is very small, and the singly-resonant diagrams become important. 

3.1.3 Classification and gauge-independent definition of the non- 
factorizable doubly-resonant virtual corrections 

Manifestly non-factorizable corrections arise from photon exchange between the final states 
of the two W bosons (ff'); between initial and final state (if), and between one of the 
intermediate resonant W bosons and the final state of the other W boson (mf'). Examples 
for these types of corrections are shown in Fig. [D] (c), Fig. ^]2| (a), and Fig. (a,b), 
respectively. In addition, there are diagrams where the photon does not couple to uniquely 
distinguishable subprocesses. These contributions can be classified into photon-exchange 
contributions between one of the intermediate resonant W bosons and the final state of the 
same W boson (mf), between the intermediate and the initial state (im), between the two 
intermediate W bosons (mm'), and within a single W-boson line, i.e. the photonic part of 
the W-self-energy corrections (mm). Diagrams contributing to these types of corrections 
are given in Fig. |T2] (c), Fig. |372| (b), Fig. [3J] (d), and Fig. |3]2| (d), respectively. Because 
the photon coupling to the W boson can be attributed to the decay or the production 
subprocesses, these diagrams involve both factorizable and non-factorizable corrections. 

In order to define the non-factorizable corrections, we have to specify how the factoriz- 
able contributions are split off. This should be done in such a way that the non-factorizable 
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Figure 3.2: Further examples of non-factorizable photonic corrections in 0(a). 



corrections become gauge-independent. In Refs. |38|, |39| this was reached by exploiting the 
fact that in ESPA the matrix element can be viewed as a product of the lowest-order ma- 
trix element with two conserved currents. Taking all interferences between the positively 
and the negatively charged currents arising from the outgoing W bosons and fermions 
gives a gauge-independent result. 

We have chosen a different definition of the non-factorizable corrections, which, how- 
ever, turns out to be equivalent to the one of Refs. |38|, [39j in DPA. Our approach has the 
advantage of providing a clear procedure how to combine factorizable and non-factorizable 
contributions to the full 0(a) correction in DPA. Because the complete matrix element is 
gauge-independent order by order, the sum of all doubly-resonant 0(a) corrections must 
be gauge-independent. On the other hand, the factorizable doubly-resonant corrections 
can be defined by the product of gauge-independent on-shell matrix elements for W-pair 
production and W decays and the (transverse parts of the) W propagators, 



M f 



1 



SM e+e ~ 



+W+W" 



M 



W+^/l/2^W-^/ 3 /4 



Born 



Born 



(3.1.11) 



+ M e e 



>w+w- 



Born 



Born 



where SAi e+e ^ w+w ; <5A / i w+ ^^ 1 ^ 2 , and SAi w -*M* denote the one-loop amplitudes of the 
respective subprocesses. We can define the non-factorizable doubly-resonant corrections 
by subtracting the factorizable doubly-resonant corrections from the complete doubly- 
resonant corrections. This definition allows us to calculate the complete doubly-resonant 
corrections by simply adding the factorizable corrections, defined via the on-shell matrix 
elements, to our results. Our definition can be applied diagram by diagram. In this way, 
all diagrams that are neither manifestly factorizable nor manifestly non-factorizable can 
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be split. Such diagrams receive doubly-resonant contributions from the complete range 
of the photon momentum q, and not only from the soft-photon region. This is obviously 
due to the presence of two explicit resonant propagators. However, after subtracting 
the factorizable contributions, all doubly-resonant terms that are not IR-singular in the 
on-shell limit cancel exactly, i.e. only the soft-photon region contributes. Consequently, 
also in this case q can be neglected everywhere except for the denominators that become 
IR-singular in the on-shell limit. As an example, we give the no n- factorizable correction 
originating from diagram (d) of Fig. |3.1| f| 



M 



w+w- 

nf 



ie 2 M 



Born 



d 4 C 



4k, k. 



(2tt) 4 q 2 [(q + k + ) 2 - M 2 ] [(q - k_) 2 - M 2 ] 
d A q 4k + k^ 
(2tt) 4 {q 2 - X 2 )(q 2 + 2qk + )(q 2 - 2qk„) 



(3.1.12) 



This example shows that the on-shell subtraction introduces additional IR singularities. 
If the IR singularities in the non-factorizable real corrections are regularized in the same 
way, they cancel in the sum. In (|3.1.12| ) an infinitesimal photon mass A is used as IR 
regulator, but we have repeated the same calculation also by using a finite W-decay width 
as IR regulator instead of A, leading to the same results in the sum of virtual and real 
photonic corrections. 

We illustrate our definition of the non-factorizable corrections also for the photonic 
contribution to the W-self-energy correction [diagram (d) of Fig. |3.2| |. The non-factorizable 
part of the W + self-energy reads 



M n{ ~ 



-ie 2 A<Born|y 



d 4 g 



4k 2 



(2tt) 4 q 2 [{q + k + ) 2 - M 2 ]{k 2 + - M 2 ) 



1 



+ 



k\ — M 2 
d 4 q 



d A q 



4k\ 



(2tt)V(<? 2 + 2# h 

4kl_ 



kl=M 2 



(2vr) 4 (g 2 - X 2 )(q 2 + 2qk^ 



k\=M^ 



(3.1.13) 



The first integral results from the off-shell self-energy diagram, the second from the corre- 
sponding mass-renormalization term, and the third integral is the negative of the on-shell 
limit of the first two integrals. The integrals in ( |3.1.13j ) are UV-divergent and can be easily 
evaluated in dimensional regularization. 

The gauge independence of the non-factorizable corrections has been ensured by con- 
struction. The consistent evaluation of gauge theories requires, besides gauge indepen- 
dence of the physical matrix elements, the validity of Ward identities. It was found in 
Ref. J7|, H B [TO] ^ na ^ ^ ne violation of Ward identities can lead to completely wrong pre- 
dictions. The procedure described above for extracting the non-factorizable corrections 
from the full matrix element does not lead to problems with the Ward identities that 
rule the gauge cancellations inside matrix elements. This is due to the fact that the non- 
factorizable corrections are proportional to the Born matrix element. Therefore, if Ward 



1 We use the sign ~ to indicate an equality within DPA, i.e. up to non-doubly-resonant terms. 
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identities and gauge cancellations are under control in lowest order, the same is true for 
the non-factorizable corrections. 

Finally, we show how our definition of the non-factorizable corrections can be rephrased 
in terms of products of appropriately defined currents. By using 



(q ± k±) 2 — M 2 q 2 ± 2qk± 



k\ 



M 2 



(q±k ± ) 



M 2 



(3.1.14) 



and the fact that in DPA k\ can be put to M(y before integration in integrals that do not 
depend on M 2 , the contribution ( |3.1.12|) can be expressed as 
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(q-k_) 2 -M 2 (q + k + ) 2 - M 2 (q - k_) 2 - M 2 



. (3.1.15) 



The other non-factorizable corrections that involve photons coupled to W bosons can 
be rewritten in a similar way. Finally, all non-factorizable virtual corrections can be cast 
into the following form: 
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(3.1.16) 
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2k 



4 M 
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2k. 



q 2 - 2qk 4 q 2 - 2qk_ J (Jfe_ - q) 2 - M 2 



k 2 — M 2 



: . (3.1.17) 



The last term in ( |3.1.16| ) originates from the Feynman graphs shown in Fig. |3.1| and 
those where the final-state fermions are appropriately interchanged [interference terms 
(ff), (mf), and (mm')]. The contributions involving the current j^ ir % „ _>w w contain the 
interference terms (if), (mf), (im), (mm), and the remaining contributions of (mf') and 
(mm'). The contribution of the W + self-energy is given, for instance, by the product of 

the two terms involving k +fl in j v tt,^T W+W an d Jvirt^^^ ' ■ 

In DPA, the q 2 terms in the denominators of (|3.1.17|) can be neglected, and the cur- 
rents are conserved. The currents jXt./T^ 2 an d jXt./T^ 4 are the ones mentioned in 
Refs. |38| , |39| . For the virtual corrections, this shows that our definition of non-factorizable 
doubly-resonant corrections coincides with the one of Refs. [^, 3I| in DPA. 
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Figure 3.3: Example of a non-factorizable real correction. 



3.1.4 Doubly-resonant real corrections 

The photonic virtual corrections discussed above are IR-singular and have to be combined 
with the corresponding real corrections in order to arrive at a sensible physical result. The 
real corrections originate from the process 

e+(p + )+e-(p_) -> W + (A4)+W- (£;'_)[ + 7(g)] 

- fi(K) + HK) + fs(k' 3 ) + HK) + 7 (g). (3.1.18) 

Note that we have marked the fermion momenta k\ by primes in order to distinguish 
them from the respective momenta without real photon emission. The momenta of the 
W bosons are k' + = k[ + k' 2 and k'_ = k' 3 + k' 4 if the photon is emitted in the initial state 
or in the final state of the other W boson, and k + = k[ + k' 2 + q or k- = k' 3 + k' 4 + q if the 
photon is emitted in the final state of the respective W boson. 

The non-factorizable corrections induced by the process ( |3.1.18| ) arise from interferences 
between diagrams where the photon is emitted from different subprocesses. A typical non- 
factorizable contribution is shown in Fig. 1373]. Including the integration over the photon 



phase space, this contribution has the following form: 

•d 3 q N(q, k[)5( P+ + p_- k[ - k' 2 - k' 3 - k\ - q) 
2g [(q + k' 2 ) 2 - m 2 2 ][(q + k' + ) 2 - M 2 ](k' 2 - M 2 ) 
{ 1 

x 



X 



[(g + k' 3 ) 2 - ml]{k% - M 2 )[(g + k'_) 2 - M 2 } 
d 3 q N(q, K)5( P+ +p.-k' 1 -k' 2 ~k' 3 -k' 4 -q) 
2g (2qk' 2 )(2qk' + + k' 2 - M 2 )(k' 2 - M 2 ) 

1 



x 



{2qk' 3 )[k' 2 - (M*) 2 ][2qk'_ + k' 2 - (M*) 



(3.1.19) 

90=\/q 2 +A 2 



where we again use a photon mass A to regularize the IR singularities, and N(q, k^) has 
the same meaning as above. Again, the doubly-resonant contributions are characterized 
by a quadratic IR singularity for k'± = M{y, Tw — » 0, and only soft-photon emission is 
relevant in DPA. For this reason, the W bosons are nearly on shell, and the on-shell width 
is appropriate. As for the virtual corrections, the introduction of the width before or after 
phase-space integration leads to the same results in DPA. As already indicated in ( |3.1.19| ), 
in the following real integrals we use M 2 with the understanding that it has to be replaced 
by M-yy after integration where possible. 
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The aim is to integrate over the photon momentum analytically and to relate the 
fermion momenta k[ to the ones of the process without photon emission, k^. Primed and 
unprimed momenta differ by terms of the order of the photon momentum: k[ = ki + 0(qo). 
In DPA we can neglect q in N(q,k' i ), leading to the replacement Ni^q.k'^ — > N(0,ki). 
Moreover, we can extend the integration region for go to infinity, because large photon 
momenta yield negligible contributions in DPA. After extension of the integration region 
the integral becomes Lorentz-invariant. 

While the correction factor to the lowest-order cross section is universal in SPA for 
all observables, the correction factor is non-universal in ESPA. In order to define this 
correction factor in a unique way, one has to specify the parameterization of phase space, 
i.e. the variables that are kept fixed when the photon momentum is integrated over. This 
fact has not been addressed in the literature so far. 

Let us consider this problem in more detail. It can be traced back to the appearance 
of the photon momentum q in the 5-function for momentum conservation. In the usual 
SPA q is neglected in this 5-function, which is sensible if the exact matrix element is a 
slowly varying function of q in the vicinity of q = 0. However, in the presence of resonant 
propagators, in which q cannot be neglected, the simple omission of q in the momentum- 
conservation 5-function leads to ambiguous results: putting q = in the 5-function and 
identifying k[ with ki in ( |3.1.19|) yields 



1 



d 3 q N{0, kj)5(p+ + p--k 1 -k 2 -k 3 -k 4 
2q (2qk 2 ) (2qk + + k\ - M 2 )(k 2 _ - M 2 ) 
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x 



(2qk 3 )[k 2 + - (M*) 2 ][2qk. + ki - (M*) 2 ] 



, • (3-1-20) 



On the other hand, eliminating k' + in the denominator of (|3.1.19| ) with the help of the 
5-function, putting q = in the 5-function, restoring k' + with the modified 5-function, and 
setting k[ — > ki results in 



1 



d 3 q iV(0, kj)5(p+ + P- ~ h- k 2 - fc 3 - fc 4 ) 
2q (2qk 2 ) (k 2 + - M 2 ) (k 2 _ - M 2 ) 

1 

x 



{2qkz)[-2qk + + k% - (M*) 2 ][2qk_ + k 2 _ - (M*) 2 ] 



. - (3-1.21) 

90=Vl 2 + A2 



Both expressions differ by a doubly-resonant contribution. The difference is in general 
confined to the W propagators and originates from the fact that not only soft photons but 
also photons with energies of the order of \k± — M^\/M^ or, after the inclusion of the 
finite width, of order Tw contribute in DPA. Since photons with finite energies contribute, 
it is evident that the integral over the photon momentum depends on the choice of the 
phase-space variables that are kept fixed. 

As a consequence, one has to choose a definite parameterization of phase space and to 
exploit the ^-function carefully, in order to define the non-factorizable corrections uniquely. 
For instance, if the vector k' + = k[ + k' 2 is kept fixed, the alternative ( 3.1.21| ) is excluded. 



However, because of momentum conservation, not all external momenta can be kept fixed 
independently. 
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It is, however, possible to keep, for instance, the invariant masses of the final-state 
fermion pairs k' 2 = (k[ + k' 2 ) 2 and k' 2 = (k' 3 + k' A ) 2 fixed when integrating over the photon 
momentum. If we require {k^ + k^ 2 = k' 2 = k'] r = (ki + k 2 ) 2 , we obtain for the denominator 
of the W + boson 

(q + k' + ) 2 — M 2 = 2qk' + + k' 2 — M 2 = 2qk + + k\ - M 2 + O{q 2 ) 

= ( q + k + ) 2 -M 2 + 0(q 2 ), (3.1.22) 

where k\ = ki + O(qo) was used. Based on the power-counting argument given above, the 
terms of order q$ can be neglected in DPA, and we find 

{q + k'f-M 2 ~ (q + k + f-M 2 . (3.1.23) 



If we choose to eliminate k' +} as done in the derivation of ( |3.1.21| ), we find, on the other 
hand, 

(q + k' + ) 2 — M 2 = (p + + p_ - k'_) 2 — M 2 = (p + + p_ - k_) 2 — M 2 + O(q ) 

= k\- M 2 + O(q ). (3.1.24) 

The O(qo) terms are relevant in DPA [and in fact given by (|3.1.22j) l. As a consequence, 
( p.l.21| ) is not correct if we choose to fix k' 2 = k\ when integrating over the photon momen- 
tum. For fixed k! 2 = h\, ( |3.1.22j ) leads to the unique result ( |3.1.20|) for the W + propagator 



in DPA, independently of the other phase-space parameters. If we choose, on the other 
hand, to fix fci = k+, which corresponds to a different definition of the invariant mass of 
the W + boson, we obtain 

(q + k' + ) 2 — M 2 = k\ — M 2 , (3.1.25) 
and thus (|3.1.21|) instead of ( |3.1.20| ). Consequently, the different approaches (|3.1.20|) and 



( p. 1 . 21|) correspond to different definitions of the invariant mass of the W + boson which 
decays into the fermion pair (/l,/^)- In order to define the DPA for real radiation, one 
has to specify at least the definition of the invariant masses of the W bosons that are kept 
fixed. In the following we always fix k\ = (k[ + k' 2 ) 2 = (k\ + k 2 ) 2 and k 2 _ = (k' 3 + k' 4 ) 2 = 
(&3 + ki) 2 , as it was also implicitly done in Refs. |37], ^8|, ^9[. Once the invariant masses 
of the W bosons are fixed in this way, the resulting formulae for the non-factorizable 
corrections hold independently of the choice of all other phase-space variables. 

We stress that the results obtained within this parameterization of phase space dif- 
fer from those in other parameterizations by doubly-resonant corrections. As already 
indicated in the introduction, in an experimentally more realistic approach the invariant 
masses of W bosons are identified with invariant masses of jet pairs, which also include 
part of the photon radiation. Since this situation can only be described with Monte Carlo 
programs, our results (as well as those of Refs. |37], |38|, |39|) should be regarded as an 
estimate of the non-factorizable corrections. 



3.1.5 Classification and gauge-independent definition of the non- 
factorizable doubly-resonant real corrections 

The doubly-resonant real corrections can be classified in exactly the same way as the 
virtual corrections. For each virtual diagram there is exactly a real contribution, which 
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Figure 3.4: Real bremsstrahlung diagram containing non-factorizable and factorizable 
contributions. 



we denote in the same way, e.g. ff refers to all interferences where the photon is emitted 
by two fermions corresponding to the two different W bosons. 

As in the case of the virtual corrections, one has to define the non-factorizable real 
corrections in a gauge-independent way. To this end, we proceed analogously and define 
the non-factorizable real corrections as the difference of the complete real corrections and 
the factorizable real corrections in DPA. The factorizable corrections are defined by the 
products of the matrix elements for on-shell W-pair production and decay with additional 
photon emission. There are three contributions to the factorizable real corrections, one 
where the photon is radiated during the production process and two where it is emitted 
during one of the W-boson decays. The corresponding matrix elements read 
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M 2 } 



(3.1.26) 



in analogy to ( |3.1.8|) . Note that the matrix elements A^ rea i,2 an d A^ rea i,3 involve an explicit 
g-dependent propagator. The factorizable corrections are given by the squares of these 
three matrix elements, and include by definition no interferences between them. 

As an example for the extraction of non-factorizable corrections from real diagrams that 
involve both factorizable and non-factorizable corrections, we consider the contribution of 



the diagram in Fig. [3.4| . After subtraction of the factorizable contribution, which originates 
2 it gives rise to the following correction factor to the square |A^Bom| 2 of the 
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lowest-order matrix element ( p. 1.8 
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(3.1.27) 
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Note that the form of the correction factor is only correct for fixed {k^ + k^) 2 and (k^ + k^) 2 . 
For other conventions the off-shell contribution changes, whereas the on-shell contribution 
stays the same. Using the relations ( |3.1.14j ) for q 2 = 0, in DPA we can rewrite ( |3.1.27| ) as 
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real 



~ — e 
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(3.1.28) 



In the same way, all other contributions that originate from a photon coupled to a 
W boson can be rewritten such that the complete real non-factorizable corrections can 
finally be expressed as the following correction factor to the lowest-order cross section, 
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(3.1.29) 



with the currents 



•e+e-^W+W- 
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(3.1.30) 



The factor (|3.1.29|) for the non-factorizable correction can be viewed as the interference 
contributions in the square of the matrix element (e denotes the polarization vector of the 
photon) 

-/^'real — - /v 'Born£ Urealu + Jreal,u + Jreal.u J; (^.J-.O-LJ 



which is just the sum of the three matrix elements M. Tea x,i, i = 1,2,3, in ESPA (includ- 
ing radiation from the external fermions and the internal W bosons). The respective 
squares of these three contributions correspond to the factorizable corrections. Note that 

is the soft-photon matrix element for on-shell W-pair pro- 
w : -/i/ 2£ MjW+-/i/2 and Mvorn hh ^J^ hh correspond to the 
soft-photon matrix elements for the decays of the on-shell W bosons, apart from the extra 
factors (A; 2 . — M 2 )/[(k± + q) 2 — M 2 ]. These factors result from the definition of the lowest- 
order matrix element in terms of k± and from the fact that we impose k: 



•' l Born fc Jreal,fj, 

duction. Similarly, M%Z?' l '*e*j££»» and M^ 3 ^^,, 



kl 



k\. 



,12 



k± and not 



Obviously, the currents ( |3.1.30|) are conserved, and the corresponding Ward identity 
for the U(l) em symmetry of the emitted photon is fulfilled. 
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3.2 Calculation of scalar integrals 



In this section we set our conventions for the scalar integrals. We describe the reduction 
of the virtual and real 5-point functions to 4-point functions and indicate how the scalar 
integrals were evaluated. More details and the explicit results for the scalar integrals can 
be found in the appendix. 



3.2.1 Reduction of 5-point functions 
Reduction of virtual 5-point functions 

The virtual 2-, 3-, 4-, and 5-point functions are defined as 

Bo(pi,rn ,mi 

Co(pi,P2,m ,m 1 ,m 2 

D {o,n} (pi ,P2,P3,m ,m 1 ,m 2 ,m 3 

E{o,rf (pi, P2, Pa, Pi, m , mi, ra 2 , m 3 , m 4 
with the denominator factors 



1 /> 1 



17T / 
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17T 
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NoNiN 2 n 3 ' 

iVoiViAWaAV 



jV = g 2 - m„ + ie, Ni = (q + pi) 2 - m 2 + ie, 



i = 1,...,4, 



(3.2.i; 



(3.2.2) 



where ie (e > 0) denotes an infinitesimal imaginary part. 

The reduction of the virtual 5-point function to 4-point functions is based on the fact 
that in four dimensions the integration momentum depends linearly on the four external 
momenta pi [13|, |44| . This gives rise to the identity 







2q 2 2qpi . . . 2qp 4 
2p x q 2p\ ... 2p x p 4 



2p A q 2p 4 p x 



2pl 



2N + Y 00 2q Pl ... 2qp 4 

m-No + Y^-Yoo 2p\ ... 2p lP4 

N 4 -N + Y i0 - Y 00 2p iPl ... 2p\ 



(3.2.3) 



with 



Yoo = 2m , Y i0 = Y 0i = m Q + m i -p t , Y {j = m i + m j - (pi - pj) , i, j 
Dividing this equation by N Ni • ■ ■ N 4 and integrating over d A q yields 



1,2,3,4. 
(3.2.4) 



= 7-5 A 
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(3.2.5) 
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Expanding the determinant along the first column, we obtain 



2A)(0)+r o£c 



2piPi . . . 2p x p A 



2p A p x . . . 2p A p A 

+ J2(-^ k {D„(k)-[D ll (0)+p All D (0)} + p 4fl D (0) + (Y k0 -Y 00 )E fl } 
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2p k _ipi . . . 2p k _ip A , (3.2.6) 
tyk+iPi • • • 2Pk+m 

2p A p x . . . 2p A p A 

where D (k) denotes the 4-point function that is obtained from the 5-point function E 
by omitting the fcth propagator . The terms involving p A)M D (0) have been added for 
later convenience. 

All integrals in ( |3.2.6| ) are UV-finite and Lorentz-covariant. Therefore, the vector 
integrals possess the following decompositions 



En = EiPin, 

i=l 

1=1 
3 

Dp{0) +p^D o {0) = E A(0)(Pi-p4V 



k = 1,2,3,4, 



(3.2.7) 



i=i 



The last decomposition becomes obvious after performing a shift q — > q— p A in the integral. 
From ( |3.2.7[ ) it follows immediately that the terms in ( |3.2.6D that involve D^{k) drop 
out when multiplied with the determinants, because the resulting determinants vanish. 
Similarly, -D^(O) + p A ^D (0) vanishes after summation over k. Finally, the term p AfM D (0) 
contributes only for k = 4, where it can be combined with the first term in Q3.2.6| ). 
Rewriting the resulting equation as a determinant and reinserting the explicit form of the 
tensor integrals leads to 



= 7-5 A 
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(3.2.8) 
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Using 



2qpj 



y%j ~ Yio ~ Yoj + *oo> 
Nj - iV + Y 0j - Yqo, 



(3.2.9) 



adding the first column to each of the other columns, and then enlarging the determinant 
by one column and one row, this can be written as 
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Equation ( |3.2.10| ) is equivalent to 
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(3.2.11) 



which expresses the scalar 5-point function Eq in terms of five scalar 4-point functions 



(3.2.12) 



where Y = (Yy), and Yi is obtained from Y by replacing all entries in the ith column 
with 1. 

In the special case of an infrared singular 5-point function we have 



^oo — 2A — > 0, Y01 
and the determinants fulfill the relations 

= det(y) - Y 02 det(Y 2 ) - Y 03 det(Y 3 ). 
= (F 24 - Y Q2 ) det(Y) + Y 02 Y U det(Yi) 
= (F31 - Y 03 ) det(Y) + Fossil det(F 1 ) 



m i - v\ = 0, 



Y 



04 



ml - v\ = 0, 



(3.2.13) 



(3.2.14) 

(Y w Yu - Y 03 Y 24 ) det(Y 3 ) + Y 02 Y U det(Y 4 ), 
(Y 03 Y 21 - Y 02 Y 31 ) det(Y 2 ) + Y 03 Y 41 det(Y 4 ), 



which allow the simplification of ((3.2.12 ). 
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Reduction of real 5-point functions 



The real 3-, 4-, and 5-point functions are denned as 



Co(pi,p 2 , \,m 1 ,m 2 ) 
^o(Pi,P2,P3, \,mi,m 2 ,m 3 ) = 3 
^o(pi,P2,P3,P4, \,m 1 ,m 2 ,m 3 ,m 4 ) 



d 3 



q 



qo <AE2q N[N^ 
d 3 q 1 



2g N[N^ 
d 3 q 1 



) 

9o=vV+ A2 



tt7 2g NIN2N3N4 



qo- 



with 



iV/ = 2gpj + p 2 — m 



(3.2.15) 



(3.2.16) 



The shift of the integration boundary of go to infinity leads to an artificial UV divergence 
in the 3-point function Cq, which is regularized by an energy cutoff AE —>■ 00. In the 
following we only need differences of 3-point functions that are independent of AE and 
Lorentz- invar iant . 

Because of the appearance of UV-singular integrals in intermediate steps, the reasoning 
of the previous section cannot directly be applied to Sq. Therefore, we rewrite the real 
5-point function as an integral over a closed anticlockwise contour C in the go plane and 
introduce a Lorentz-invariant UV regulator A: 

If 1 —A 2 

S (p 1 ,p 2 ,p 3 ,p 4 ,\,m 1 ,m 2 ,m 3 ,m 4 ) = lim — / d 4 g — — — — — (3.2.17) 

A^oo m 2 Jc NqN[ ■ ■ ■ N4 q 2 — A 2 

with 

Nq = q 2 — A 2 . (3.2.18) 

The contour C is chosen such that it includes the poles at g = Vq 2 + A 2 and go = 
Vq 2 + A 2 , but none else (see Fig. |3~5Dp|. 

The integral ( 3.2.17Q can be reduced similarly to the virtual 5-point function. Owing 
to the different propagators ( |3.2.3j ) leads to 



A lim ~ 

A— >oo m z 



d 4 g 



1 



-A 2 



N^N{ ■ ■ ■ Ni g 2 - A 2 



2N' Q 



2qp 1 
N[ + Y 10 2p\ 

K + ^40 2p 4 Pl 



2gp 4 
2pip 4 

2pl 



(3.2.19) 



instead of ( 3.2.5 ), with Yij from ( |3.2.4| ), and A 2 can be set to zero in all Yij, in particular 
we have loo = 2 A 2 — > 0. After expanding the determinant along the first column and 



2 It is straightforward to check that the naive power counting for the UV behaviour in ( 3.2.17 ) is valid 
for time-like momenta pi ; light-like pi can be treated as a limiting case. The basic idea for the proof is 
to deform the contour C to the vertical line Rejgo} = | q | — e with a small e > 0, which is allowed for 
sufficiently large |q|, more precisely when all particle poles appear left from the line Re{go} = |q| — e - 
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Im{<jo} 







90 (P4) | X ^ 


^\ c 






[ 1 / 90 (Pi) 




• 


• • 


\+* — 1 


1 * J *" 


-\J q 2 +A 2 


• 

<?0(P2) 


|Vq 2 +A 2 


v/q 2 +A 2 



Re{g }=|q|- e 



Figure 3.5: Illustration of the contour C of ( |3.2.17| ) in the complex go plane. The open 
circles indicate the "particle poles" located at qo(pi) = (2qpi — pj + m 2 ) / (2p i0 ) . 



using the Lorentz decompositions of the integrals where N' k is cancelled, we see that these 
terms vanish, and we are left with 



If 1 

0= * lin L —2 L d< iTuw- 



-A 2 



A^oo ivr 2 Jc * N^N[ ■ ■ ■ N' A q 2 - A 2 



2N' Q 2q Pl 
Yio 2pi 

Y m 2p 4 pi 



2qp A 



2p 2 4 



(3.2.20) 



Using 



ZpiPj — Yij — — Yoj, 
2q Pj = N; + Y 0j , 



(3.2.21) 

adding the first column to the other columns and extending the determinant leads to 



If 1 
= lim — / d 4 q — Ar> 

A-oo ITT 2 Jc N^N[ ■ 



-A 2 



• N' iq 2 - A 2 



1 Y 01 

2N' Q N[ + 2N' Q + Foi 

Y w Yxx — Y i 

Y40 Y41 — Yq\ 



^04 

N' 4 + 2N' G + F 04 

Y L4 — 1q4 



I44 — 1q4 



(3.2.22) 
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Subtracting the first row from the second, adding the first row to the other rows, and 
exchanging the first two rows, we arrive at 



= lim 



1 



-A 2 



a - x , ^ Jc dq N^N[ ■ ■ ■ N' A q 2 - A 2 



-1 2N' Q N[ + 2N 
1 Koi 
1 no Y n 



1 Ya 



10 



Y A 



K + 2N 

Yqa 

Y u 



Ya 



14 



(3.2.23) 



Now we perform the contour integral using power counting for A — > oo (see footnote ||]). In 
the contribution of the pole at go — Vq 2 + -a 2 we have Nq = in the numerator; the limit 
A — > oo can be trivially taken and reproduces usual bremsstrahlung integrals, as defined 
in ( |3.2. 15[ ). In the contribution of the pole at go — Vq 2 + A 2 the term containing Nq in the 
numerator survives and will be calculated below, but all other terms vanish after taking 
the limit A — ► oo. 
Thus, we find 

-So P (0) V (1) + V Q (0) V Q (2)+f> (0) P o (3)+P o (0) V (A)+V (0) 



1 


Yqo 


ni 


Y02 


Y03 


Yq4 


1 


Y l0 


Y n 


Y 12 


Y 13 


Y U 


1 


Y20 


Y 21 


Y22 


Y23 


Y24 


1 


no 


ni 


Y32 


n 3 


Y u 


1 


no 


ni 


Y42 


Y43 


Y44 



or 



£0 



det(y) 



det(Fo) V (0) + £ det (n) \v (i) + V (Q) 



i=l 



Here 



Z> (0) = lim 



d 4 g 



-A 2 



(3.2.24) 
(3.2.25) 

(3.2.26) 



A^oo m z J C ' N{ ■ ■ ■ N' 4 q 2 - A 2 ' 
and the 4-point bremsstrahlung integrals T>o(i), i = 1,2,3,4, result from So by omitting 
the ith denominator N[. The result ( |3.2.24|) differs from ( 3.2.10Q only by the extra Pq(0)' s 
added to the Po(0' s - 

The integral X>o(0) stems from the terms involving in the numerator in ( |3.2.23|) and 
can be expressed as follows: 



£0(0) = — 2 [ dV 1 a 1 

m z Jc 2q'pi ■ ■ ■ 2q'p 4: q' z — 1 



(3.2.27) 



where the contour C surrounds q' Q 
yields 



q' +1. Performing the contour integral over dq' Q 



f> (0) 



d 3 q' 1 
2q' 2q'p 1 ■ ■ ■ 2q'p 4 ^ 



(3.2.28) 
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Now the vector q' is time-like. Since also the vectors Pi are time-like (or at least light-like), 
the scalar products q'pi cannot become zero. After redefining the momenta, 



Pi = &iPi 



Oi = ±1, 



Pio > 0, 



(3.2.29) 



and extracting the signs Oi, this integral can be evaluated by a Feynman-parameter rep- 
resentation and momentum integration in polar coordinates resulting in: 



X\>(0) = — o\020^a^ \ dxi dx 2 dx 3 dx 4 5 I 1 — 



i=l 



-i -2 



^2%iP, 



l x i=l 



. (3.2.30) 



This is just the Feynman-parameter representation of a virtual 4-point function such that 
we finally obtain 



(3.2.31) 



£>o(0) = -oi a 2 a 3 a±D [p 2 -pi,Pa ~ Pi,Pi ~ Pi, VPi, VPl, VPs, VP* ■ 



Explicit reduction of the virtual 5-point function for the photon exchange 
between f 2 and f 3 

For the photon exchange between f 2 and /3 the following scalar integrals are relevant: 



E 


= Eq{ 


-fa, 


—k-, k + , fa, A, m 3 , M, M, m 2 ), 


A)(0) 


= A>( 


-fa, 


k + + fa, fa + fa, 0, M, M, 0), 


A>(1) 


= A>( 


-k- 


k+,k 2 ,0,M, M,m 2 ), 


D (2) 


= A>( 


-fa, 


k + ,k 2 ,\,m 3 ,M, m 2 ), 


A) (3) 


= D ( 


-fa, 


-k-,k 2 ,X,m 3 ,M, m 2 ), 


-Do (4) 


= A>( 


-fa, 


-k_,k + ,0,m 3 ,M, M). 



(3.2.32) 

Since we neglect the external fermion masses, the last two relations ( p. 2. 14 ) simplify to 



= (s 23 + s 24 ) det(y) + K_s 23 det(Yi) - K + {s 23 + s 24 ) + K_M\ det(K 



3), 



= (sia + s 23 ) det(y) + K + s 23 det(y 4 ) - [K-(s 13 + s 23 ) + K + M 2 _] det(F 2 ). (3.2.33) 
These relations allow us to eliminate det(Yi) and det(>4) from (|3.2.12|) , resulting in: 

det(Yo) 



E (-fa, -k-, k + , fa, X, m 3 , M, M, m 2 ) 



det(y) 



A)(0) 



+ a TF^TT {\ K -i 3 ™ + s ^ + K + M 2 ]D (4) + K + s 23 D (2)\ 
det(F )K + s 23 i J 

-^±^D (4) - ^^A(l). 



(3.2.34) 
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The matrix Y reads 
/ 

Y = 



-K_ 

* M 2 

* * 2M 2 



-K + 
{-K + - Si 3 - s 23 ) 
2M 2 - s 
2M 2 
* 





-S23 

~ s 23 ~ s 2i) 

M 2 




(3.2.35) 



Neglecting terms that do not contribute to the correction factor in DPA, the corresponding 
determinants are given by 

det(F) ~ 2s 23 \K + K.s 1A s 23 - (K + M^ + K„s 13 )(K„M^ + K + s u ) 



-K 



det(F ) . 
det(yi) ~ Kt 



M^(s 23 — s 2i) + (^23 + •524)('Sl3S24 — ^14523) 

+ K_M^(-M^ + 2s 13 s 23 + s 13 s 24 + sus 23 ), 
det(y 2 ) ~ -s 23 [K + (M^ + s 13 s 24 - s u s 23 ) + 2K_M^s 13 

det(y 3 ) = det(y 2 )| x+ ^_ iSi3 _ 24 , 
det(y 4 ) = det(y 1 )| x+ ^_ >si3 _ 24 , 

where the shorthand /t\y is defined in Appendix |A.1[ 



(3.2.36) 



Explicit reduction of the real 5-point function for the photon exchange between 
f 2 and / 3 

In this case the integrals appearing in (|3.2.25|) read 

£ = £ (k 3 ,k_,k + ,k 2 ,\,m 3 ,M*,M,m 2 ), 
V (0) = -D (k 4 , k + - k 3 , k 2 - k 3 , 0, M_, M+, 0), 
P (l) = V (k.,k + ,k 2 ,0,M*,M im2 ), 
V (2) = V {k 3 , k + , k 2 , A, m 3 , M, m 2 ), 
V Q (3) = V (k 3 , k-, k 2 , A, m 3 , M*, m 2 ), 

P (4) = A)(&3, fc-, 0, m 3 , M*,M). (3.2.37) 

In analogy to the virtual 5-point function, we can express the real 5-point function with 
the help of the relations (we denote the matrix Y for the real 5-point function with a 
prime, in order to distinguish it from the one for the virtual 5-point function): 

= (s 23 + s 24 ) det(r) + K*_s 23 det(Y() - [k + (s 23 + s 24 ) - K*_M 2 + ] det(K0, 

= (si 3 + s 23 ) det(r) + K + s 23 det(F 4 ') - [k*_(s 13 + s 23 ) - K + M 2 } det(K,'), (3.2.38) 

by 

det(Y') ~ 

£ (k 3 , fc_, k+, k 2 , A, m 3 , M*, M, m 2 ) = ^-^lv (0) 

52 



+ 



+ 



det(KQ 



det(Y')Kls 23 

det(F 2 ') 
det(Y')K + s 23 

Sl3 + s 23 



K+(s 23 + s 2A ) - KlMl P (l) + A)(0) + K*_s 23 X? (3) + ^o(O) 



Po(4)+P (0) 

A + S 2 3 

For the matrix Y' we find 



iC(si3 + s 23 ) - K+M!] [P (4) + P (0) 

S23 + s 24 



K*_S 



2:-! 



X>o(l)+X> (0) 



F' 







-K* 



(M*) 2 
* 2(M*) 2 



-2K, 



-K+ 
-K + + s 13 + s 23 ) 
2iT + s - (M*) 2 - M 2 ) 

2M 2 
* 



K + s 23 [V (2)+V (0) 

(3.2.39) 





S23 

-JT* + s 23 + s 24 ) 
M 2 

J 



(3.2.40) 



Replacing —K*_ by K- and (M*) 2 by M 2 and multiplying the second and third columns 
and rows by —1, this becomes equal to ( |3.2.35| ) in DPA. 

In DPA, Z> (0) can be neglected in the terms V (i) +V (0) in ( p.2.25|) and ( p.2.39|) , and 
the reduction of the real 5-point function becomes algebraically identical to the reduction 
of the virtual 5-point function, apart from the differences in signs of some momenta. As 
a consequence, the results for the virtual corrections can be translated to the real case if 
we substitute K_ — ► — K*_ in all algebraic factors such as the determinants and i?o — * £0, 

A)(0) -> P o (0), D (l) -> -X> (1), D (2) -> -P (2), ,Do(3) -> X> (3) and ,Do(4) -> Z> (4). 
In particular, the determinants are related by 



det(F') 
det(F/) 
det(F 3 ') 



+ det(F) 
-det(F; 
+ det(F 3 ; 



K- 



K- 



-K* 



-K* 



-K* 



det(Fj) ~ +det(F ) 
det(F 2 ') ~ -det(F 2 ; 
det(F 4 ') ~ +det(F 4 ) 



K- 



K- 



-K* 



-K* 



(3.2.41) 



3.2.2 Calculation of 3- and 4-point functions 



The scalar loop integrals have been evaluated following the methods of Ref. fll3" |. Our 



explicit results are listed in Appendix [A.3| . For vanishing W-boson width they agree with 
the general results of Ref. |43j . For finite W-boson width the virtual 4-point functions are 
in agreement with those of Refs. [38, 39] in DPA. This shows explicitly that the q 2 terms 



in the W-boson and fermion propagators are irrelevant in DPA. 

An evaluation of the bremsstrahlung integrals, which follows closely the techniques for 
calculating loop integrals, is sketched in Appendix [A.2| . The final results in DPA are listed 
in Appendix [A.3| , and the 4-point functions agree with those of Refs. 0, . We have 
analytically reproduced all exact results for the occurring bremsstrahlung 3- and 4-point 
integrals by independent methods. In addition, we have evaluated the IR-finite integrals 
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£>o(l)> £>o(4), an d £o — by a direct Monte Carlo integration over the photon 

momentum, yielding perfect agreement with our exact analytical results for these integrals. 
Note that this, in particular, checks the reduction of the bremsstrahlung 5-point function 
described in the previous section. 

Our results for the 3-point functions cannot directly be compared with those of 
Refs. [38|, , because different approaches have been used. While our results are IR- 
singular owing to the subtracted on-shell integrals, the results of Refs. |5SL are artifi- 
cially UV-singular owing to the neglect of q 2 in the W propagators. However, when adding 
the real and virtual 3-point functions the two results agree. This confirms that our defi- 

in DPA. 



nition of the non-factorizable corrections is equivalent to the one of Refs. [3£ 
Thus, it turns out that in DPA the subtraction of the on-shell contribution is effectively 
equivalent to the neglect of the q 2 terms in all but the photon propagators. 



3.3 Analytic results for the non-factorizable correc- 
tions 

3.3.1 General properties of non-factorizable corrections 



In Ref. [|37| it was shown from the integral representation that the non-factorizable correc- 
tions associated with photon exchange between initial and final state vanish in DPA. This 
was confirmed in Refs. |38], |3~9| . Via explicit evaluation of all integrals we have checked 
that the cancellation between virtual and real integrals takes place diagram by diagram 
once the factorizable contributions are subtracted. In this way all interference terms (if), 
(mf), (im), and (mm) drop out. Examples for the virtual Feynman diagrams contributing 
to these types of corrections are shown in Fig. |3.2| . 

The only non- vanishing non-factorizable corrections are due to the contributions (ff), 
(mf), and (mm'). The corresponding virtual diagrams are shown in Fig. S7I , apart from 



permutations of the final-state fermions. Two of the corresponding real diagrams are 
pictured in Figs. [ITS] and |3.4j. Since these corrections depend only on s-channel invariants, 



the non-factorizable corrections are independent of the production angle of the W bosons, 
as was also pointed out in Refs. [3~F 



3.3.2 Generic form of the correction factor 

The non-factorizable corrections do~ n f to the fully differential lowest-order cross-section 
duBorn resulting from the matrix element ( p. 1.8 ) take the form of a correction factor to 
the lowest-order cross-section: 

dffnf = S n{ ddBorn- (3.3.1) 

Upon splitting the contributions that result from photons coupled to the W bosons ac- 
cording to 1 = Qw+ — Qi ~ Qi an d 1 = —Qw- — Qi ~ Q3 into contributions associated 
with definite final-state fermions, the complete correction factor to the lowest-order cross- 
section can be written as 

<W = E E (-l) a+fe+1 QaQb - Re{A(k + , k a , fc_, k b )}. (3.3.2) 

a=l,2 6=3,4 71 
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In the following, only A = A(k + , k 2 ; k_, k 3 ) is given; the other terms follow by obvious 
substitutions. As discussed above, A gets contributions from intermediate-intermediate 
(A mm /), intermediate-final (A mf /), and final-final (Aff) interactions: 

A = A w + A mf + A f f. (3.3.3) 

The quantity (A mm / is independent of the final-state fermions. The individual contribu- 
tions read 



A r 



(2M& - s){ C (k + , -k.,0, M, M) - C (k + , -k.,\, M w , M w ) 



k\=Ml, 



- C (k + , k_, 0, M, M*) + C (fc+, fc_, A, M w , 



w. 



(3.3.4) 



A mf ' ~ -(s 23 + s 24 ) K+D (l) - K+V (l) - (s 13 + s 23 ) tf_A)(4) - iT P (4) , (3.3.5) 



Aff ~ —K+s 



2.3 



K_£ - iT£ 



(3.3.6) 

with the arguments of the 5- and 4-point functions as defined in ( |3.2.32| ) and ( |3.2.37| ). 

The sum A m f + A f f can be simplified by inserting the decompositions of the 5-point 
functions (gggp and (|3.2.39|). In DPA this leads to 



A m f + Aff 



det(F) uw det(Y') 

- {[K+(s 23 + ^24) + K„Ml\D (l) + K„s 23 D (3)} 

- ^^y^ {[^-(si3 + s 23 ) + AT + M^]A,(4) + K + s 23 D (2)} 
+ { [K+{S23 + S2i) " K * M w]^o(l) + K\s 23 V (3)} 

+ {[^-(^13 + * 23 ) - ir + M^]P (4) + ir +S23 P (2)} .(3.3.7) 

Note that A m f is exactly cancelled by the contributions of the last two terms in ( |3.2.34| ) 
and (1372739| ). 

Inserting the expressions for the scalar integrals from Appendix |A.3| and using the 
first of the relations ( p. 2. 14 ), various terms, notably all IR divergences and mass-singular 
logarithms, cancel between the real and virtual corrections, and in DPA we are left with 



A mm ' ~ 



sp- 



w 



Ci- 



K~4 



,x w ) + Ci 2 [ — , x w x + Ci 2 



K 4 



K*_ 



x w 



- £i 2 [ -— ,x w + 27riln —2 



K 4 



+ imaginary parts, (3.3.8) 



A m f + Aff 



K + Ks 23 det(Y ) K + det(Y 3 ) ^ K„det(Y 2 ) „ 
— L>o{U) — -^7T7~\ — -^3 1 ; /^^x — -T2 



+ 



det(Y) 
K + K*_s 23 det(Y>)~ 



det(F) 



det(Y') 
+ imaginary parts, 



Vo{0) + det(y) Tz + 



det(Y) 
/Cdet(y 2 / ) 
det(Y') 



(3.3.9) 
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with -Dq(O) and t> (0) given in ( |A.18|) and ( |A.26| ), respectively, and 



■la-A ^ • 

- Li 2 



K_ 

■524 
S23 



2 In 



ie ) + E 

T=±l 



^2 f ^ Zw) - Ci 2 ( r + k ' X W 



•§23 



ie ) In 



M 2 W 



-In 5 



S23 + s 24 



— ie 



+ 2vri 



21n|l-^ g23 +/ 24 Vlnfl+ A ' 



+ ln 



iif* 



VS23 + s 24 



In 1 + 



S23 + S24 



-F9 — F% 



T 2 



K- 



-In 1 



Sl3 + 523 



+ In 



^+S 2 3 



iM^(s 13 + s 23 ) 



(3.3.10) 



The variables /?w, xw an d the dilogarithms Li 2 , Ci 2 are defined in Appendix |A.1| . 

The above results contain logarithms and dilogarithms the arguments of which depend 
on K±. It is interesting to see whether enhanced logarithms of the form hi(K±/Myj) 
appear in the limits K± — > 0. It turns out that such logarithms are absent from the 
non-factorizable corrections, irrespective of the ratio in which the two limits K± — > are 
realized. 

Moreover, the non-factorizable corrections vanish in the high-energy limit. This fea- 
ture of the non-factorizable corrections has been checked by analytical and numerical 
calculations. 



3.3.3 Non-factorizable virtual corrections 

In Chapter |] the factorizable and non-factorizable radiative corrections to four-fermion 
production are investigated. In order to overcome the problem of overlapping resonances 
for the DPA of the real corrections discussed in Section |4.1.3| , the complete real process, 
involving resonant and non-resonant diagrams, is evaluated by Monte Carlo integration. 
On the other hand, only the doubly-resonant part of the virtual corrections is taken into 
account. Therefore, the results of the non-factorizable virtual corrections are required 
separately. 

As in Section |3.3.2j , the non-factorizable virtual corrections factorize to the lowest-order 



cross section: 

d<i rt = ^nf rt d^Bo m (3.3.11) 

with 

= E E (-l) a+fe+1 QaQ b - Re{A virt (A; + , k a , k b )}. (3.3.12) 

a=l,2 6=3,4 71 
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The explicit results of the term A virt = A virt (A;_|_, k 2 ; k-, fa) are listed in the following: 

A virt = A££ + A£? + A£? + A£* + A£ rt + A™2 + A£, rt . (3.3.13) 



with the individual contributions of the diagrams shown in Fig. 3.1 and Fig. 3.2 



A 



virt 



2M| 
k\ — M 2 

-2M^ 



B (k + ,0,M)- B (k + ,0,M) 



k 2 + =M 2 



B' (k +} \ } M W ] 



k\=M^ 



+ {k + <-> 



A virt 



k\=M^ 



a; 



virt 
im 



A 



virt 
mm' 



|C (fc 2 , k+, 0, m 2 , M) - [C (k 2 , k+, A, m 2 , M w 
+ (k 2 <-> fa, k + <-* fc_), 

2A; + p + |c (p + , 0, m e , M) - C (p+, fc + , A, m e , M w ) 
+ (p+ P-,k + <-> - <-> - (p + <-> 

-2A:+A;_ f C (A: +J -k_, 0, M, M) - \C (k+, A, M w , M w 



(3.3.14) 

(3.3.15) 

(3.3.16) 
,(3.3.17) 



a virt 



A£, rt 



2p + k 2 K + D (p + , k + , k 2 , A, m e , M, m 2 ) 

- (p+ «-> P-) - (fo fca, ^+ ^ fe-) + (^a ^ fcs, fc + <-+ k_,p + p_), (3.3.18) 

-2k 2 k_K + D (l) - 2k 3 k + K_D (A), (3.3.19) 

-2k 2 faK + K_E . (3.3.20) 



The other terms, A virt (A; + , can be obtained by obvious substitutions. 

Contributions of the diagrams (if), (im), (mf), and (mm) of Fig. |3.2| cancel in the 
sum of virtual and real corrections in DPA but remain for the virtual corrections. After 
combining several diagrams the result becomes relatively simple: 



A 



virt 



-2 In 



A 



virt 



+ 2 In 



-2 In 
K 



K_ 



4 In 



A 



Ci 2 



K 4 



X W + £«2 



M w 
( K - 



4. 



\K 4 



(3.3.21) 



Li 2 (l - i^) + 7r 2 + ln 2 (-x w ) 



VAM W 



ln(xw) — 27riln(l — x 



+ imaginary parts, 



(3.3.22) 



and 



A m f + Afi rt + A J rt + A$ + A£, rt ~ In 



A 2 



In 



-^23 ~ 1€ 



K + K_ S23 det(y ) n , n s det(y 3 ) J^det(F 2 ) F 



+ 2 In 



det(y) uv y det(y) 

AM w J m U +2 (M^-u) 1 - Ul 



det(y) 



AM- 



w 



/ M+3 (M^-t) \ 



57 



t+2 / V t-3 



-Li 2 1-- ^ -Li 2 1- ^ . (3.3.23) 

where -D (0) can be found in ( |A.18| ) and the variables /?w, an d the dilogarithms Li 2 , 
Ci2 are defined in Appendix [A.l| . 

All mass singularities vanish in the non-factorizable virtual correction factor 5^. In- 
frared singularities remain which have to be cancelled by the corresponding real correc- 
tions. 



3.3.4 Inclusion of the exact off-shell Coulomb singularity 

For non-relativistic W bosons, i.e. for a small on-shell velocity /3\y ; the long range of the 
Coulomb interaction leads to a large radiative correction, known as the Coulomb singular- 
ity. For on-shell W bosons, this correction behaves like l/Ay near threshold, but including 
the instability of the W bosons the long-range interaction is effectively truncated, and the 
l//?w singularity is regularized. Therefore, for realistic predictions in the threshold re- 
gion, the on-shell Coulomb singularity should be replaced by the corresponding off-shell 
correction. The precise form of this off-shell Coulomb singularity reveals that correc- 



tions of some per cent occur even a few W-decay widths above threshold. As explained 
above, the DPA becomes valid only several widths above threshold. Nevertheless, there 
exists an overlapping region in which the inclusion of the Coulomb singularity within the 
double-pole approach is reasonable. 



The Feynman graph relevant for the Coulomb singularity is diagram (d) of Fig. |3TT 
The non-factorizable corrections contain just the difference between the off-shell and the 
on-shell contributions of this diagram. Therefore, the difference between off-shell and 



on-shell Coulomb singularity is in principle included in A mm ', as defined in ( |3.3.4|) . The 



genuine form of A mm / in DPA, which is given by ( 3.3.8|) , does not contain the full effect of 



the Coulomb singularity, because in both Cq functions of (|3.3.4j) the on-shell limit K± — > 



was taken under the assumption of a finite /?w I n order to include the correct difference 
between off-shell and on-shell Coulomb singularity in A mm /, the on-shell limit of the Co 
functions of ( |3.3.4| ) has to be taken for arbitrary /3\y- The full off-shell Coulomb singularity 
can be included by adding 



2M^-s 



2th //3 + Am-A 2vri / K+ + K- + sp w A M 



P \P + A M + PJ Av V 2(3^s 



(3.3.24) 



to the genuine DPA result (|3.3.8| ) for A mm /. The quantities /?, /?, and Am are defined in 
( |A.1| ). After combination with the factorizable doubly-resonant corrections, all doubly- 
resonant corrections and the correct off-shell Coulomb singularity are included. The on- 
shell Coulomb singularity contained in the factorizable corrections is compensated by a cor- 
responding contribution in the non-factorizable ones. Note, however, that this subtracted 
on-shell Coulomb singularity appears as an artefact if the non-factorizable corrections are 
discussed separately. 
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3.3.5 Non-local cancellations 



In Ref. [|36| it was pointed out that the non-factorizable photonic corrections completely 
cancel in DPA if the phase-space integration over both invariant masses of the W bosons 
is performed. This cancellation is due to a symmetry of the non-factorizable corrections. 
The lowest-order cross section in DPA is symmetric with respect to the "reflections" 
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(3.3.25) 



Therefore, A can be symmetrized in 



-K* + or K_ 



-K*_ if the respective invariant 



mass is integrated out. For instance, if k_ is integrated out, A can be replaced by 
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imaginary parts, 



(3.3.26) 



with 2 defined in ( |A.20| ). Note that this expression is considerably simpler than the full 
result for A. In particular, all dilogarithms have dropped out. 

Symmetrizing (|3.3.26|) also in K + <-> — leads to further simplifications if the rela- 
tions ( |3.2.41| ) for the determinants are used. Under the assumptions that (s 13 + s 23 ) > 
(s 23 + s 24 ) > — M^%, and that k w is imaginary, the real part of the result 
vanishes. These assumptions are fulfilled on resonance, k\ = M^; off resonance, there are 
regions in phase space where the assumptions are violated. The volume of those regions of 
phase space is suppressed by factors \k± — M$f\/Myf and thus negligible in DPA. Therefore 
we can use the above assumptions and find that A vanishes in DPA after averaging over 



the four points in the (k+, kl) plane that are related by the reflections ( |3.3.25[ ) 
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(3.3.27) 
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This explicitly confirms the results of Ref . [Bbj . In particular, the non-factorizable correc- 
tions vanish in the limit \k± — <C r^M^r, i.e. for on-shell W bosons. 

The above considerations lead to the following simplified recipe for the non-factorizable 
corrections to single invariant-mass distributions, i.e. as long as at least one of the invariant 
masses k± is integrated out: the full factor A can be replaced according to 



A 



-|A + A 



-K* 



K 4 



-K* 



(3.3.28) 



where the first term on the r.h.s. is given in ( |3.3.26|) , and the second follows from the first 
by interchanging K + «-> and S13 «-> S24. Note that no double counting is introduced, 
since the first (second) term does not contribute if k\ {k 2 _) is integrated out. In order 
to introduce the exact Coulomb singularity, one simply has to add the additional term of 
( EHT24D to ( EHT26D and §]Of). 



3.3.6 Non-factorizable corrections to related processes 

Since all non-factorizable corrections involving the initial e + e~ state cancel, the above 
results for the correction factor also apply to other W-pair production processes, such as 
77 — > WW — > 4 fermions and qq — > WW — > 4 fermions. 

The presented analytical results can also be carried over to Z-pair-mediated four- 
fermion production, e + e~ — ^ ZZ — >■ 4 fermions. In this case, A f f yields the complete non- 
factorizable correction, where all quantities such as M w and T w defined for the W bosons 
are to be replaced by the ones for the Z bosons. The fact that Q\ = Q2 and Q3 = Q4 
has two important consequences. Firstly, it implies the cancellation of mass singularities 
contained in A ff / when all contributions are summed as in (|3.3.2j ). Secondly, it leads to the 
antisymmetry of 5 n f under each of the interchanges k\ <-> ki and k% ^ k^. It is interesting 
to note that (|3.3.2| ) with (|3.3.3| ) are directly applicable, since A m f and A mm / cancel in 
the sum of (|3.3.^ ). Therefore, a practical way to calculate the non-factorizable corrections 
to e + e~ —>■ ZZ — > 4 fermions consists in taking A f f + A m f from (|3.3.9| ) and ( p.3.10| ), and 
setting A mm / to zero. 



3.4 Numerical results 

We used the parameters 
a' 1 = 137.0359895, M z = 91.187 GeV, M w = 80.22GeV, T w = 2.08 GeV, (3.4.1) 



which coincide with those of Ref. [BS , for the numerical evaluation. 

In order to exclude errors, we have written two independent programs for the correction 
factor ( |3.3.2| ) and compared all building blocks numerically. These subroutines are imple- 



mented in the Monte Carlo generator EXCALIBUR p3[0 as a correction factor to the three 

3 Since the Monte Carlo program of Chapter was programmed after the results of the non-factorizable 
corrections were published, the Monte Carlo program EXCALIBUR was used for the results of this chapter. 
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(doubly-resonant) W-pair-production signal diagrams. In all numerical results below, only 
these signal diagrams are included, and no phase-space cuts have been applied. 

If not stated otherwise the results for the figures have been obtained from 50 million 
phase-space points using the histogram routine of EXCALIBUR with 40 bins for each figure. 
For each entry in the tables, 10 million phase-space points were generated. 



3.4.1 Comparison with existing results 

The non-factorizable photonic corrections have already been evaluated by two groups. As 
was noted in the introduction, the authors of Refs. |38|, [39| have not confirmed^ the results 
of Ref. [ 37] . Therefore, we first compare our findings with the results of these two groups. 



Melnikov and Yakovlev give the relative non-factorizable corrections only to the 
completely differential cross section for the process e + e~ — > WW — > v e e + e~v e as a function 
of the invariant mass M + of the u e e + system for all other phase-space parameters fixed. 
All momenta of the final-state fermions lie in a plane, and the momenta k.2 and It3 of the 
final-state positron and electron point into opposite directions. The angle between the 
W~ boson and the positron is fixed to 6\y- e + = 150°, and the CM energy is chosen as 
y/s = 180 GeV. The invariant mass of the e~u e system takes the values M_ = 78 and 
82 GeV. The other parameters are a = 1/137, M w = 80 GeV, and T w = 2.0 GeV. In 



Fig. [3.6| we show our results for the non-factorizable photonic corrections for this phase- 
space configuration. The intermediate-intermediate (mm') corrections agree reasonably 
with those of Ref. [p7| , but the other curves differ qualitatively and quantitatively^. We 



mention that Fig. 3.6 has been reproduced H7| by the authors of Refs. pq, p9j within the 



expected level of accuracy. While the individual contributions shown in Fig. ^6| are at 
the level of 10% owing to mass singularities, the sum, which is free of mass singularities, 
is below 1.2%. 

Beenakker et al. pB| have evaluated the relative non-factorizable corrections to the 
distributions da/dM + dM_, da/dM + , and da/dM av , where M av = (M_ + M+)/2. Our 



corresponding results for the set of parameters given in ( p.4.1|) are shown in Fig. |3T7 



for the single invariant-mass distributions and in Table EO for the double invariant-mass 



distribution. The deviation between the distributions da/dM + and dcx/dM_ in Fig. |3TT| , 
which should be identical, gives an indication on the Monte Carlo error of our calculation. 
The single and double invariant-mass distributions agree very well with those of Ref. |38 



The worst agreement is found for small invariant masses and amounts to 0.03%. In fact, 
the agreement is better than expected, in view of the fact that our results differ from those 
of Ref. |38J by non-doubly resonant corrections. In the numerical evaluations of Ref. 



the phase space and the Born matrix element are taken entirely on shell [|7|. Moreover 



4 While the result of Refs. J3^, ^] (as our result) for the complete non-factorizable correction is free of 
mass-singular logarithms, the result of Ref. p7| contains logarithms of ratios of fermion masses. However, 



the authors of Ref. []37| have informed us [|46| that the results of Ref. [|37| and Refs. [|38| |39| agree for equal 
fermion pairs in the final state. 

5 Although not stated in Ref. [j37| , mass-singular parts have been dropped there in the numerical eval- 
uation jl6| rendering a thorough comparison of the (mf') and (ff') parts impossible. Comparing the sum 
of all three contributions, i.e. the complete non-factorizable correction factor, our result differs from the 
sum read off from Ref. [R7| ■ 
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Figure 3.6: Relative non-factorizable correction factor to the differential cross section 
for the phase- space configuration given in the text. The curves labelled mm', mf', and 
ff correspond to the curves A, B, and C of Ref. P7J, respectively. 
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Figure 3.7: Relative non-factorizable corrections to the single invariant-mass distributions 
d<x/dM ± and da/dM av for the CM energy yfi = 184 GeV. 
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Table 3.1: Relative non-factorizable corrections in per cent to the double invariant-mass 
distribution dcx/dM + dM_ for the CM energy yfs = 184 GeV and various values of M± 
specified in terms of their distance from Mw in units of IV, i.e. A± = (M± — Mw)/Tw- 



the scalar integrals are parameterized by scalar invariants different from ours, leading to 
differences of the order of \k\ — M^|/M^. 

In Ref. fl38f , additionally, the decay-angular distribution d<r/dM_dM + dcos^w+e+ has 
been considered, where #w+e+ i s the decay angle between k + and k 2 in the laboratory 
system. Our results for this angular distribution are shown in Fig. [3.8| . The cross section 
is small for cos#w+e+ ~ — 1, where the corrections are largest. Unfortunately the corre- 
sponding figure in Ref. |38j is not correct ^7|. The authors of Ref. |38| have provided a 



corrected figure, which agrees reasonably well with Fig. |3.8| , but does not show the kinks 
in the curve for M± = 78 GeV. The kinks are due to a logarithmic Landau singularity in 
the 4-point functions. If one employs the on-shell parameterization of phase space, as in 
Ref. []5S| , the Landau singularities appear at the boundary of phase space. Although no 
kinks appear in the physical phase space in this case, the Landau singularities still give 
rise to large corrections for cos#w+e+ ~ — 1- Since the kinks appear in a region where the 
cross section is small, they are not relevant for phenomenology. The issue of the kinks is 
further discussed in Section p. 4. 3 . 



3.4.2 Numerical results for leptonic final state 

In Fig. [3]9] we show the non-factorizable corrections to the single invariant-mass dis- 
tribution do"/dM + for various CM energies. While the corrections reach up to 1.3% 
for y/s = 172 GeV, they decrease with increasing energy and are less than 0.04% for 
y/s = 300 GeV. Note that the shape of the corrections is exactly what is naively expected. 
If a photon is emitted in the final state, the invariant mass of the fermion pair is smaller 
than the invariant mass of the resonant W boson, which is given by the invariant mass 
of the fermion pair plus photon. Since we calculate the corrections as a function of the 
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Figure 3.8: Relative non-f actor izable corrections to the decay-angular distribution 
dcr/dM_dM + dcos6'w+e+ for fixed values of the invariant masses M± and the CM energy 
v/i = 184 GeV. 
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Figure 3.9: Relative non-factorizable corrections to the single invariant-mass distribution 
da/dM + for various CM energies. 
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Figure 3.10: Relative non-factorizable corrections to the angular distribution 
d<r/dM_dM_|_dcos(/> for fixed values of the invariant masses M± and the CM energy 
v/i = 184 GeV. 
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Table 3.2: Shift of the maximum of the single invariant- mass distributions dcr/dM + in- 
duced by the non-factorizable corrections at various CM energies. 



invariant masses of the fermion pairs, the cross section tends to increase for small invariant 
masses and decrease for large invariant masses. 

The non-factorizable corrections distort the invariant-mass distribution and thus lead 
to a shift in the W-boson mass determined from the direct reconstruction of the decay 
products with respect to the actual W-boson mass. This shift can be estimated by the 
displacement of the maximum of the single-invariant-mass distribution caused by the cor- 
rections shown in Fig. |3.9| . To this end, we determine the slope of the corrections for 
M + = Mw, multiply this linearized correction to a simple Breit-Wigner factor, and de- 
termine the shift AM + of the maximum. The smallness of the correction allows us to 
evaluate AM + in linear approximation, leading to the simple formula 



AM, 



d5- 



nf 



dM, 



r 2 

1 w 



M + =M W 



(3.4.2) 



Extracting the slope from our numerical results we obtain the mass shifts shown in 
Table 1X21. 



In Figs. |3.10] and |3.11| we show the effect of the non-factorizable corrections on various 
angular distributions. Since the non-factorizable corrections are independent of the 
production angle of the W bosons, it suffices to consider distributions involving the angles 
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Figure 3.11: Relative non-factorizable corrections to the angular distribution 
da/dM_dM_|_dcos# e + e - for fixed values of the invariant masses M± and a CM energy 
y/s = 184 GeV. 



of the final-state fermions. We define all angles in the laboratory system, which is the 
CM system of the production process. The distribution over the angle <fi between the two 
planes spanned by the momenta of the two fermion pairs in which the W bosons decay, 
i.e. 

(ki x k 2 )(k 3 x k 4 ) 

cos0 = - — —, (3.4.3) 

|ki x k 2 ||k 3 x k 4 | 

is presented in Fig. |3.10| . The corrections are of the order of 1% or less. Like the distri- 
bution, the distribution over the angle between positron and electron ^g+e- (Fig. |3.11| ) is 
symmetric under the interchange of M + and M_. As for the # e +w+ distribution (Fig. |3.8|) , 
the corrections reach several per cent in the region where the cross section is small. 

The distribution in the electron energy E e - is considered in Fig. 3.12 . The corrections 
are typically of the order of 1%. Again the corrections become large where the cross section 
is small. 

In Section |3.3.4j we have introduced a correction term that includes the full off-shell 
Coulomb singularity. The results for the non-factorizable corrections with this improve- 
ment are compared with those of the pure DPA in Fig. 3.13| . For ^/s = 184 GeV the 
additional terms shift the non-factorizable corrections by up to 1.4% for dcr/dM av and by 
up to 0.8% for da/dM + for small invariant masses, whereas for large invariant masses there 
is practically no effect. The difference originates essentially from the differences between 
1/(3 and l//3w m (P-3.24J) . For large invariant masses, the explicit logarithms in (|3.3.24|) 
are small, i.e. the Coulomb singularity correction is minuscule, and this difference practi- 
cally makes no effect. For small invariant masses, the logarithms are approximately m and 
the difference causes the effect seen in Fig. |3.13| . In Table |3^ we show the non-factorizable 
corrections to the double invariant-mass distribution, as in Table |3.1|, but now with the 
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Figure 3.12: Relative non-factorizable corrections to the electron-energy distribution 
da/dM_dM + dE e - for fixed values of the invariant masses M± and a CM energy y/s = 
184 GeV. 
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Figure 3.13: Relative non-factorizable corrections to the single invariant-mass distributions 
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Table 3.3: Same as in Table |3.1| but with improved Coulomb- singularity treatment. 



improved Coulomb-singularity treatment. We find a difference of up to half a per cent 
for small invariant masses but no effect for large ones. We mention that the difference 
between the entries in Tables [O] and |3.1| is directly given by the contribution ( p.3.24| ) to 



A mm i, without any influence of the phase-space integration. 



3.4.3 Discussion of intrinsic ambiguities 

In the results presented so far, all scalar integrals were parameterized by s, S23, S13, S24, and 
k± (parameterization 1). In DPA, however, the parameters of the scalar integrals are only 
fixed up to terms of order k± —M^. We can for example parameterize the scalar integrals in 
terms of s, s 2 3, S123, s 234 , and k\ (parameterization 2) instead. As a third parameterization, 
we fix all scalar invariants except for k± by their on-shell values, corresponding exactly 
to the approach of Ref. 



38 



The results of these three parameterizations differ by non- 
doubly-resonant corrections. 

The difference between parameterizations 1 and 2 is illustrated in Fig. 3.14 for the 
single invariant-mass distribution. The difference amounts to ~ 0.1%. Note that for an 
invariant mass M + = 70 GeV we have a\M^ — k+\/Myj ~ 0.002 and would thus expect 



absolute changes in the non-factorizable corrections at this level. 

For the non-factorizable corrections to the angular distributions, uncertainties of the 
same order are to be expected. The only exceptions are the distributions over the decay 
angles #w+e+ and #w-e-- Let us explain this for 6*w+e+ in more detail: the non-factorizable 
correction contains the term 27ri ln[l+xw^w/( s i3+ s 23)], which can be evaluated by taking 
( s i3 + S23) directly or (s 12 3 — M^) as input. This parameterization ambiguity can lead to 
larger uncertainties, because the above logarithm can become singular, and the location of 
this Landau singularity is shifted by the ambiguity. Since there is a one-to-one correspon- 
dence between S123 and #w+e+ for fixed s and k+, this logarithmic singularity is washed 
out if the angular integration over #w+e+ is performed, but appears as a kink structure in 
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Figure 3.14: Relative non-factorizable corrections to the single invariant-mass distribution 
dcx/dM_|_ for the CM energies 172, 184, and 192 GeV and two different parameterizations. 
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Figure 3.15: Relative non-factorizable corrections to the decay-angular distribution 
dcr/dM_dM + dcos6'w+e+ for fixed values of the invariant masses M± and a CM energy 
a/s = 184 GeV using three different parameterizations, as specified in the text. 
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the angular distribution over this angle. Fig. 3.15 shows the non-factorizable corrections 
to this angular distribution for the three parameterizations. For M + = 78 GeV we still 
have a\k+ — M^\/M^ ~ 0.0004. Apart from the regions where the Landau singularities 
appear, this is indeed of the order of the differences between the three parameterizations. 
When considering the Landau singularities, one should realize that the parameterization 
ambiguity of the locations of the singularities is not suppressed by a factor a, i.e. different 
parameterizations shift the locations at the level of \k+ — M^\/M^ ~ 0.05. However, the 
impact of the corresponding kinks on observables is again suppressed with a\k^_—M^\/M^- 
if the angles are integrated over, since the singularities can only appear near the bound- 
ary of phase space and disappear from phase space exactly on resonance. Since the cross 
section is small where the Landau singularity appears, the effect is phenomenologically 
irrelevant. 



3.4.4 Comparison between leptonic, semi-leptonic, and hadronic 
final state 

The non-factorizable corrections to the invariant-mass distributions are different for dif- 
ferent final states and in general also for the intermediate W + and W~ bosons.0 The 
invariant-mass distributions to the intermediate W bosons coincide only if the complete 
process is CP-symmetric. In this context, CP symmetry does not distinguish between 
the different fermion generations, since we work in double-pole approximation and neglect 
fermion masses; in other words, the argument also applies to final states like i , e e + fi"u^ 
and udsc, which are not CP-symmetric in the strict sense. Thus, we end up with equal 
distributions for the bosons in the purely leptonic and purely hadronic channels, 
respectively, but not in the semi-leptonic case. 

Fig. |3.16| shows the non-factorizable corrections to the single-invariant-mass distri- 
butions for leptonic, hadronic, and semi-leptonic final states at various centre-of-mass 
energies. We observe the same qualitative features for all final states; the corrections are 
positive below resonance and negative above. Quantitatively the differences between the 
corrections to the different final states are small; we note that the slopes of the correc- 
tions on resonance, which are responsible for the shift in the maximum of the distribution, 
are maximal for the leptonic final state. Therefore, we conclude that the W-boson mass 
determination by invariant-mass reconstruction at LEP2 is not significantly influenced by 
non-factorizable corrections. 

The authors of Ref. E8], have also calculated [HTI] the non-factorizable corrections to 



the single-invariant-mass distributions shown in Fig. |3.16| for ^fs = 172 GeV and 184 GeV. 
They find good agreement with our results for positive invariant masses. However, their 
corrections are antisymmetric and therefore differ from our results for negative invariant 



6 In Ref. (3^, |3{| and in the preprint version of Ref. |Q it has been argued that the relative non- 
factorizable corrections to pure invariant-mass distributions are identical for all final states in e + e~ — ► 
WW — > 4 fermions and vanish for Z-pair-mediated four-fermion production. This was deduced from the 
assumption that (up to charge factors) the non-factorizable corrections become symmetric under the 
separate interchanges k\ <-> ki and k^ <-> ki after integration over all decay angles. Although the function 
A for the relative correction has this property, this assumption is not correct, because the differential 
lowest-order cross section is not symmetric under these interchanges. 
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Figure 3.16: Relative non-factorizable corrections to the single-invariant-mass distributions 
da/dM± for e + e~ — > WW — > 4fermions with different final states for various centre-of- 
mass energies. 
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masses. The differences are of the order of non-doubly-resonant corrections and due to 
different parameterizations of the corrections. 

In the previous sections, we investigated the influence of the non-fact orizing corrections 
on various angular and energy distributions with fixed invariant masses for the final- 
state fermion pairs. We have repeated this analysis for hadronic and semi-leptonic final 
states and found corrections of the same order of magnitude, viz. of typically 1% at LEP2 
energies. 



3.4.5 Numerical results for Z-pair production 

For the production channels via a resonant Z-boson pair, e + e~ — > ZZ — > 4fermions, we 
have fi = / 2 and f 3 = / 4 . Owing to Bose symmetry the lowest-order cross section 
dcr Bo rn is invariant under the set of interchanges (ki,k 2 ) «-> (^3,^4). This symmetry, 
which is respected by the non-factorizable corrections, implies that the single-invariant- 
mass distributions to each of the final-state fermion pairs of the two Z-boson decays are 
equal. CP invariance leads to the additional symmetry with respect to (p + ,ki,k 2 ) <-> 
&4, k 3 ); after integration over the Z-pair production angle this substitution reduces 
to (k\,k 2 ) (k4,k 3 ). In view of non-factorizable corrections it is also interesting to 
inspect the behaviour of do"Born under the replacements k\ <-> k 2 and k 3 <-> /c 4 separately, 
since terms in dciBorn that are symmetric in at least one of these substitutions do not 
contribute to d<r n f if all decay angles are integrated over. This is a direct consequence of 
the antisymmetry of <5 n f in each of the substitutions k\ <-> k 2 and k 3 <-> k±, which follows 
from ( |3.3.2j ) and Qi = Q 2 , Q 3 = Q A . 

In order to study the behaviour of do"Born under the replacements k\ k 2 and 
k 3 <-» /C4, it is convenient to consider the helicity amplitudes for the two signal diagrams 
for e + e~ — > ZZ — ► 4fermions, which contain two resonant Z-boson propagators. These 
amplitudes are proportional to the right- and left-handed couplings gf = t> j =F a i of each 
fermion /j = fi, f 3 to the Z boson. As can be seen from the explicit form of the ampli- 
tudes, the substitution k\ <-> k 2 transforms the helicity amplitudes to those with reversed 
helicities of the fermions f\ and f 2 = f\ apart from changing the couplings gf into gf. 
Therefore, the differential lowest-order cross section, i.e. the squared helicity amplitudes 
summed over all final-state polarizations, can be split into two parts: one is symmetric 
in k\ «-> k 2 and proportional to [(gf) 2 + (g\ ) 2 ]/2 = v\ + a 2 , the other is anti-symmetric 
and proportional to [(gi) 2 — (gf) 2 ]/2 = 2v\d\. The analogous reasoning applies to the 
substitution k 3 <-* k^. After performing the angular integrations, we finally find that the 
lowest-order cross section is proportional to (v\ + af)(v 2 + a^), and the non-factorizable 
correction proportional to 4QiViaiQ 3 v 3 a 3 , where the charge factors Qi stem from the cor- 
rection factor 5 n f. Comparing pure invariant-mass distributions for different final states, 
the ratios of the non-factorizable corrections should be of the same order of magnitude as 
the ratios of the corresponding coupling factors, 



4QiViai<33t>3a 3 



(vl + a\)(vl + al) 



(3.4.4) 
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Figure 3.17: Relative non-factorizable corrections to the single- invariant-mass distributions 
der/dMi^ for e + e~ — > ZZ — > 4fermions with different final states for y/s = 192 GeV. 



The factors F take the following values: 
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(3.4.5) 



where £, u, d generically refer to leptons, up-type quarks and down-type quarks, respec- 
tively. The reason for the smallness of the factors F is different for leptons and quarks: 
for leptons the suppression is due to the small coupling Vi to the vector current, for quarks 
the factor F is reduced by the relative charges Qi. 

Fig. p. 17] shows the non-factorizable corrections to the single-invariant-mass distribu- 
tions du/dMi^, where Mi^ denote the invariant masses of the first and second fermion- 
anti-fermion pairs, respectively. The ratios of the different curves are indeed of the 
order of magnitude of the ratios of the factors F given in (|3.4.5|) . For equal signs of 
Qi and Q3 the shape of the corrections is similar to the shape of the corrections to 
e + e~ — > WW — > 4fermions, for opposite signs of Qi and Q3 the shape is reversed. The 
corrections by themselves are very small and phenomenologically unimportant. The small- 
ness of these corrections can be qualitatively understood by comparing the factors F of 
( 3.4.4 ) for the ratios of the couplings with the corresponding one for the W-pair-mediated 
processes. For e + e~ — > WW — > 4 leptons we simply have F = 1, because in the LEP2 
energy range the purely left-handed t-channel diagram dominates the cross section, and 
no systematic compensations are induced by symmetries. Therefore, the factors in ( |3.4.5| ) 
should directly give an estimate for the suppression of S n f for e + e~ — > ZZ — > 4fermions 
with respect to four-lepton production via a W-boson pair. Comparing the corrections 
for energies with the same distance from the respective on-shell pair-production thresh- 
olds, i.e. the curve for y/s = 184 GeV in the W-boson case (Fig. |3.16| ) with the curves 
for a/s = 192 GeV in the Z-boson case (Fig. |3.17|) , we find reasonable agreement with 
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Figure 3.18: Relative non-factorizable corrections to the angular distributions 
da/dMidM 2 dcos0 and dcr/dMidM 2 dcos6' M + r ~ in e + e~ — > TIL — > h~h + t~t + for fixed 
values of the invariant masses M 12 and yfs = 192 GeV. 



our expectation. The authors of Ref. |3^, |39| have reproduced the corrections shown in 
Fig. |3.17| with good agreement [57[ . 

Finally, we inspect the impact of non-factorizable corrections to some angular distribu- 
tions in Z-pair-mediated four-fermion production for fixed values of the invariant masses 
Mi 2- Since the presence of the suppression factor F relies on the assumption that the 
phase-space integration is symmetric under k\ <-> k 2 and k 3 <-> fc 4 , this suppression in 
general does not apply to angular distributions. However, partial suppressions occur, e.g., 
if the integration is still symmetric under one of these substitutions and, in particular, 
for quarks in the final state because of their smaller charges. Two examples for angular 
distributions without any suppression are illustrated in Fig. [3.18] for the purely leptonic 
final state yT fi + r~r + . The angle <fi is defined by the two planes spanned by the momenta 
of the two fermion pairs in which the Z bosons decay, 

(ki x k 2 )(k 3 x k 4 ) 

cos0 = — — -, (3.4.6) 

|ki x k 2 ||k 3 x k 4 | 

and 6fj,+ T ~ denotes the angle between the momenta of the /i + and the r~, respectively. The 
shapes of the curves in Fig. 3.18j , specifically the curves for the distribution in cos <ft, nicely 



reflect the approximate anti-symmetric behaviour in the angular dependence, which leads 
to the suppression in the invariant-mass distributions. The size of the corrections turns 
out to be at the level of a few per cent, i.e. they are not necessarily negligible in precision 
predictions. Note, however, that the cross section for Z-pair production is only one tenth 
of the W-pair production cross section. 
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Chapter 4 



Radiative corrections to 
W+W" 4/ 



As discussed in the introduction, the cross section for e + e~ — > W + W~ — > 4/ should be 
known with an accuracy of 1% or better in order to cope with the experimental precision 
at LEP2. This requires the calculation of the 0(a) corrections to the processes e + e~ — > 
W+W- -> 4/. 

The leading radiative corrections, like the running of the electromagnetic coupling 
constant, universal corrections associated with the p parameter, the Coulomb singularity, 
and the resummed leading-logarithmic corrections from initial-state radiation (ISR), are 
already implemented in several event generators (see Ref . [ 30| and references therein) . The 



neglected non-leading corrections can be estimated from the on-shell W-pair production 
, where the full one-loop calculation differs from the one including only these leading 
effects by about 1-2% at LEP2 energies and 10-20% in the TeV range. Hence, the leading 
radiative corrections are in general not sufficient to match the experimental accuracy of 
LEP2. 

The calculation of the full 0(a) corrections is extremely complicated since the number 
of diagrams is of the order of 10 3 -10 4 including one-loop integrals with up to six propa- 



gators [Q. The numerical computation of these corrections is highly non-trivial because 
the tensor reduction to scalar integrals and the calculation of the one-loop integrals with 
five and more propagators cause numerical instabilities. Furthermore, the numerical inte- 
gration via Monte Carlo techniques is rather slow for such complex calculations. 

Approximations for radiative corrections beyond the leading level have been calculated 
by two groups. A first calculation of the complete doubly-resonant O{o) radiative correc- 
tions to the four-fermion processes e + e~ — > W + W~ — > 4/ has been discussed in Ref. [50|. 



There a semi-analytic approach has been used with different matrix elements for different 
phase-space regions and for different observables. Moreover, the four-fermion phase space 
has been factorized into the phase space of the on-shell W-pair production and of the 
on-shell W decays and the integrations of the two invariant masses of the intermediate 
W bosons. For the numerical discussion, only results for leptonic processes have been 
included. The authors of Ref. |50| have found a large shift of the Breit-Wigner line shape 
due to logarithms of the form \n(mj/s) resulting from final-state radiation (FSR). These 
logarithms are due to the absence of collinear photons in the definition of the invariant 
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mass of the W boson. In realistic observables, collinear photons cannot be resolved from 
charged fermions, except for muons, and have to be included in the reconstructed invariant 
mass of the corresponding W boson. 

In Ref. ||51|| , a four-fermion event generator has been presented, including all 0(a) 
radiative corrections to the on-shell W-pair production, e + e~ — > W + W~, with exponenti- 
ation of the universal corrections from photon radiation off initial-state e ± and off the W ± 
bosons. Leading-logarithmic corrections from FSR have been included in this Monte Carlo 
generator in Ref. [|52|. However, non-leading 0(a) corrections to the W decays and the 
non-factorizable corrections are missing. The results of Ref. [|5tJ have been qualitatively 
confirmed by the calculations of Ref. |52| . 



Hence, a Monte Carlo generator including the complete doubly-resonant 0{a) correc- 
tions is needed in order to match the accuracy of LEP2 and to take into account realistic 
experimental situations. 



4.1 Strategy of the calculation 

In this chapter, we consider the virtual corrections to the processes 

e + (p + )+e-(p_) - W+(fc + )+W-(M - fi(h) + f 2 (k 2 ) + f 3 (k 3 ) + U(h) (4.1.1) 

and the complete bremsstrahlung processes 

e + (p + )+e-(p_) - fi(k 1 )+f2{k 2 )+Mk 3 ) + U{kt)+>y{q), (4.1.2) 

where the relative charges of the fermions are represented by Qi with i = 1, ... ,4. 
The masses of the external fermions are neglected, except where this would lead to mass 
singularities. For the virtual corrections, the momenta of the intermediate W bosons read 

k + = ki + k 2l k_ = k 3 + k A . (4.1.3) 

Furthermore, the square of the complex W-boson mass is defined by M 2 = — iM-^T^, 
and the center-of-mass (CM) energy is ^/s. 



4.1.1 Doubly- resonant virtual corrections 

The diagrams of four-fermion production can be classified into the doubly-resonant, singly- 
resonant, and non- resonant diagrams according to the number of resonant W-boson prop- 
agators. Hence, the amplitude of the virtual corrections can be written in the following 
way after implementation of the finite W-boson width: 

M- ^±Mi^S | g±Migij) | R^K^e) 

MvlTt ~ (4_M2)(jfcl -M2) + kl-M* + kl-M* + »W+>*-> & b v • • / 

v ' ' non-resonant 

doubly-resonant single-resonant 

where the variable 9 symbolizes all phase-space variables, except for the invariant masses 
fcl of the W bosons. 
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Since the matrix elements of non-doubly-resonant diagrams are suppressed by a factor 
oT w ln(. . .)/(7rM w ) ~ 0.1%, a reasonable approach is to include only the doubly-resonant 
diagrams into the calculation. However, the naive inclusion of the doubly-resonant dia- 
grams yields gauge-dependent results and the reliability of this approximation is unclear. 

In order to separate the doubly-resonant corrections in a gauge-invariant way, the pole 
scheme has been proposed in Ref. [[11], |12[1 , where the matrix element is expanded about 
the complex W-boson masses. Since the complete matrix element is gauge- invariant, the 
single terms of the pole expansion are also gauge-invariant. If only the leading term is 
kept, the expansion is known as the double-pole approximation (DPA): 

<DPA i? + _(M 2 ,M 2 

(k\ — M 2 )(k 2 _ — M 2 ) 

Note that the residue is taken at the complex pole resulting in complex kinematical 
invariants. In order to avoid the calculation of one-loop integrals for complex invariants, 
the W width is neglected in the numerator of Q4.1.5|) , where the neglected terms are 
suppressed by a factor r^/Mw and thus negligible in DPA. 

For the application of the DPA on the virtual corrections, a set of eight independent 
phase-space variables including k\ has to be chosen, which determine the momenta of 
the final-state fermions uniquely. For several choices of the parameterization of the phase 
space, the DPA differs only in non-doubly-resonant contributions. Note that in general 
events can be outside of the physical phase-space boundaries for on-shell W bosons (see 
Section 3.4.3|) . A relatively simple choice that generates only events within the physical 



^virt A = TT^kni ' U - (4-1.5) 



phase-space domain is given in Appendix JTJ. 

The DPA is not valid near threshold since the phase-space region where both W bosons 
are resonant is suppressed by the kinematical factor Az (s, k\ , fc 2 .) (see ( |C.8| )). Therefore, 
the singly-resonant corrections are as important as the doubly-resonant corrections at 
threshold. On the other hand, for some processes the non-doubly-resonant contributions 
can be enhanced by nearly on-shell virtual photons. This enhancement can be avoided by 
introducing appropriate phase-space cuts in the calculation. 

In DPA the virtual corrections can be classified into factorizable and non-factorizable 



ones 



12[, |42]| - The square of the matrix element of the virtual corrections reads in DPA 
l-Mvirt A | 2 = \M t | 2 + \M Bom \ 2 6 v j\ (4.1.6) 

where A-^f, A^Bom, and 5^ rt are the matrix elements of the factorizable corrections, the 
Born matrix element defined in ( |3.1.8| ), and the non-factorizable correction factor, respec- 
tively. 

Factorizable virtual corrections 

The factorizable virtual corrections are defined by the product of the on-shell matrix 
elements of the W-pair production and the W decays and the (transverse parts of the) 



W propagators (see Fig. 



+ M£:; n - W+W -5M W+ ^MZrn M ' 4 + MtTn^ W+W ~MZn fj2 SM W -^ 
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On-shell production 



On-shell decays 



Figure 4.1: Factorizable corrections to e + e — > 4/ 



^^^ 2 , and 5A4 W ^ 3 -^ 4 denote the one-loop amplitudes of 
The sum runs over the physical polarizations A± of the 



where 5M e+e ~^ w+w ~ , 5M W+ 
the respective subprocesses. 
W 1 * 1 bosons. 

Since the on-shell W-pair production and the on-shell W-decays are gauge-invariant 
processes, the factorizable corrections are also gauge-invariant. The results of the on-shell 
W-pair production and the on-shell W-boson decay are explicitly given in Ref. |53| and 
Ref. [p4"|, respectively, and can be implemented in the Monte Carlo program, where we 



have to take care of the spin correlations of the W bosons in the production and decay 
subprocesses. The virtual corrections are build up by scalar form factors, which include 
all one-loop integrals, and standard matrix elements, which depend on the polarization 
vectors and spinors of the external particles. In DPA the form factors depend exclusively 
on the W-production angle and can be evaluated very fast by an expansion in Legendre 



polynomials [21 



Non-factorizable virtual corrections 

The non-factorizable virtual corrections are defined as the difference between the virtual 
corrections to the complete four-fermion process in DPA and the factorizable virtual cor- 
rections. A representative set of diagrams contributing to the non-factorizable corrections 
are shown in Figs. |3.1| and |3.2| . The matrix element of the non-factorizable virtual correc- 
tions .Mnf rt factorize to the Born matrix element in DPA (see Chapter |]for more details): 

Mt? = M BoT J™\ (4.1.8) 



where the correction factor ^ rt is relatively simple and explicitly given in Section |3.3.3| . 



Since the factorizable and the complete virtual corrections are gauge-invariant, the non- 
factorizable corrections are also gauge- invariant. 

Note that the non-factorizable corrections involve logarithms of the form \a(k\ — M 2 ), 
which become singular in the limit k± — > M 2 . Therefore, these corrections are calculated 
for off-shell phase space, i.e. k± ^ M 2 , while the off-shellness (k± — M 2 ) is neglected 
whenever possible. 
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4.1.2 On the definition of the reconstructed W-boson mass 



In Chapter the non-factorizable virtual and real corrections have been calculated in 
DPA, where the integration over the photon momentum of the real corrections is performed 
analytically. As for the virtual corrections, the non-factorizable real corrections depend on 
the choice of the independent phase-space parameters in DPA. These parameters are fixed 
while the integration over the photon momentum is performed. To apply the DPA on the 
non-factorizing real corrections, two of these parameters are identified with the invariant 
masses of the W bosons. Owing to the presence of the bremsstrahlung photon the definition 
of these invariant masses is not unique. For different definitions the result for the real non- 
factorizable corrections differs by doubly-resonant contributions, as discussed in detail in 
Section |3.1.4j . The reason is that the resonant W-boson propagator l/[(fc+ + q) 2 — M 2 ] 
is constant, if the invariant mass is defined by (k + + q) 2 , but depends on the photon 
momentum q for the invariant mass k+. Therefore, the calculable observables are restricted 
to the actual parameterization of the phase space in a semi-analytic calculation. For 
instance, only invariant-mass distributions, where the W-boson masses are defined by k 2 ±1 
can be calculated with the results of Chapter [|. 

A DPA including all 0(a) corrections for four-fermion production has been worked out 



in Ref. ||50|| . As in Chapter ^ a semi-analytic approach has been used in Ref. ||50|| , where, 
for the invariant mass distributions, the W-boson masses have been defined by invariant 
masses of the fermion-antifermion pairs, i.e. M\ = k±, resulting in large corrections from 
collinear radiation of bremsstrahlung photons off final-state fermions. These originate 
from logarithms of the form ln(m^/s) and yielding large distortions of the peak position 
of the Breit-Wigner line shape at the CM energy 184 GeV of -77MeV, -38MeV, and 
— 20MeV for e + z/ e , p^v^ and t + u t , respectively. In practice, the invariant masses of W 
bosons have to be reconstructed in a more realistic way. 

In realistic experimental situations, the final-state quarks are observed as jets. Photons 
radiated collinear to these jets cannot be resolved, and the photon momentum is included 
in the jet momentum. On the other hand, the momentum of the neutrino can only be 
calculated from the missing momentum of the final-state particles. In this way the neutrino 
momentum includes also the momenta due to emission of photons collinear to the beam 



55] . Note that the reconstruction of the neutrino momentum is possible for semi-leptonic 
final states, but not for purely leptonic final states, where two neutrinos are involved. 
Furthermore, it is very difficult to separate collinear photons from electrons. However, 
photons can be resolved from muons even when they are collinear. Therefore, only a few 
observables, like the invariant-mass distribution d<r/dM_ with = h 2 _ of the process 
e + e~ — > u d ii~v i x- l are directly sensitive to the distortion due to large corrections from FSR. 
In any case, a Monte Carlo generator is required in order to take into account realistic 
experimental situations. 



4.1.3 Overlapping resonances in the bremsstrahlung process 

The definition of a suitable approximation for the real corrections is made difficult by over- 
lapping resonances of the W-bosons. In the bremsstrahlung process, the W-boson propaga- 
tors differ depending on whether the bremsstrahlung photon is emitted from initial-state 
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Figure 4.2: Diagram with overlapping resonances 

particles or from the decay products of the W boson. Both types of propagators are 
present in Fig. |4.2j where the poles of the propagators are located at k 2 _ = M 2 , k\ = M 2 , 
and (k + + q) 2 = M 2 . Therefore, the expansion of the cross section about k 2 ^ = M 2 and 
k 2 _ = M 2 , like in the case of the virtual corrections, is not suitable. 

The matrix element of the bremsstrahlung process can be written in the following way 
(see Section |3.1.5|) : 



Tl{k + ,k_,q) H+{k + ,k_,q) 

•I** Teal / i O a *-n \ / i Q * *-r>\ l" 



(fc+ — M 2 )(k 2 _ — M 2 ) [(k + + q) 2 — M 2 ](k 2 — M 2 ) 

where M includes the matrix elements of all non-doubly-resonant diagrams. The factor- 
izable corrections of the cross section are the squares of the single terms corresponding to 
TZ, 7Z + , and 1Z-, while the non-factorizable corrections are the interferences between these 
terms. 

In the hard photon region, E 1 3> Tw, the two resonances of the W + bosons located 
at k 2 + = M 2 and (k + + q) 2 = M 2 are well separated in phase space. Only the factorizable 
corrections contribute in this region in DPA, and the non-factorizable corrections vanish 
in DPA. For soft photons, E 1 -C Tw, the resonances almost coincide, and the photon 
momentum can be neglected in the resonant propagators. However, for semi-soft photons 
with an energy E 7 = (9(T W ), the resonances overlap. Thus the definition of an appropriate 
approximation is rather complicated and, hence, the reliability becomes unclear. 



4.1.4 Inclusion of the real corrections 

In order to avoid the problem of overlapping resonances, the bremsstrahlung process is 
taken into account exactly. 

The bremsstrahlung processes are already studied in Chapter ||] including all diagrams. 
For the numerical discussion, the processes e + e~ — > v^ii + t~v T) e + e~ — ► ud/i'z/^, and 
e + e~ — > udsc are considered. These are the processes with the smallest number of 
diagrams within the leptonic, semi-leptonic, and hadronic process classes. Note that the 
radiative corrections for any leptonic process are equivalent to the radiative corrections for 
e + e~ — > v^^t'Vj. in DPA, in the absence of logarithms ln(m 2 /s) from FSR. The same 
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is valid for e + e -> ud/j and e + e -> udsc concerning semi-leptonic and hadronic 
processes, respectively. 

4.2 Subtraction method 

When calculating radiative corrections for four-fermion production, one has to take care 
of the singularities of the real and virtual radiative corrections. According to the Bloch- 
Nordsieck theorem R| the infrared singularities cancel between the virtual and real part. 



Furthermore, collinear singularities, i.e. mass singularities, appear in the radiative cor- 
rections and are regularized by the fermion masses which are neglected otherwise. These 
singularities originate from collinear photon emission from light external fermions and show 
up as large logarithms \n{mj/s). For inclusive enough observables, these collinear singu- 
larities cancel between the real and virtual radiative corrections owing to the Kinoshita- 



Lee-Nauenberg theorem fl57|, except for the mass singularities remaining from initial-state 
radiation (ISR) or resulting from renormalization. Two methods are exist for the treatment 
of the singularities in a Monte Carlo generator: phase-space slicing and the subtraction 
method. 

Phase-space slicing splits the integration domain by a small separation cut into a part 
which includes all singularities and a finite part. For infrared singularities the splitting 
is usually done by a cut on the photon energy, for collinear singularities by a separation 
angle between the photon and the particle which emits the photon. The integral of the 
singular part over the photon momentum is analytically performed with an appropriate 
regularization, and terms of order of the cut parameter are neglected. In this way the 
singularities are extracted and can be added to the virtual part of the radiative corrections. 
Both the singular and the finite part of the radiative cross section depend logarithmically 
on the cut parameter. This dependence must cancel in the complete result. The logarithms 
of the separation cut are determined in the singular part by analytical integration. On 
the other hand, they must be compensated by numerical calculation of the corresponding 
logarithms included in the finite part of the cross section. Therefore, the cut has to be 
small enough that terms of order of the cut parameter can be neglected in the singular 
part, but large enough to obtain numerically stable results. 

In contrast to phase-space slicing, the subtraction method requires no cut parameter. 
We use the subtraction method since we expect a better convergence behaviour. This 
method has been mostly applied for NLO predictions in massless QCD [J58, [59|][] where 



the singularities are usually regularized dimensionally. In Electroweak Standard Model 
processes it is more convenient to introduce an infinitesimal photon mass as a regularization 
parameter for infrared singularities and to use small fermion masses for the regularization 
of the collinear singularities. In order to apply the subtraction method to four-fermion 
production, this method is formulated for mass regularization. 

In the following, only singularities for photon emission off nearly massless external 
fermions are considered. Other singularities like collinear singularities of diagrams where 
a virtual photon decays into two external fermions have to be excluded by appropriate 



lr The subtraction method for massive particles can be found in Refs. 
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cuts. The subtraction method described in the following is based on the dipole formalism 
of S. Catani and M.H. Seymour |58|. The results of this section have been worked out 
independently by S. Dittmaier ||60|| . Comparing both results we find full consistence. 

In the dipole formalism, various mappings are constructed from the five-particle phase 
space of the bremsstrahlung process e + e~ — > 4/7 into the four-particle phase space of 
the non-radiative process e + e~ — > 4/. These mappings are required to obtain a process- 
independent formulation of the subtraction method. Each mapping corresponds to a 
certain singular behaviour of the differential cross section. According to these mappings, 
the five-particle phase space $^ 5 - ) splits into a four-particle phase space $^ 4 - ) and a remaining 
one-particle phase space, which contains the singularities: 



d $(5) _ J dxidx 2 J d$ (4) (xipi,x 2 p 2 ) j d$ ifc 



(4.2.1) 



with i, k = 1, ... 6. In the following, particles 1, 2 are initial-state fermions and 3, . . . , 6 
final-state fermions. For the discussion of the subtraction method, all particle momenta 
are denoted by pi and the photon momentum by q. 

The arguments of d$^ 4 ^(xij>i, X2P2) indicate that the final-state momenta are calculated 
for incoming momenta Xipi and X2P2- The one-particle phase space is decomposed into 
integrations over and over the momentum fractions of the incoming momenta X\ and 
X2- In this way, all different splittings of the five-particle phase space fit in this formula. 
The variables x\ and X2 are partly fixed by ^-distributions included in definition of 

For each mapping a subtraction term Vik is constructed in such a way that it matches 
the singular behaviour of the cross section of the bremsstrahlung process in a certain 
phase-space region and that it can be analytically integrated over The subtraction 
term is subtracted from the real corrections and added to the virtual corrections after 
analytic integration over 



d$ (5) 
d$ (5) 



- /w Born 



O (^) + I d*«2 Re {M^Mgt} O ($ 



(4) 
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(5) 
Born 



i,k=l 
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+ ^ dx!dx 2 j d$ (4) (x lPl , x 2 p 2 ) O ($ (4) ) 

x {2 Re (M^ t M { £: n ) 5(1 - Xl )5(l - x 2 ) + E f / d$ lfe V., 

A 1 1 -J 



i,k=l 



(4.2.2) 



where M^l^, Mg] m , and M.^ are the matrix elements of the tree-level process, the 
bremsstrahlung process, and the virtual corrections, respectively. The experimental situa- 
tion, e.g. cuts, are included in the definition of the observable O. Note that the observable 
O depends on the momenta of the four-particle or five-particle phase spaces. All integra- 
tions are performed numerically, except for the integration over which is performed 
analytically. 

In this way, a part of the real corrections including the singularities is transferred to the 
virtual corrections. The difference between the cross section of the bremsstrahlung process 
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and the subtraction terms includes no singularities and hence can be integrated over the 
whole phase space. Note that it is not allowed to exclude the soft-photon region from 
the integration domain where the cross section of the bremsstrahlung process becomes 
infrared singular, since the infrared singularities have to cancel between the real and 
virtual corrections. This is in accordance with the fact, that it is not possible to separate 
experimentally photons with infinitesimal small energies from charged particles. 

The situation for collinear singularities is different since collinear singularities show 
up as large logarithms of small but finite fermion masses. If fermion masses are used as 
regularization parameters, the phase-space region where the photon becomes collinear to 
an external fermion has to be included in the integration domain. 



4.2.1 Behaviour of the cross section for collinear photons 

As a first step, the collinear limit of the cross section for ISR is considered. It is convenient 
to calculate the collinear limit of the cross section in an axial gauge, where the photon 
polarization sum runs only over transverse polarizations and hence interference contribu- 
tions do not involve collinear singularities. A photon with momentum q is emitted from 
an unpolarized initial-state fermion with momentum p a , mass m a , and relative charge Q a . 
The relevant part of the cross section multiplied by the polarization sum of the photon 
reads 



ie(5a)7^(#a + m a )ieQ a 'y u 
1 - x)4 a - 4 ( x 




1 



x 




XI, 



PaQ, 



with n 2 = 0, qn ^ 0, x = 1 



(nq)/(np a ), and the splitting function is given by 
1 + x 2 



P(x) 



X 



(4.2.3) 



(4.2.4) 



The arbitrary vector n is introduced to define the polarization of the photon. 

The factor (1 — x)^ a — $ vanishes in the collinear limit. Only the last term on the 
right-hand side proportional to x^ a is singular in the collinear limit q — > (1 — x)p a . Hence, 
the collinear limit for photon radiation off unpolarized initial-state particles yields 



M {5) 
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(4.2.5) 



The behaviour of the cross section for polarized particles in the collinear limit can be found 

e 2 Ql 



in Ref. 62 
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(4.2.6) 
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Kf = + 




Figure 4.3: Value of Kf for the different photon radiation off fermions 



with the helicity a a of the initial-state fermion a. In the following all results are calculated 
for helicity eigenstates of the external fermions. 

For massive fermions the collinear limit involves a helicity-flip term proportional to the 
square of the fermion mass, called finite-mass term. Note that if the cross section of the 
bremsstrahlung process is calculated for vanishing fermion masses the finite-mass terms 
are missing and have to be added to the cross section of the bremsstrahlung process or 
skipped in the corresponding subtraction terms in ( |4.2.2 ). 

In order to be able to combine the collinear and soft limit in Section [4.2.3| , we use 
charge conservation 



(4.2.7) 
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to obtain 
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M^Lixpa,-^) 2 .(4.2.8) 



The sign K t refers to charge flow of fermion i into or out of the diagram, respectively, which 
is illustrated in Fig. |4~3| . 

For FSR, q has to be replaced by — q: 



■M collinear (Pi) a ii Q) 
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P(z) 



1 + 2r m\ 
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-^Born — ,<?% 
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+ 



(4.2.9) 



with z = (npi)/(npi + nq), where the particle i is a final-state particle. 
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Figure 4.4: Diagrams for different combinations of initial-state or final-state emitter and 
spectator 



4.2.2 Soft-photon approximation of the cross section 

The infrared limit is given by the well known soft-photon approximation 
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After applying partial fractioning 

PiPk 



PiPk 
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(4.2.10) 



(4.2.11) 



(PiQ) (Pkq) [(pi + Pk)q] (pkQ) [fa + Pk)q] faq) 

and distributing the finite-mass term with i = k with the help of charge conservation 



( 4.2.7 ) to the terms with i ^ k, the soft limit reads 
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Born 



(4.2.12) 



In this way, the soft-photon approximation ( |4.2.12j ) is written in a similar way as the 
collinear limits flPTgp and (jXg). 



4.2.3 Construction of the subtraction terms 

The subtraction terms have to be constructed in such a way that they reproduce both 
the soft and collinear limit of the cross section. The soft-photon limit ( [4.2. 12| ) depends 
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on two momenta Pi and p k . Particle i is called emitter and k spectator, since the soft- 
photon approximation is singular if the photon becomes collinear to the emitter, but not 
if the photon is collinear to the spectator. Since the mass of the spectator regularizes 
no singularities, it is neglected in the following calculations. The diagrams corresponding 
to the four possible combinations of initial-state or final-state emitter and spectator are 



shown in Fig. 4.4 



Although the following results are formulated for four-fermion production, they can be 
applied to all processes where the external particles are light charged fermions and any 
other massless or massive neutral particles. Note that it is some freedom in the definition of 
the subtraction term and the separation of the one-particle phase space. In the following, 
all initial-state particles are denoted by the indices a or b with a ^ 6 and a,b = 1,2. 

Final-state emitter with final-state spectator 

Firstly, the case is considered where the emitter and spectator are final-state particles. A 
suitable subtraction term is defined byQ 
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with the momenta of the four-particle phase space 

Vik 



Pi = Pi + q 
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(4.2.13) 
(4.2.14) 



and the variables 



Z-ik 



PiPk 



Vik 



PiQ. 



PiPk + Pkq PiPk + Piq + Pkq 

where the indices i and k mark the emitter and the spectator, respectively. 



(4.2.15) 



The subtraction term obeys the soft limit (|4.2.12|) for z ik — > 1 and — » 0, which can 
be verified with the help of 



Piq + Pkq + PiPk 



(4.2.16) 



1 - z ik {l - y ik ) (Pi + Pk)q 

Furthermore, if the momentum n after (|4.2.9| ) is identified with the momentum of the 
spectator p k , the subtraction term reproduces the cross section in the collinear limit, i.e. 
y ik -> 0, z ik -> z, and pi -> Pi/z, ( |4.2.9|) . 

The subtraction term is proportional to the Born matrix element A^Bom which depends 
on the momenta pi and p k of the four-particle phase space. These momenta are fixed for 
the analytical integration over The momenta pi and p k are chosen in such a way that 
they fulfil momentum conservation 

= K + pi+p k + q = K + pi+p k (4.2.17) 

2 Note that the actual definition of the subtraction term of Ref. |60| differs by a factor 1 — yi k . 



N6 



and the on-shell conditions of the four-particle phase space 

m 2 = m 2 + A 2 ^0, m 2 k = 0. (4.2.18) 

The momentum K is the sum of the momenta of the external particles, except for the 
emitter, the spectator, and the photon. The infinitesimal photon mass is denoted by A. 

The next step is to relate the momenta of the five-particle phase space to the momenta 
of the four-particle phase space and the variables z ik and y^ k of the one-particle phase 
space. Therefore, we apply the mappings 
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(4.2.20) 
(4.2.21) 



1 - ZikPiPk 

with k±pk = k±q = and find the splitting of the five-particle phase space into the four- 
particle phase space determined by the momenta pi and pk and the remaining one-particle 
phase space: 

d 4 pj d 4 p/c d 4 g 
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with the integral over the separated one-particle phase space 
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(4.2.23) 



2ifc(l - Zik)2pip k 

The additional integrations over X\ and x 2 are introduced in order to match the definitions 
of ( |4.2.lD . The masses A and rrii in the ^-distribution of (|4.2.22| ) can be neglected since they 
do not regularize singularities. The step function 9 in ( |4.2.23j ) results from the requirement 
k\ < of ( [4.2.21|) and regularizes all singularities. It can be verified in the rest frame 
of Pi + p k that k± is a space-like vector. The variable <pk ± is the solid angle of k± with 
respect to pk in the rest frame of pi + p k . Since the subtraction term does not depend on 
<pk ± the integration over <pk ± yields 2tt. Note that the Born cross section included in the 
subtraction term Vik depends exclusively on the momenta of the four-particle phase space, 
which are fixed for the integration over z\ k and yi k . 

Finally, the subtraction term is integrated over the one-particle phase space yielding 
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(4.2.24) 
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with 
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2pip k ) 2 \mf ) 2 \2pip k j 
where the photon mass is neglected with respect to the fermion masses, i.e. A <C mi. 

Final-state emitter with initial-state spectator 

Next, the subtraction term for final-state emitter and initial-state spectator is considered. 
The subtraction term with the correct collinear and soft limit is defined by 
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with the momenta of the four-particle phase space 

Pi = Pi + q ~ (1 - X ia )p a , p a = X ia p a , 

and the variables 

PiPa PiPa + p a q - Piq 
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Similarly to the previous section, the momenta pi and p a fulfil momentum conservation 
and the on-shell requirement for the external particles of the four-particle phase space. 

In order to separate the integral over the one-particle phase space from the integral 
over the five-particle phase space, the following mappings are used: 
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with k±p a = k±q = 0. They result in the splitting of the phase space 
d 4 pi d A q 
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(4.2.32) 



where 4>k ± is the solid angle with respect to p a in the rest frame of p a + p^ Since the 
spectator is an initial-state particle, the four-particle phase space becomes x a - dependent 
and the one-particle phase space is included in the ^-integration. As in the previous 
case the ^-function originates from k\ < of the mapping ( 4.2.30Q and regularizes all 
singularities. Furthermore, the masses rrii and A can be omitted in the 5-distribution of 

Since the subtraction term does not depend of <pk±i the integration over <fti :± yields 
2/T. Moreover, the Born cross section included in the subtraction term does not depend 
on Zi a , but is a function of Xi a . In order to avoid the analytic integration over the Born 
cross section, the subtraction term is integrated for fixed Born cross section at x ia = 1. 
The Xi a - dependence of the Born cross section is taken into account with the help of the 
-l—distribution, which is defined as usual by 
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where g(x) is an arbitrary test function. 

The integration of subtraction term yields 
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with ACirij and the function £ defined in (|4.2.25| ). The mass and infrared singularities 
of Q4.2.34p have a similar form as in ( |4.2.24| ). Note that the -(--distribution acts on the 
matrix element, the observable O (see ( [4.2.2|) ), and the momenta pi and p a . Here and in 
the following case, the momenta of the four-particle phase-space are implicitly understood 
to be functions of x a with x a 



Xi 



Initial-state emitter with final-state spectator 



Borrowing the previous definitions of the momenta of the four-particle phase space ([4.2.27|) 
and of the variables Xi a and Zi a ( |4.2.28| ), an appropriate subtraction term is defined by 
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As in the previous case, the mappings 
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with k±p a = k±pi = result in the splitting of the phase-space integral 
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The variable 0fc ± is the solid angle between p a and in the rest frame of p a + Pi and the 
masses A and m a can be omitted in Q4.2.38 ). The ^-function results from k\ < 0. 

The integration of the subtraction term over Zi a and xi a for fixed p a and p« provides 
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with A <C m a and the function £ defined in ( [4.2.25 ). The x a -independent singularities are 
similar to the previous integrated subtraction term. Extra x a -dependent mass singularities 
appear in (f4.2.40|) and are proportional to the splitting function P(x a ). 

While in Ref. |5IJ] the momentum p a has been fixed for the :r; a -integration, here we fix 
the momentum p a . Although the actual form of the subtraction terms are different, the 
results for integrated observables are the same. One can reproduce the result of Ref. 
with 
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Initial-state emitter with initial-state spectator 

Finally, the subtraction terms for initial-state emitter and spectator are considered. The 
subtraction term reads 
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with the momenta of the four-particle phase space 

Pa = XabPa, Pb = Pb 

and the variables 
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Since the momenta of the four-particle phase space have to fulfil momentum conservation, 
the momenta Pi of the final-state particles are also modified: 
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(4.2.45) 



with following auxiliary momenta 
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Equation (4.2.45) is a proper Lorentz transformation of the momenta Pi into pi in the 
massless limit .0 

The mapping from the five-particle phase space into the four-particle phase space reads 



q i-> x ab ,v ab ,k ± : q 



V ab - (1 -X ab ~ Vab) ~ 

PaPb 



p b + (l- x ab - v ab )p a + k± (4.2.47) 



with k±p a = k±pb = 0. 

Hence, the splitting of the phase-space integral takes the form 



n 



d 4 pi d A q 
(2tt) 3 (2tt) 3 



S(p 2 - m 2 )6{p lfl )5{q 2 - \ 2 )6(q )5 {4 \pa + p b - EiPi ~ q) 



f o Axa&XbXY J ^5(pj - m 2 )6(p lfi )6^ (p a +Pb - EiPi) J d$ab (4.2.48) 



3 The mapping of the final-state particles coincides with Ref. |p8fl for vanishing fermion masses. 
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with 



ab 



PaPb 
87T 2 J(> 



2 77 



d(p k 



l—x„ 



2tt 



dx ab / dv ab 5(x a - x ab )5(l - x b ) 



xO\v ab - 



ZpaPb 

(1 -x ab -v ab ) 2 m 2 a + X 2 ' 



(4.2.49) 



[1—Xab- V ab )2p a p b 

The variable <pk ± is the solid angle between p a and k± in the rest frame of p a + p b . The 
masses of the final-state particles of the four-particle phase space, 

m 2 - [(1 - x 2 ab )ml + A 2 ] 

, if 2 + Kir / 



m 2 



x 



+ 4 



Kpi\ 2 _ A {Kp l + Kp l )Kp l 



K 2 



(K 2 + KK)K 2 



(4.2.50) 



can be replaced in ( |4.2.48|) by the masses rrii of the final-state particles of the five-particle 
phase space after neglecting m a and A. The ^-function is due to the requirement k\ < 0. 
Finally, the integrated subtraction term is given: 



d^abVab 



VaQl 

8vr 2 



x 



6{l-x b ) 
£(>?,m 2 a,2papb) + r - 



6(1 - x a ) + [P(x a )l 



•^Bornfe, VaiPb) ~ + { ^6(1 - X a ) + [1 - X a 



hi 
1 



ZPaPb 



mi 



M { vl rn {Pa 7 -Va 7 Pb) 



+ 0{\) + 0{m a ) 



(4.2.51) 



with A <C m a and the function C defined in ( 4.2.25 ). The momentum p a is implicitly 
understood to be a function of x a with x a = x ab . 



4.2.4 Remarks to four-fermion production 

Since the virtual corrections are calculated in DPA, the integrated subtraction terms that 
are added to the virtual radiative corrections have to be treated in the right way to achieve 
the correct cancellations of the singularities. However, it is not possible to apply the DPA 
to the whole integrated subtraction terms because the x a -dependent parts are evaluated 
for reduced CM energies ^x a s, which are below the W-pair production threshold for 
small x a , where the DPA is not possible. Therefore, the integrated subtraction terms have 
to be divided into two parts, one including all x a -dependent terms and a second which is 
evaluated in DPA including all infrared and mass singularities, except for the x a -dependent 
mass singularities resulting from ISR. Here, the following decomposition is used: 

/ ^ k V lk ~ y° PA + V°r hel \ (4.2.52) 

irDPA K-iKke 2 QiQk x , \ f /, \/V\2 2 O ~on~on\ 

V ik = S{l-xi)8{l-x 2 )C{\ ,m i ,2p < p k ) 

.(4) .... 2 



X 



M^UPT.^pT) (4-2.53) 
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with i, k = 1, . . . , 6, and the on-shell momenta defined in ( |B.1| ). The function £ is defined 
in ( }4.2.25|) . Different choices of the decomposition (|4.2.52|) of the subtraction terms differ 
only in non-doubly-resonant contributions. 

The singularities of (}4.2.53|) cancel with the singularities of the virtual corrections, while 
x a -dependent mass singularities remain in V^~ shc11 from ISR. Moreover, in the integral over 
the five-particle phase space of Q4.2.2|) , the cuts applied to the bremsstrahlung process 



and to the individual subtraction terms are different, because the observable O depends 
on the five-particle phase space for the bremsstrahlung process, but on the four-particle 
phase space in the case of the subtraction terms. Hence, one should take care of the 
correct implementation of cuts in the phase-space generators, if the subtraction terms are 
involved in the calculation. 



4.3 Numerical results 

For the numerical discussion, the radiative corrections are evaluated as described in the 
previous sections, i.e. the virtual corrections are calculated in DPA, and the complete 
bremsstrahlung process e + e~ — > 4/7 is included for the real corrections. Although polar- 
ized cross sections can be calculated in our approach, only unpolarized cross sections are 
considered in this thesis. 

The fixed-width scheme and following input parameters are used for the numerical 
discussion: 



a 
M w 

r z 



1/137.0359895, 
80.26 GeV, 
2.4971 GeV, 



G„ = 1.16639 x 10- 5 GeV" 



w 



2.08174 GeV, 



M H 
M z 



300 GeV, 
91.1884 GeV, 



(4.3.1) 



except for Section [4.3.3| where the input is given explicitly. If not stated otherwise, the 
weak mixing angle is defined by c w = M w /M z , s^, = 1 — c^, and the following fermion 
masses are used: 



m e = 0.51099906 MeV, m M 
m u = 47 MeV, m c 
m d = 47 MeV, m s 



105.658389 MeV, m T = 1.7771 GeV, 
1.55 GeV, m t = 165.26 GeV, 

150 MeV, m b = 4.7GeV, 



(4.3.2) 



where the light quark masses are adjusted in such a way that the experimentally measured 
hadronic vacuum polarization is reproduced. 



The input parameters coincide with those of Ref. ||50|| , except for the additional Z- 
boson width. The finite Z width is required, since the momenta of the bremsstrahlung 
process are generated for off-shell W bosons. For vanishing Z-boson width, the matrix 
element of the bremsstrahlung process becomes singular if the real photon has the energy 
E 1 = -y/i — Mz- This singularity is due to diagrams where the incoming electron and 
positron annihilate into a virtual Z boson after radiation of a bremsstrahlung photon. 
The virtual corrections are evaluated for vanishing Z-boson width. 
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If not stated otherwise, we use the parameterization, where the Fermi constant 
and the fine structure constant a are related by 

Gn = X — . (4.3.3) 

The symbol Ar denotes the radiative corrections to the muon decay. The cross section 
in G^ parameterization can be obtained from the results in the on-shell renormalization 
scheme in the following way: 



da"" =da£ rn (l + ^) 



with 



d 0. _ da Born eG M _ d <irt + d <cal _ . * i-loo P I A o A\ 

The quantities crg orn , cr" irt , and a^eai denote the tree-level cross section, the virtual, and 
the real corrections, respectively, calculated in the one-shell renormalization scheme with 
the free parameters a(0), Mw, M%, Mh, and m^. The value of Ar can be calculated from 
( [4.3.3| ) and the one-loop corrections to the muon decay yield Ar 1_1 °° p = 0.0373994 for the 



set of input parameters (|4.3.1|) and ( |4.3.2|) . 
4.3.1 Total cross section and angular distributions 



In order to compare our results with the results of Ref. [p0|1 , the doubly- resonant elec- 
troweak 0{ot) radiative corrections are calculated for the leptonic process e + e~ — > z/ M/ u + r _ P r 
and the LEP2 energy 184 GeV. In Ref. ||50|| , the DPA is applied to both the matrix el- 
ement and the four-fermion phase space, except for the Breit-Wigner propagators. The 
integrations of the Breit-Wigner propagators over the invariant masses k\ are extended 
to the full range (— oo, +00), resulting in 

1 7T 

dfc2 t — (4 3 5) 

, ± 4-M^ + iM w r w M w r w ' v ' ' ; 



In contrast to Ref. [5D|, the radiative corrections in our approach are evaluated as 
described in the previous sections with the exact bremsstrahlung process, and the virtual 
corrections are taken in DPA. Whereas in Ref. ||50|| , the momenta are generated for on- 



shell W bosons, we use the exact off-shell phase space for the calculation of the radiative 



corrections. The DPA-Born cross section is evaluated in the same way as done in Ref. |50 
in particular with on-shell phase space, in order to have a common normalization of the 
relative corrections. Moreover, the real photon is recombined with a charged final-state 
fermion if their invariant mass is smaller than all other invariant masses 771(7, /)> where 7 
denotes the photon and / an initial- or final-state charged fermion. 

CM-energy dependence of the total cross section 

In Fig. [4.5| the total cross section and the corresponding relative correction factor 5 are 



given as a function of the CM energy. The radiative corrections are large and negative, 
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Figure 4.5: CM-energy dependence of the total cross section and the relative corrections 
for the process e + e~ — > v^^ + t^v t 



especially close to the W-pair threshold. This effect is due to real-photon ISR, which ef- 
fectively reduces the available energy of the W-pair production subprocess, combined with 
the fact, that near the W-pair threshold the cross section is rapidly decreasing with de- 
creasing energy. The large corrections result from diagrams of the bremsstrahlung process 
where the real photon is emitted from an initial-state electron or positron. 

The DPA-Born cross section agrees very well with Fig. 8 of Ref. |5(J. The curve for 
the cross section including 0(a) corrections agrees for high energies, but differs for low 
energies. For instance, the cross section at 165 GeV is about 7% larger than the cross 

5(J. Above 190 GeV the relative corrections agree very 



section taken from Fig. 8 of Ref 
well with the results of Fig. 9 in Ref 



Note that the dominating corrections resulting 
from ISR are evaluated in Ref. |K| in DPA with on-shell phase space, i.e. k± = M^, but 
calculated in our approach with the full off-shell kinematics. This could account for the 
deviation to Ref. fl50l. 



Production-angle distribution 



The production-angle distribution is shown in Fig. P~6[ This distribution is, in particu- 
lar, important in order to get more strict bound on the non-standard triple-gauge-boson 
couplings. The radiative corrections are negative, and increase in size with decreasing 
production angles. The origin of the distortion of the distribution can be traced back to 
hard initial-state photonic corrections. Hard-photon emission boosts the CM system of 
the W bosons, causing a migration of events from regions with large cross section in the 
CM system (e.g. forward direction) to regions with small cross section in the laboratory 
system (e.g. backward direction). 
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Figure 4.6: Production-angle distribution for the CM energy s/s = 184 GeV and the 
process e + e~ — > v^h + t~u t 



Our result based on the DPA-Born cross section agrees very well with Fig. 10 of 
Ref. [j50| . The relative corrections are about one per cent larger than the results of Ref. [S^ 
for small scattering angles and agree for large scattering angles . Note that the total cross 
section of Ref. |50| is about one per cent smaller than our result at y/s = 184 GeV. 



Decay-angle distribution 

The decay-angle distribution of the final-state fi + with respect to the W + boson in the 
laboratory frame can be found in Fig. 4/7. As in the case of the production-angle distri- 
bution, the large corrections are mainly due to hard-photon boost effects. The radiative 
corrections are negative and large for small decay angles. 
The differential cross sections of Fig. |47 



agree very well with Fig. 16 of Ref. |50 
Deviations of about one per cent are visible from the relative radiative corrections shown 
in Fig. 17 of Ref. 



50 



4.3.2 Invariant-mass distribution 

For the definition of realistic observables the momentum of collinear photons have to be 
recombined with the momentum of the nearest fermion, except for a muon in the final 

Otherwise logarithms of the form ln(m^/s) remain 



state, as discussed in Section 4.1.2 



from FSR in the invariant-mass distributions. These logarithms are not calculable in our 
approach, which is explained in the following: we use the subtraction method, where all 
mass singularities, which appear for inclusive observables, are transferred from the real 
part to the virtual part of the radiative corrections with the help of subtraction terms 
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Figure 4.7: Decay- angle distribution of the fi + with respect to the W + in the laboratory 
frame for the CM energy y/s = 184 GeV and the process e + e~ — > z/ M /i + r~z/ T 



and the fermion masses are neglected everywhere, except for the mass singularities. For 
observables, which are not inclusive in the collinear region, like invariant-mass distribu- 
tions without recombination cuts, the phase-space integration over the difference of the 
bremsstrahlung cross section and the corresponding subtraction terms in ( 4.2.2|) becomes 
divergent for vanishing fermion masses. Therefore, a comparison with the invariant-mass 
distribution of Ref. |50 is not possible, since the calculation is performed without recom- 
bination cuts. 

The following photon recombination procedure is used. All photons within a cone of 
5° around the beams are discarded. If the event is within the recombination cuts, the 
photon is recombined with the nearest final-state charged fermion, more precisely, with 
the fermion that has the smallest invariant mass with the photon. Finally, all events are 
discarded if a charged final-state fermion is within a cone of 10° around the beams. The 
recombination cuts read: 

recomb a: m(f, 7) < 5 GeV, 

recomb b: m(f, 7) < 25 GeV, 

where m(f, 7) denotes the invariant mass of a final-state charged fermion and the photon. 

For the invariant-mass distributions, the Born cross section is calculated with the 
complete matrix element and off-shell kinematics. The radiative corrections are evaluated 
as in the previous section. 

The invariant-mass distributions are shown in Figs. 4.8 , 4.9 , |T0, and |4.11| for two 
recombination cuts and two CM energies: the LEP2 energy 184 GeV and a possible linear 
collider energy 500 GeV. The results include the leptonic process e + e~ — > v il ii r T^v T ^ the 
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Figure 4.8: Invariant-mass distributions for different photon-recombination cuts and en- 
ergy y/s = 184 GeV 



semi-leptonic process e + e~ — > ud/i"^, and the hadronic process e + e~ — > udsc. Note 
that the radiative corrections of all leptonic processes are equivalent in DPA if no mass 
singularities remain from FSR. The same is true for the semi-leptonic and hadronic process 
classes. 



As discussed in Section [4.3.1| , the large negative radiative corrections for 184 GeV orig- 
inate from hard-photon ISR, which reduces the available CM energy for the W-pair pro- 
duction subprocess, combined with the fact that the total cross section of the W-pair 
production subprocess is steeply decreasing for decreasing CM energy near threshold. 
Since the linear-collider energy 500 GeV is far away from threshold and the cross section 
of the W-pair production is slowly decreasing for increasing energies for \/s < 500 GeV, 
the relative corrections are small and positive. 

Apart from the large reduction of the cross section for 184 GeV due to ISR, a distortion 
of the Breit-Wigner line shape is caused by FSR, which can be explained as follows: if 
the event is outside the recombination cuts, the invariant W-boson masses are defined by 
the momenta of the corresponding fermion-antifermion pairs: 

M\ = k\ = (ki + k 2 ) 2 , M 2 _ = k 2 _ = (k 3 + h) 2 , (4.3.6) 

where k%, . . . , k± are the momenta of the final-state fermions. For FSR, one of the two 
resonant W bosons decays not only into a fermion-antifermion pair, but also in a brems- 
strahlung photon. Thus, the corresponding W-boson propagator, 

depends in addition on the photon momentum q and leads to a shift of the Breit-Wigner 
line shape to smaller invariant masses. These effects are especially large without recombi- 
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Figure 4.9: Invariant-mass distributions for different photon-recombination cuts and en- 
ergy v/i = 500 GeV 



nation cuts since mass singularities remain from FSR ||50|| . With recombination cuts, the 
mass singularities disappear and are effectively replaced by logarithms depending on these 
cuts. A small recombination cut leads to large distortions of the Breit-Wigner resonance 
shape, while large recombination cuts yield relatively small effects. 

If the photon is inside the recombination cuts, it is combined with the nearest charged 
fermion and hence is included in the invariant mass of the corresponding resonant W 
boson: 



Ml = (k + + qf 
The resonance of the W-boson propagator, 

p(4) = 



or 



M_ 



[k 2 ± - M 2 } 



(4.3.8) 



(4.3.9) 



which do not depend on the photon momentum, is located at larger invariant masses, i.e. 
M± > Mw, where M± is defined in (|4.3.8Q . This leads to a shift of the Breit-Wigner line 
shape to larger invariant masses. 

These effects all visible in Figs. O, O, |4.10|, and |4.11|. For M± < M w , the relative 



corrections increase for decreasing invariant masses due to propagators defined in ( 4.3.7|) 
and events, which are outside the recombination cuts. This effect is large for small recombi- 
nation cuts, as expected from the previous discussion. On the other hand, for M± > M w 
the relative corrections increase for increasing invariant mass owing to the propagators 



defined in ( |4.3.9|) and events, where the photon is recombined with the nearest fermion. 
This effect is visible for the large recombination cut b. Both effects result in a shift and a 
broadening of the Breit-Wigner line shape. 
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y/s/ GeV 


189 


500 


2000 


10000 


a/fb 


1284.9(5) 


570.4(3) 


83.4(6) 


6.52(1) 



Table 4.1: Comparison with Table 3 of Ref. |22[: Total cross section of the process e + e — > 
udjtt"^ with the photon mass m 7 = 10~ 6 GeV 

Note that the radiative corrections from FSR include the squares of the charges of the 
final-state fermions, which emit the bremsstrahlung photons. Thus, the FSR corrections 
are proportional to Qf± = 1 for invariant-mass distributions of W-bosons, which decay into 
leptons, and proportional to Ql + Ql = 5/9, if the W-boson decays into a quark-antiquark 
pair. 

4.3.3 Test of the subtraction method 



In Table 3 of Ref. p2| , the sum of the soft-photon cross section with E y < to and the 
hard-photon cross section with photon energy > uj of the process e + e~ — > ud/i~z/ M 7 
is given for several values of the separation cut u>, in order to show the independence on 
the parameter u. Since the soft-photon cross section depends on the photon mass m 1 = 
10~ 6 GeV, it is not a physical observable. However, it is a good test for the implementation 
of the subtraction method. 

In Ref. | ]22[| , the fixed- width scheme and the following input parameters have been 
used: M w = 80.23 GeV, T w = 2.085 GeV, M z = 91.1888 GeV, T z = 2.4974 GeV, m e = 
0.51099906 MeV, m u = 5MeV, m d = lOMeV, and si = 0.22591. All couplings have been 
parameterized by «(Mw) = 1/128.07, except for the couplings of the bremsstrahlung 
photon which have been parameterized by a = 1/137.0359895. 



Our results for these input parameters are shown in Table |4.1| and agree well within 



the statistical error, apart from the CM energy lOTeV, where the numerical integration 
is most complicated. 
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Figure 4.10: Invariant-mass distributions for different photon- recombination cuts, different 
processes, and energy yfs = 184 GeV 
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Appendix 



A Calculation of real and virtual corrections 
A.l Useful definitions 

In the main text, we have already used the following short-hand expressions, 



R 1 4M W j_ ■ 

/3 W = \ 1 + ie, 



s 



p = \ 1 



4M 2 





= \J A(Mw, S13S24, S14S23) - ie, 
IP _ p I 

A M = ^ (A.1) 



where 



X(x, y, z) = x 2 + y 2 + z 2 - 2xy - 2xz - 2yz. (A.2) 
The evaluation of one-loop n-point functions naturally leads to the usual dilogarithm, 

U 2 (z) = - — ln(l-t), |arc(l - z)\ < tt, (A.3) 



/o t 

and its analytically continued form 

Ci 2 (x,y) = Li 2 (l - xy) + [ln(xy) - ln(a;) - \n(y)] ln(l - xy), 

I arc(x)|, I arc(y)| < tt. (A. 4) 

A.2 Calculation of the real bremsstrahlung integrals 

In this appendix we describe a general method for evaluating the bremsstrahlung 3- and 
4-point integrals defined in (|3.2.15|) . We make use of the generalized Feynman-parameter 
representation ||]| 

-r(n)/ dx 1 ---dx n , (A.5) 



n? =1 N> "'Jo (ZLiXfr 



103 



where the real variables aj > are arbitrary, but not all equal to zero simultaneously. The 
sum Yli=i N'iXi must be non-zero over the entire integration domain. 

We first consider the 3-point integral Co- We only need the difference between the 
general IR-finite integral and the corresponding IR-divergent on-shell integral, 



C {pi,P2,0,m 1 ,m 2 ) - C (pi,P2, A, \/pJ, yjpj 



d 3 



q 



7iq I (2piq + pj- mj)(2p 2 q +p%- mf) (2piq) {2p 2 q) 



Vq 2 +A 2 ' 



(A.6) 



which is UV-finite and Lorentz- invariant. Extracting the signs <7j of p^ via definition 
( |jj.2.29 ), and assuming for the moment o"jRe(p^ — mf) > 0, the Feynman-parameter rep- 
resentation (|A.5| ) can be applied to each term in the integrand of (|A.6| ). The integration 
over d 3 q can be carried out and yields 



Co(p 1 ,p2,0,m 1 ,m 2 ) - C (pi,p 2 , A, \JpJ, \fp^ 2 



(A.7) 



dxidx 2 5[ 1 — aiXi 
i=i 



ln(A 2 ) + In 



i=l 



r 2 



2 In 



i=i 



Putting for instance ai = and a 2 = 1, and using 

1 1 



da; 



X + Z\ x + z 2 J 



ln(l + 



zx) 



£i 2 (zi, z) - £i 2 (z 2 , z) 



(A.8) 



where | arcf^i^, z, z zi y2 )\ < vr, the remaining one- dimensional integration and the analyti- 
cal continuations to arbitrary complex (pf — mf) are straightforward. 
Next we consider the IR-finite 4-point function 



T> (p 1 ,p 2 ,p3,0,'m 1 ,m 2 ,m 3 ) 
"d 3 q 



7rq (2p 1 q + p\ - m\ ){2p 2 q +p| - rn%)(2p 3 q + pi - ra| 



(A.9) 



go=|q| 



For <7i Re(pf — mf) > we can proceed as above and find 



T> (p 1 ,P2,P3, 0, mi, m 2 , m 3 ) = <J\ a 2 cr 3 / d^id^dxs 

Jo 



5(l - ^aixA 
^ i=i ' 



2 2\ 

m 4 )Xi 



i=i 



(A.10) 

Again, the two-dimensional integration over the Feynman parameters and the analytical 
continuations in (pf — m\) are straightforward. 

Finally, we inspect the IR-divergent 4-point function 



D„( Pl ,p 2 ,P3,A, v ^,m 2 , V ^) = / 



1 



d 3 q 

irq {2piq) (2p 2 q +p\~ rn$) (2p 3 g) 
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?0 : 



Vq 2 + A2 



• (A.n; 



Instead of applying the Feynman-parameter representation directly, it is more convenient 
to extract the infrared singularity by subtracting and adding an IR-divergent 3-point 
function containing the same IR structure. Since the difference between the 4- and 3- 
point functions is IR-fmite, we can regularize the IR divergence in these two integrals in 
a more convenient way Therefore, we write 



£>o(pi,P2,P3, A, \fp\, m 2 , 
1 



lim 

n 2 — >p 2 



2 -C (pi,P3,0,m, slpl) + — 



^o(pi,P2,P3,0, 



Co[pi,P3, A, 



m, m 2 , 



m 



(A.12) 



i.e. we regularize the IR divergence in the 4-point function by the off-shellness p\ — m 2 ^ 0. 
Both the 4-point function as well as the difference of 3-point functions can be treated 
as above, yielding straightforward Feynman-parameter integrals according to (|A.7| ) and 
( |A.10| ), respectively. The limit m 2 — > p\ can easily be taken after the integrations have 
been performed. 



A. 3 Explicit results for scalar integrals 
Loop integrals in double-pole approximation 

The following one-loop integrals are required for the calculation of the non-factorizable 
corrections in Sections 3.3.2 and 3.3.3 . The 2-point and 3-point functions of the (mm), 
(mf), and (im) corrections read in DPA: 



k\-M 2 
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(A.13) 
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(A.14) 



In 



m e M w 



In 



+ In 



#4 



A 2 



+ In 



Mi 



w 



+ 



7T 

6~ 



(A.15) 



In the (mm') correction the following combination of 3-point functions appears in the 
strict DPA: 



*4= M w 



C (fc + , -fc_, 0, M, M) - [C (fc+, -A;-, A, M w , M w ) 

^— |-£i 2 ^^,a;w) +£i 2 (-^ =1 ) ^w) + Li 2 (l - a;^) + vr 2 + ln 2 (-x w ) 



s(3- 



+ 21i1 (aM^) ln(:Ew) ~ 27riln(1 _ 



(A.16) 



105 



In order to include the full off-shell Coulomb singularity, one has to add 
2vri f(3 + A M - (3\ 2vri , ( K+ + K- + s(3 w A M 



— In 



s/3 \p + A M + pj sp w V 

to the r.h.s. of ( |P6| ). 

For the (if), (mf ), and (ff') corrections we need the following 4-point functions: 



(A.17) 



D (p+, fc 2 , A, m e , M, m 2 ) 
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(A.18) 



(A.19) 
(A.20) 
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(A.21) 



A)(2) = D (-k 3 , k+, k 2 , A, m 3 , M, m 2 ) 
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(A.22) 



^(3) = D (-k 3 , -k-,k 2 , A, m 3 , M, m 2 ) = D (2) 



^ A'_ ,?Tl2<->?723,Si3^S24 



(A.23) 



D (4) = J D (-A; 3 ,-A ; _,fc + ,0,m 3 ,M,M) = D (i; 
Bremsstrahlung integrals in double-pole approximation 



(A.24) 



In the (mm') interference corrections the following combination of 3-point functions ap- 
pears: 



C (k+, k-,0, M, M*) - [C (fc+, fc_, A, M w , Af w ) 

C o (A; + ,-A;_,0,M,M) - \C Q (k+, A, M w , M w ) 
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2m , r K+ + K*_x w 
m — 

s/3 w LiAM w (l — x 



(A.25) 
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For the (mf') and (ff') interference corrections the following 4-point functions are re- 
quired: 



P o (0) = -D (h, k + - k 3 , k 2 - k 3 , 0, M_, M+, 0) 
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X> (2) = £>o(&3, k + , k 2 , A, m 3 , M, m 2 ) 
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iXm 2 (s 13 + s 23 ) 
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X> (3) = V (k3,k_,k2,\,m3,M*,m 2 ) 
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P (4) = Z> (fc 3 , 0, m 3 , M*, M) 



-D (4) 
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KI(s 13 + s 23 )-K + M 
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In 1 + 
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-In 1 + 



x w M^ 

Sl3 + s 23< 

(A.30) 



We note that the logarithmic terms on the r.h.s. of ( |A.26| ) yield purely imaginary contri- 
butions if («i 3 + s 23 ) > -M^,x w , (s 2 3 + S24) > -M^i w , and k w is imaginary. These 
conditions are fulfilled on resonance. Off resonance, this is no longer true near the bound- 
ary of phase space. But since these conditions are only violated in a fraction of phase 
space of order \k± — M^|/M^, which is irrelevant in DPA (compare the discussion of the 
relevance of the Landau singularities in Section |3.4.3|) , it would even be allowed to replace 
Vq(0), which is real, by — Re{-D (0)}, as it was done in Refs. |58|, j39| . 



B Four-fermion momenta for on-shell W bosons 



In order to define a DPA for the virtual corrections of four-fermion production, one has to 
project the four-fermion momenta for off-shell intermediate W bosons to the ones, where 
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the W bosons are on shell. Therefore, eight independent parameters have to be chosen, 
which determine the four-fermion phase space uniquely and include the invariant masses of 
the resonant W bosons k±. One possible choice is to fix the directions of the particles W + , 
fi, and fz in the CM frame of the initial-state e + e~ pair while taking the limit k± — > M^. 
This results in the on-shell projection: 

V£ i_on _ ^+ , 

2 ' + IkJ 2 



M 2 M 2 

l.on _ 7, W Lon _ jLon _ r on r on _ j W (.on _ .on _ i on / R i \ 

A/i — Ajl . , An — rt i rui ■ A q — A Q . , A ,i — Ay A/q I JJ.X J 

1 1 2k 1 k^' 2 + 1,3 i 2k 3 k??' 4 - 3 v ; 

with fc + = k\ + A;2 and k- = h% + ^4. The momenta of the initial-state e ± are denoted by 
p± and the momenta of the final-state fermions by hi, i — 1, . . . , 4. The on-shell momenta 
are marked by the superscript 'on'. 

C Construction of phase-space generators 

The scattering amplitude of a given process includes usually a lot of different propagators 
corresponding to the different diagrams. The different propagators peak in different phase- 
space regions and, therefore, slow down the convergence of the numerical integration or 
cause numerical instabilities. If the differential cross section varies too much, the result 
and, in particular, the statistical error of the Monte Carlo integration can not be trusted. 
As described in Section |2.2|, the multichannel approach is applied to obtain accurate 



predictions for the processes e + e~ — > 4/ ( + 7), where each channel is responsible for a 
diagram or a class of diagrams with the same propagators. A channel is defined by 
the mappings from random numbers to the momenta of the event. For each diagram, a 
phase-space generator is constructed in such a way that the corresponding propagators are 
smoothed by the Jacobian of the specific mapping. To do so, the inverse of the Jacobian 
of the channel has to simulate the cross section in those regions of phase space where the 
cross section becomes large. The construction of the phase-space generators is explained 



in the following (see in addition Refs. p6| , |64| , |65|1). 



C.l Smoothing propagators depending on time-like momenta 

Diagrams usually involve several propagators depending on time-like or space-like mo- 
menta. In this section, diagrams which depend only on propagators with time-like mo- 
menta are considered. It is appropriate to decompose the integral over the phase space 
into several integrals over the invariants of the propagators and the remaining integrals 
over the angles in appropriate rest frames. Each invariant is mapped to a random number 
r such that the integral takes the following form: 

ds = / 2 (C.l) 

Jo g s (s(r),m 2 , is, s min , s max ) 

with < r < 1 and the density g s defined by 

1 dh(r,m 2 ,v, s min , s max ) 



g s (s(r),m 2 , is, s min , s max ) dr 



(C.2) 
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The mapping of the random number r into the invariant s 



s(r) = h(r, in 2 , u, s 



mm ; ° max ; 



(C.3) 



has to be chosen such that the density g s simulates the behaviour of the integrand in 
the region where the integrand is large. This is called importance sampling, because 
more events are sampled in the important integration region of s in which the integrand 
becomes large. Thus, it is a reasonable approach to require that the density has to include 
the inverse of the propagator. 

The mappings belonging to the different propagator types read: 



Propagator with vanishing width [^(J : 



a oc l/\s — m 
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ln(s max - in 2 ) - ln(s min - m^Jj (s - m^j 

for v = 1. 

Breit-Wigner propagator: a oc l/[(s — My) 2 + M v r v ] 
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(y 2 - yi )[(s - M 2 y + M 2 T 2 y 



S min/max 

. M v r v , 



(C.6) 



(C.7) 



These mappings are applicable not only for the Mandelstam variable s, but also for the 
absolute value of the Mandelstam variables \t\ and \u\. 

The parameter v can be tuned to optimise the Monte Carlo integration and should be 
chosen > 1. The naive expectation v = 2 is not necessarily the best choice, because the 
propagator poles of the differential cross section are partly cancelled in the collinear limit. 
Note that the mappings ( |C.4[ ) and ( |C.5| ) are undefined for v > 1 and s min < in 2 . 

Since all fermion masses are neglected in our calculations for four-fermion production, 
m can be omitted in ( |C.4j ) and (|C.5|) . However, I found it convenient to use a small 
negative mass parameter in 2 = —a to avoid numerical problems for s m i n = 0. 
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Isotropic decay into two particles 



Phase-space generators of diagrams, which include only propagators depending on time- 
like momenta, can be composed by isotopic particle decays. Therefore, the isotopic decay 
is described in the following. 

A particle with momentum p decays into two particles with momenta k\ and fc 2 . The 
momenta of the initial-state particle and the masses of the final-state particles are fixed. 
The polar angle <ft* and azimuthal angle 9* in the rest frame of the decaying particle are 
suitable integration variables: 



d 4 fc 1 d 4 fc 2 5(fc 1 2 - ml)e{k 10 )S{k 2 2 - m 2 2 )8{k 20 )5 (4) {p -k x - k 2 ) 



— / d(p dcos6> 

8jr Jo J-i 



(C. 



with the function A defined in ( A.2 ). For the numerical integration using the Monte 
Carlo technique, the angles have to be expressed by random numbers r\ and r 2 with 
< 7"i, r 2 < 1: 

1 r l 

d$ d {p,ml,ml) = ——. — = — 2T / dr i dr 2 (C.9) 



(C.IO) 



g d (p 2 \m\,m\) Jo 
with = 27T7"i, cos#* = 2r 2 — 1, and the density 

2p 2 



I 2 2 2 

g d (p ,m 1 ,m 2 



\2 (p 2 , m 2 , vp^)tx 

Since the laboratory frame usually does not coincide with the rest frame of the decaying 
particle, the proper Lorentz transformation is introduced. The Lorentz transformation of 
momentum k into the rest frame of the particle with momentum p is defined by 

k' = B( Po ,p)k (C.ll) 

with the explicit form (see e.g. Ref. |Q ) 

bk 



K = lh + bk, k' = k + b- 



bfcn 



(C.12) 



1-7 

where b = — p/m, 7 = p /m and m = \fp 1 . The inverse Lorentz transform is obtained by 
substituting p into —p. For instance, if the momentum k is identical to p, the particle is 
transformed into its rest frame and vice versa: 

p* = B(po, p)p, V = B(po, -p)p* (C.13) 

with p* = (m,0,0,0) T . 

Since the decay is isotropic, the orientation of the coordinate system can be arbitrarily 
chosen and the momenta of the outgoing particles read 

I p 2 +fc 2_ fc 2 ^ 

2^ 




h =B(p ,-p)K(<f>*,cos6*) 





A^(p 2 ,fc2,fc|) 

V 2^ J 



k 2 =p-k 1 



(C.14) 
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Figure 4.12: Diagram with two resonant W bosons of the process e + e — > z/ M /i + r v T 



with the explicit rotation 



K{(j)*,cos6*) 



( 







o o\ 



cos0* sin0* 
— sin0* cos0* 
\0 I) 



(l \ 

cos 6* sin#* 
10 
\0 - sin6* cos9* J 



(C.15) 



Example 



As an example, a phase-space generator for the diagram of Fig. ^4.12| is constructed. For 
both intermediate particles, photon or Z boson, only one generator is required, since the 
corresponding propagator is fixed by the momenta of the initial-state particles and, hence, 
no mapping is needed. The phase-space integral can be decomposed of the integrals over 
the invariants of the intermediate W bosons, the W-pair production, and the two W 
decays: 



d$ 



p 

dk 2 
o Jo 



X 



d$ d (k + ,k 2 ,k 2 2 ) J d$ d (k_,klkl) 



(C.16) 



with p = p + + p_, k + = k\ + k 2 , and k. 
propagators are determined by 



ka + kd. The invariant masses of the W-boson 



k\ = /i(r l5 M^-iMwr w ,2,0,p 2 ), 
k 2 _ = h(r 2 , Ml - iM w r w , 2, 0, (yp2 - Jk 

with the function h defined in (U.6). The total density reads 

gtat = 9s(k 2 + , - iM w r w , 2, 0,p 2 )g s (k 2 _, M 2 , - iM w T w , 2, 0, ( 
xg d (p 2 , k 2 + , k 2 _)g d (k 2 + , kl, kj)g d (k 2 _,kl k\). 

including the two W-boson propagators. 



(C.17) 
(C.18) 



- \k 



(C.19) 
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C.2 The t-channel diagram of a 2 — > 2 particle process 



So far, only phase-space generators for diagrams with propagators depending on time-like 
momenta are considered. The simplest diagram including virtual particles with space-like 
momenta is the t-channel diagram of a 2 — > 2 particle scattering process which is outlined 
in the following. 

Two particles with momenta p + and p_ transform into two particles with momenta k\ 
and k 2 , where the momenta of the initial-state particles and the masses of the final-state 
particles are fixed. Furthermore, the diagram includes a propagator depending on the 
Mandelstam variable t = (jp + — ki) 2 which has to be smoothed. 

The calculation is performed in the rest frame of p = p + + p_ with the 3-axis in p + 
direction. The momenta of the external particles take the form 



f p 2 +p±~p% \ 

2^ 





V 



o 



k\ = K{<t>*, cos 9* 



( p 2 +kl-kl \ 







\?(p 2 ,k 2 ,k%) 
\ 2^ 



(C.20) 



and k 2 = p\ + p*_ — k\. Since the invariant mass of the propagator t depends exclusively 
on the azimuthal angle cos 9* , 

2 2 {p 2 + k\ - k\ ) {p 2 + p\ - p 2 _ ) - A 5 {p 2 , k\ , k\ ) A i {p 2 , p\ , pi ) cos 0* 

t - K 1 +p + — , 

the phase-space integral can be converted into an integral over the polar angle <fi* and over 
the absolute value of the Mandelstam variable t: 

J d$ p (p + ,p-,ml,ml) 

d%d 4 k 2 5^\p- Pl -p 2 )5 (h\ - m\) 6(k w )5 (k 2 2 - mj) 9(k 20 ) 

d\t\, (C.22) 



4Aa (p 2 ,p+,p?_) Jo 



2- 



^max 



where the integration boundaries can be calculated from ( U.21 ) with — 1 < cos 6* < 1. 
Like in the previous section, the polar angle <fi* and the invariant t of the propagator are 
determined by 



27T7"i, 



\t\ = h(r 2 , m 2 , u, —t 



max; "mm 



resulting in 



Jd$ p (p +1 p- 1 k%,k%) = J drijT dr, 
with the density 



9p{P^i P+i P— i ti m2 , t m im t max ) 



(C.23) 



(C.24) 
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7T 



A 2 (p ,p + ,p_)g s (-t,m ,v, 



~^min)- (C.25) 
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Figure 4.13: An example for a multi-peripheral diagram contributing to the process 
e + e" — > u e u e fi~fi + 'y 



To obtain the momenta of the final-state particles in the laboratory frame, a rotation 
and a proper Lorentz transformation are performed: 

ki = B(po,-p)K(-4>,coB§)k; (C.26) 

with i — 1,2 and the angles 

arctanf?^) , p+i > * p,o 

' ^ + , cos# = f^- (C.27) 
arctan (|±^J + tt , p +i < |P+I 

with p+ = B(p ,p)p+. 

In order to make the numerical integration more efficient, cuts can already be intro- 
duced in the generation of the event, i.e. an upper limit on t or on cos#*. Since the 
Mandelstam variable t is Lorentz invariant the cut on t is valid in arbitrary frames. An- 
gular cuts have to be transformed into the rest frame of the incoming particles. 

C.3 Phase-space generator for a multi-peripheral diagram 



The multi-peripheral diagram of Fig. |4.13| is investigated as an example, because it in- 



volves all types of propagators. In the following, all fermion masses are neglected and the 
following definitions are used: 

&23 = fa + fa, fa 3 4 = fa + fa + fa, fa 23 4 = fa + fa + fa + fa, 

q+ = p+-fa, q-=p--fa, p = p++p-. (C.28) 



The infrared and collinear singularities are excluded by a lower limit on fc| 3 > k 2Z 
The phase-space integral 

2 2 2 

d$ = r d/4 r dfa 234 r dk 2 12U f d%{ P+ ,p_,k 2 1234 ,ki) 



linn 



x J d$ p (p + , g_, kl k 2 2U ) J d$ d (fa u , k 2 23 , k\) J d<S> d {k 23 , kl kl) (C.29) 
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is decomposed of two scattering processes and two particle decays. The intermediate 
particles of this decomposition are the virtual particles and an additional fictitious particle 
with momentum ^1234- The invariant masses of the external particles of the scattering 
processes and particle decays have to be determined first: 

k h = h ( r i,0,u u kl 3min} p 2 ), (C.30) 

k% 3 4 = h ( r 2,0, "2, kj 3 ,p 2 ), (C.31) 

fci234 = h(r 3 ,0,0,k 2 2U ,p 2 ). (C.32) 

As a second step, the momenta of the final-state particles are calculated: 



P+ + P- -> ^1234 + h 



The initial-state electron and positron transforms into the final-state u e and a fictitious 
particle with momentum ^1234- The invariant mass of the W-boson propagator is deter- 
mined by 

\q 2 _\ = h(r 4 ,M^-iM w T w ,2,0,p 2 -k 2 12U ), (C.33) 
where the boundaries are taken from ( |C.21 ). The momenta can be calculated according 



to Section IC.2 



p+ + q- — > h + k 



23 i 



The incoming particles are the initial-state positron and the incoming virtual W boson 
which couples to the initial-state electron. Note, that the momentum of the incoming W 
boson is fixed by the momentum of the initial-state electron and the already calculated 
momentum of the final-state v e . The invariant mass of the W-boson propagator reads 

\q\\ = h(r 5 , il4 - iM w r w , 2, 0, (kf 2U - k 2 2U )(k 2 l2U - q 2 _)/k 2 123A ) (C.34) 

with the boundaries calculated from (|C.21|) . In this way, the momenta of the final-state 
v e and the virtual photon are determined. 



^234 ~~ * ^23 + ^4 



The virtual photon decays isotropically into the final-state /i + and the virtual \i according 
to Section |C.1| . 



^23 -> k 2 + k 3 



Finally, the virtual fi decays isotropically into the final-state \x and the final-state 7. 

The total density reads 

fl'tot — 9s{k 23) 0, i>i, k 23tmin ,p 2 )g s (k 2U , 0, v<2, k 2 3l p 2 )g s {k\ 23A) 0, 0, kl u ,p 2 ) 
xg p (p 2 , 0, 0, ql, - iM w T w , 2, 0,p 2 - k 2 12U ) 

x 9 P (k 2 1234 , 0, q 2 _,ql, - iM w r w , 2, 0, (kf 234 - kl 34 )(k 2 1234 - ql)/kf 2U ) 
xg d (k 2 2U ,k 2 23 ,0)g d (k 2 23 ,0,0) (C.35) 

and includes all propagators of this diagram. 
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Figure 4.14: The generic phase-space generators symbolized by the topologies of the 
corresponding diagrams 

C.4 Remarks on four-fermion production 

In this way, five generic phase-space generators for e + e~ — > 4/ and ten for e + e~ — > 4/7 are 
constructed. The topologies corresponding to the phase-space generators for the processes 
e + e~ — > 4/ are shown in the first row of Fig. |4.14j and the topologies corresponding to 
e + e~ — > 4/7 can be found in the second and third row of Fig. |4.14| . For each generic 
generator the order of the external particles and the mappings for each propagator can be 
chosen. In particular, flat mappings of a subset of propagators correspond to topologies 
where these propagators are contracted. With the help of these generators all propaga- 
tors of the cross section are smoothed. Special phase-space generators for interference 
contributions to the cross section are not constructed. 

For the process e + e~ — > e + e~e + e~7 with 1008 diagrams 928 different channels are 
included in the calculation. The difference is due to diagrams where the initial-state e + e~ 
pair couples to a virtual photon or Z boson. Since the invariant masses of the virtual photon 
and Z-boson are fixed by the CM energy, only one phase-space generator is required for 
both diagrams. For instance, the diagram of Fig. |4.12| corresponds to the first topology of 
Fig. |4.14j, where the mapping of the propagator which couples to both incoming particles 



is flat, i.e. z/j = 0. 

Phase-space generators for diagrams with a quartic-gauge-boson vertex are also in- 
cluded. This diagrams have one propagator less than the other diagrams. Therefore, one 
propagator in Fig. |4.14j is omitted and the mapping of the invariant of this propagator is 



flat. 

The order of the particle decays, the 2 — > 2 particle processes, and the determination 
of the invariants influence the convergence behaviour of the numerical integration, since 
the integration domains of the single mappings depend on the already fixed invariants 
and momenta. In particular, the invariant masses of decaying particles in a decay chain 
are not independent of each other. To improve the numerical stability, the invariants 
of the virtual massless particles that couple to external particles should be calculated 
first, applying already the cuts in the generation of the momenta. In the here discussed 
Monte Carlo program, invariant-mass cuts are used for the propagators of the virtual 
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particles with time-like momenta. For the virtual particles with space-like momenta only 
one angular cut is applied to the first 2 — > 2 particle process. Note, that the angular cut 
does not influence the boundaries of the integrations over the invariants of the virtual or 
fictitious particles with time-like momenta, which are calculated at the beginning. 

For the numerical integration of the subtraction terms, additional phase-space gener- 
ators are constructed, which take into account the propagators of the Born-cross section 
included in the subtraction terms. Since the Born cross section depends on the four-particle 
phase space, the four-particle momenta have to be generated first with the phase-space 
generators for the processes e + e~ — > 4f. These momenta are then mapped into the five- 
particle phase space according to the different mappings defined Section |4.2.3| . 

Besides momentum conservation and the requirements on the momenta of external 
particles, i.e. kf = mf and k i0 > 0, the phase-space volume is a very helpful variable to 
test the phase-space generators. For massless external particles, the phase-space volume 
can be calculated analytically as follows 



J n ^h^eik^ip-T.Uh) = y T(n ) Tin _ iy ( c - 36 ) 



where n is the number of final-state particles (see e.g. Ref. [66]). 
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